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Measurement of Instantaneous Vector Air 
Velocity by Hot-Wire Methods’ 


CARL E. PEARSONT 
Harvard Unwersity 


GENERAL SUMMARY 


One of the major problems encountered in improving the de- 
sign of jet-engine compressors has been to examine experimentally 
the airflow present in any given compressor stage. Because of 
the rapid and large fluctuations in velocity, the usual types of 
instrumentation will give a (doubtful) measure of the average 
yelocity alone. However, a hot-wire probe has been developed 
which permits measurement of instantaneous, as well as average, 


direction and magnitude under the above conditions of airflow 


Details of the construction, theory, and use of this hot-wire 
instrument are presented. Each wire used need be calibrated 
oily once; although the characteristics of a given wire change 
fapidly and somewhat erratically, a correction method has been 
such changes 


devised which automatically compensates for 


{ This correction method, based on the hypothesis that the changes 


(in characteristics are due to a combination of oxidation and plas 
tic flow of the hot-wire material, has been firmly validated by ex 


tensive experimental work 

In the measurement of instantaneous air velocity, it might ap 
pear that, because of the predominant effect of the heat capacity 
of the wire, the instantaneous velocity should be proportional 
to the rate of change of voltage across the hot-wire. This is in- 
deed the case, but to prove this rigorously (for large velocity 
changes), as well as to determine the constant of proportionality, 
fequires considerable algebraical manipulation 
ways in which the present theory may be 
in general, the agreement is excellent 


There are 
checked experimentally ; 
A typical example is given in which the instantaneous air ve 
locity behind an axial compressor is measured in three different 


many 


directions 
Received May 25, 1951 Revised and received October 6, 
1951. 
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The instrument has so far been examined only for Mach Num 
bers less than 0.7. For M < 0.3, an accuracy of better than 2 
per cent in average and 10 per cent in instantaneous velocities is 
claimed; at J = 0.7, the error in instantaneous velocity meas 
urement may be 25 per cent 

It appears that hot-wire instruments of the present type will 
extend considerably the range of experimental work possible in 
compressors and other mechanisms. The instrument as here dis 
cussed has now been used satisfactorily over a period of about 1 
year in the analysis of the airflow in the axial and centrifugal com 


pressors, both in and out of pulsation. 


(I) INTRODUCTION 


I A THIN WIRE IS HEATED by an electric current and 
inserted in an air stream, its rate of heat loss will 
depend on the velocity of the air and on the temper 
ature difference between wire and air. Since about 
1910, this fact has been used to a small extent in the 
measurement of mean velocities of steady airflows and 
to a considerable extent in turbulence research. 

The usual turbulence circuit is as shown in Fig. 1. 
The battery B supplies the hot-wire // with heating 
current through a series resistance FR sufficiently large 
that the current is almost independent of hot-wire 


resistance. The voltage across the wire is measured 
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Fic. 2. Probe construction. 


by the potentiometer P; the mean air velocity is then 
a function of this potentiometer voltage, the air tem- 
perature, and the characteristics of the wire. 

Since the air stream is, in general, not perfectly uni- 
form, the voltage across the wire will possess an a.c. 
component that is amplified by A and presented on an 
oscilloscope O. The magnitude of this voltage affords 
a measure of air-stream turbulence; to simplify cal- 
culations, it is customary to add a special circuit to the 
amplifier which will automatically compensate for the 
way in which the hot-wire response is modified by its 
own heat capacity. This compensating network is 
valid only for small velocity fluctuations and can there- 
fore not be used in the present investigation. 

The object of the hot-wire work at Harvard has been 
to develop an instrument that will give an accurate 
reading of the average and instantaneous vector ve- 
locity of airflow behind the impeller of a compressor. 
For this purpose, it was necessary to modify the cir- 
cuit of Fig. 1 in several ways. In addition, it became 
necessary to find a correction method applicable to 
the irregular changes in operating level to which all 
hot-wire materials are subject. 

Various minor troubles, such as dirt collection on the 
upstream surface of the hot-wire, existence of variable 
contact potentials, hunting in electronic power supplies, 
etc., were overcome by filtering, careful soldering of 
contacts, use of battery power supplies alone, etc. 

Both theory and instrumentation went through sev- 
eral changes as the various difficulties arose; the final 
forms of each will now be discussed in detail. 


(II) PROBE CONSTRUCTION AND MEASURING 
APPARATUS 


To avoid the temperature dependent voltage drop 
in the leads between probe and measuring apparatus, 
separate current and potential leads are used (since the 
potential leads do not carry current, the voltage drop 
across them is always zero). To measure accurately 
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the average air temperature at the probe, a small cop. 
per-iron thermocouple is built into the probe just be. 
low the hot-wire. 

The arrangement is as in Fig. 2. A total of gj 
wires is required, each of which must be insulated fron 
the others and from the metal tubing containing them, 
The most satisfactory hot-wire material is tungsten; jt 
is copper plated at its ends and soldered onto the sup. 
ports. During soldering (carried out under a micro. 
scope), the supports are pressed slightly together so a 
to hold the wire straight after release. It is important 
that the wire be well soldered at each end, withoy 
bends. 

The apparatus is shown schematically in Fig. 3. The 
amplifier is direct-coupled so as to avoid phase shift in 
coupling condensers. D is a differentiating circuit, 
The resistors R’ and R” are noninductive standard dial 
boxes. The switching arrangement shown allows the 
potentiometer P to measure the voltage across either 
R” or the wire itself; in practice, R’ is adjusted so as 
to keep constant the voltage across R” (and therefore 
also the hot-wire current). 

In the sequel it will be shown that the resistance of 
the hot-wire must be measured frequently; to do this 
a low-voltage dry cell is substituted for the battery B 
so that the resulting current will be too small to raise 
the wire temperature appreciably above air tempera- 
ture. The oscilloscope used is a Du Mont 304H (pos. 
sessing an internal d.c. amplifier); it is synchronized 
to the compressor by means of a tiny rotating magnet: 
coil system. 


(III) PROPERTIES OF A HOT-WIRE IN STEADY FLOW 


In an inviscid incompressible steady flow of fluid, the 
equation governing the conduction of heat from a hot: 
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MEASUREMENT BY 
wire 1S 
Cc (dT /dt) = Ky V?T (1) 
where 
c = fluid heat capacity per unit volume 
K,) = thermal conductivity (assumed constant) 
rT = temperature of fluid 


dT/dt = material derivative of T = (O7/Ox)u, + 
(OT /Oy)u, + (O7/0z)u, in steady flow 
where u,, etc., are the components of 
velocity at the point in question 


The boundary conditions are that 7 = 7, at in- 
finity and 7 = 7), at the wire (it is assumed that the 
wire has constant temperature throughout its mass). 
It is difficult to find an analytic solution of this bound- 
ary value problem, but one property of the solution is 
easily obtained——viz., that the rate of heat flow H/ from 
the wire has the form 


H = (T. — T,)f(V) (2) 


where {(’) is a quantity depending only on the velocity 
V at infinity. Suppose that a solution of Eq. (1) were 
known for a certain value of 7), 7T,; then, since 
Eq. (1) is linear, it is clear that, if 7, — 
bled, the new temperature gradient would be found by 
merely doubling the value of the gradient of the pre- 
vious solution at every point, and therefore the total 
rate of heat loss (proportional to the gradient at the 
A gener- 
It may be 


T, were dou- 


surface of the wire) would also be doubled. 
alization of this reasoning leads to Eq. (2). 
noted that, since air is neither incompressible nor in- 
viscid, this result is indicative only. 

The accuracy of the empirical formulas for f(V) ar- 
rived at experimentally by various authors is not suffi- 
cient for our purposes, and it is necessary to find /(V) 
for each wire used by calibration in a wind tunnel. 

The hot-wire material used is unplated tungsten wire 
with a measured diameter of 0.00885 mm. and a mass 
of 0.234 mg. per 200-mm. length. The resistivity at 
room temperature is of the order of 8.2 K 10~* ohm- 
cm., the deviation from the handbook value of 5.5 X 
10-* ohm-cm. being due to the drawing process. The 
temperature coefficient of resistance is 0.0045 per °C., 
as found experimentally by use of a constant-tempera- 
ture distilled water bath. At the normal operating 
current of 100 ma., the current density is 163,000 amp. 
per sq.cm.; because of this high density, it was neces- 
sary to check the validity of Ohm’s law, and this was 
done at various temperatures by means of the water 
bath. No significant deviation was observed. 


EFFECT OF TEMPERATURE VARIATION 
THROUGHOUT MASS OF WIRE 


(IV) 


Consider first the temperature variation in a direction 
normal to the axis of the wire. The rate of heat loss 


irom the wire surface is about 0.02 watt; using a sur- 


HOT-WIRE 


METHODS 


face area of 0.0003 sq.cm. and a coefficient of heat con- 
ductivity for tungsten of 0.4 cal. per sec. cm. °C., this 
heat-loss rate requires a temperature gradient of the 
order of 40°C. per cm. Since the wire diameter is only 
8.85 & 107-4 em., it is reasonable to assume that the 
wire temperature is substantially constant in a direc- 
tion perpendicular to the axis. This assumption is 
lent further strength by the fact that if any portion of 
the wire were at a lower temperature, the resistivity of 
this portion would be correspondingly lower and the cur- 
rent distribution would automatically readjust itself, 
the result being to increase the temperature of this por- 
tion. 

Consider now the effect of temperature variation 
along the length of the wire. It will be assumed that 
the two ends of the tungsten wire are at air temper- 
ature 7,; this assumption has been verified experi- 
mentally by measurement of thermal e.m.f.’s. 

Assume the heat loss per unit length to be given by 


H = Ko(T» — T,)¢(M) (3) 


where Ky is the thermal conductivity of air; 7, — T,, 
as before, is the difference between wire and air tem- 
M is the rate of mass flow of air; and ¢is an 
empirical function of M. Let the wire have length L, 
area A, and resistivity p, at air temperature 7,. The 
resistance at 7, is then Ry = p4(L/A). Under oper- 
ating conditions, the wire temperature 7), will vary 
along the wire; the resistivity at a particular value of 


T,, will be given by 
p = pall + a7, 


where a is the temperature coefficient of resistance. 
The wire resistance will then be given by 


perature ; 


T,) | 


pac a a 
R=R,+ (T. — T,4) dx (4) 
A Al 
Since the wire is in equilibrium, the rate at which 
energy enters an elemental wire portion of length dx 
must equal the rate at which it leaves; if / is the cur- 
rent and K is the thermal conductivity of tungsten, 
the resultant equation reads 


Ip. 


r 
A ad) + 


[1 + a(Ty, al 


1. 
AK = 


= (T,, — T4)Kod (5) 
dx* 


which may be solved subject to the end conditions 
T, = T, atx = Oand x = L to give 
I*ps 


i _— Z, _ a \ 
' bKvA —_ alps 


( sinh r(x — L) — sinh =) 
1+ (6) 


sinh rL 


where 


r? = (@Ky)/AK) — (al*p,/A*K) (7) 
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Numerical values show that rL = 8 at 150 ft. per sec.; 
since the e~”’ portions of cosh rL and sinh rL are then 

ae ni P ‘. rL . : —_ 
negligible in comparison with the e’” portions, Eqs. 


(4) and (6) give 
I*p, a =*( 2 ) 
— Ry = | — (S) 
tiaes A*Kr? rL 


which is a cubic in r. To find @¢, the value of r as ob- 
tained from Eq. (8) is merely substituted in Eq. (7). 
The energy lost through the ends of the wire may be 


dT, 2 
= I*p. : 
-_ r 


found by 


2KA (9) 


dx L 

For a typical wire at 150 ft. per sec., the value of dA» 
is 4.36 X 10~* watts per cm. °C., and r = 61.6 per em. 
Then Eq. (11) gives the rate of heat loss through both 
ends as about 23 per cent of a total energy delivery of 
0.01919 watt. 

(V) CORRECTIONS FOR CHANGE IN CONDITIONS 

Even after elimination of contact resistance, provision 
of better air filtering to reduce the effect of dirt collec- 
tion, implementation of the various circuit changes 
previously discussed, and consideration of local air tem- 
perature, the wire potential was still subject to erratic 
and sometimes frequent changes in potential, usually 
Discriminatory ex- 
other 


in the direction of an increase. 
periments showed that the cause of this effect 
than the easily corrected for dirt collection—was a 
combination of oxidation and plastic flow; it there- 
fore became desirable to develop a rapid correction 
technique to account for these processes. 

Essentially, the method of correction for oxidation 
and plastic flow is to designate a certain set of conditions 
as standard and to then calculate what a given reading 
would be under these conditions. 


Combining Eqs. (10) and (11) so as to eliminate ¢(.)/) and using Eq. (14) to put ZL. Z; in terms of Ryo/ R41, some 


2 | = ey 2 2 | 
4. S = 
L Rs I, ry Yolo 
2 |+ Rel (7)'s| 
rol» Rael rE , 


algebraic manipulation yields 


[s(7!) 


where 


‘ Koop ore ry i/e Rye > 
S = LL + a@(TMa2 — Ta1)] \G (16) 


Ko q il 


For generality, 7; and /» have until now been assumed 
different. Almost always, however, /; = /2; using 
the fact that S > 1 and 1/71; & 1/rele, Eq. (15) then 


reduces to (approximately) 


AERONAUTICAL 


ea)l -s(7), 


SCIENCES 


FEBRUARY, 


1952 


In the following, standard conditions, to which al 
observations must be reduced, will be denoted by th 
subscript “‘1.’’ Thus, at standard conditions, the air 
temperature will be denoted by 7°,:, the wire resistang 
at 74; will be Ry, and the wire will have length L, ay 


diameter D;. Then Eq. (8) will read 


TVpnaR,; 


R, ™ Rat = s- me ( = — (1{ 
OKA, 7 aly? pai r1L, 


For the same mass flow rate ./, but under a new nop 
standard condition that will be indicated by subscrip; 
‘2”” (e.g., due to plastic flow, wire length is now | 


etc.), Eq. (8) will read 
L2?p42aR 42 (1 2 1) 
poKwAs = als*pye rols ‘ 


Ke = Ras = 


where the heat loss rate ¢ has been replaced by po bx 
cause of the modifying influence of the change in diam 
eter from D, to D, (due to oxidation and plastic flow 
on the heat-loss rate. The factor / is close to unity 
for otherwise the wire resistance would change much 
more drastically than measurements indicate. 

If the wire was subject to plastic flow without oxida 
tion, there would be no change in volume; because oj 


possible oxidation, 
L,D;; = qgL2.D-" 12 


where g is also close to unity. The air-temperatur 


resistances are connected by 
L, D,? is i ‘ 
Ra = Ray yp, U + a(T 42 — Ti) \td 


Using Eq. (12), Eq. (13) may be written 


Ly _ (F")( l ) 
L, WRal\ell + Tes — Tad) 7) ™ 


R, —= Rs (Rai Re) (17 


If the conditions “subscript 2"’ and ‘“‘standard”’ are 
considerably different, Eq. (15) may be used directly 
However, 


(g alone is unknown; may be set = 1) 


Eq. (17) is almost always sufficiently accurate. 
are measured, 


e.m.f.’s rather than resistance 


fk, = Rh, ete.; then Eq. (17) becomes 


Since 


define 
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By = E2(Eai/E a2) (18) 


(Of course, /;, etc., are less than the measured e.m.f.’s 
by the voltage drop across the previously mentioned 
2 cm. of needle supports. However, for rapidity of 
calculation, the e.m.f.’s of Eq. (18) may be replaced by 
the actually measured e.m.f.’s without appreciable 
error.) 

An example may clarify the method. Suppose that a 
potential reading of 0.1972 volt is obtained at 150 ft. 
per sec. at an air temperature of 21.7°C. The resistance 
at this air temperature is 1.673 ohms. Standard con- 
ditions for this wire are chosen as 7,4; = 20°C., Ry: = 
1.650 ohms. J; = J, = 0.1 amp.; then E, is given by 
Eq. (18) 


E, = 0.1972 (0.1650/0.1673) = 0.1945 


In order to check the correction theory described 
above, a large number of tests were carried out under 
widely varying conditions. Two typical curves are 
shown in Figs. 4 and 5. Both curves include measure- 
ments taken over a 3-day period, the wire being occa- 
sionally cleaned and washed. It will be noticed that 
most of the experimental results lie above the corrected 
values; this is because the initial resistance of the wire 
was used as a basis for the corrected curve and the re- 
sistance gradually increased with time. In general, 
the resistance varies sufficiently slowly that it need be 
measured only every 10 min. or so; its value at inter- 
mediate times may be obtained by linear interpolation 
(taking into account air temperature, of course). 


(VI) VARIABLE AIRFLOW 


As before, subscripts ‘‘1”’ and ‘‘2”’ will refer to stand- 
ard and test conditions, respectively. Let c be the 
specific heat of tungsten (joules per Gm. per °C.) and 
o its density. Then, for an elemental portion of the 
hot-wire, of length dx, the rate of accumulation of energy 
is 


Substituting Eqs. (22) into Eq. (21) gives 


A 2° oC OR,’ 


re) 
= T,? (Roo -o R,') = 2K A> . (T 1,20 + 
paca Of Ox 


AL SCIENCES—FEBRUARY, 1952 


(A» dx) ac(OT,,2/Ot) 


where O77,2/0f is the time rate of change of temperatyy, 
of this portion. Equating this rate of accumulation t) 
the difference between rates of input and loss of energy 

8 


gives 
OT 1p p2 Ol we 
. Peay jue oe re a 
or ” Ao . "ioe? 


Koupd [M(t)] (Tvs sind T 42) (19 


A2oc 


where the notation ¢[1/(t) ] means that @ depends on the 
mass flow rate M, which in turn depends only on the 
time /—1.e., the scale of the fluctuation is large in com. 
parison with the wire size. 

Integration of Eq. (19) with respect to x yields 


d Le 
Asoc / T 2 dy = I,? Ro = i 


Or 
a Le od Fr 
KA» =] — Kwopdod / (Te — T42) dx (20 
Ox 0 J0 
Using Ty2 — T42 = (p2 — par)/paca, Eq. (20) becomes 
Asoc OR» ‘ - Olu po 
- Ts?Ro +KA>» = 
pao Ot Ox tI, 
Koz y A» 
— (Re — Ry) (21 
4 P22 
The term 
OT we | 
KA, ] 
Ox Io 


represents the instantaneous heat-loss rate through the 
ends of the wire. Now, the quantities R2, 7,2, and ¢ 
may each be written as the sum of a steady component 
(denoted by subscript ‘O0’’) and a fluctuating com- 
ponent (denoted by a prime) 


R, as Ro» + R.’ 
Tp2 = T 20 + loa (22 
i) = go + ¢’ ) 


Koo : 0 f A» 
Tw2") |x os —_ ==) - (Roo + R,’ — R42) (23 
Pa2e 


Taking mean values by integrating Eq. (23) in time from 0 to 7 and dividing by 7 gives, after using 


” OR,’ 
dt = R.'(r) = R,'(0) 
0 ot 


and allowing 7 to approac: 


0 
( 2*Ro — 2KA2 — (Tu) ]x=0 — 
) ie a Ox ( Ie=o 


where the averaging process carried out for the last term on the right-hand side is permissible only if R,’ is ex- 


KospooA 2 ( 
Paz 


Rep _ Rae) (24 


tremely small in comparison with Rx» — R42; for the present, this is assumed. 
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Substitution of Eq. (24) into Eq. (23) gives 
‘ ra) =e A»*oC OR,’ K 2 bayA 2 
—2KA, (Tw2") | = = f R.’ + T.?R»' 
$’ Ox paca Ot p42Q 5) 
= . (25 
(KoepAoe Pa2a) (Rep + R» -_ R42) 
With no greater error than about 2 per cent, the denominator may be simplified by neglecting R,’ with respect to 


Rw — Raz; then, using Eq. (24) again gives 


» OT 2" A220¢ OR,’ KeopooA 2 
ry —2KA; ai i R,’ + I,°R’ 
Kod Ot 1. os pasa Ol P42 »6 
= . —_ Fe : (20) 
Koodo 12*Ro» om Ts*Rao (2/ropL2) 
It is now necessary to examine the numerical magnitudes of the various quantities involved. From experi- 


ment, the denominator has a magnitude of the order of 0.015 watt, p = 1, Ro’ = 4 X 10-* ohms, Koodo & 4 X 
10-* watts per cm. °C., A2/p42 = 0.07 cm. per ohm, and a = 0.0045 per °C. Then, 





(Ko2pooA 2/ps2a) Ro’ = 0.0003 watt 


which is sufficiently small in comparison with the denominator that it may be omitted from Eq. (26), for the vari- 
ations in Koo’ / Kod of interest in this report are of the order of 0.2. Writing 


a. 1d. 
( Tw2 |. ae ~ pe 
Ox ap42 Ox 


where, as before, pz = px + p»’, it may be reasonably expected that the order of magnitude of 


1 oO 1 ” (5) _ 2R2’ 
aps2 Ox os x=0 tg? AP 42 tz, 2 m 2 — aR 


therefore the first term in the numerator has a magnitude of about 10~* and may consequently be neglected. The 


Il? 


result is 
Kod’ (A»s*oc ‘pa2a) (dR,’ dt) 27) 
ne (2 
Kood To*Ro — I2?Ra2 (2/reol.2) 
Kno’ and Kod) may now be put in terms of the standard resistance R,. From Eq. (8), 
oo + ¢’ — [Ra/(Ri — R4i)] [1 — (2/nL1)] + 1 
do [Rai (Rio —s Rai) ] (1 — (Zz rio1-1) | + | 
which may be solved for R, to give 
Rai [1 — (2/rnL1)] 
R, =R 28 
a ae (: 2 )+ eh 4 _ Ra (: 2 )] ” 
Rw» — Ray i0h01 do Rw —_ Ray Pipl} 
Substituting for ¢’/@») from Eq. (27), Eq. (28) may be written 
2 
l —_ 
i+ Bn te Ce rill . 
~ - _ wed | 2 aa = Rai(2 atl 
riol-1 apaz dt LRai — Rai(2/riol) (29) 
To?Ro mast To*Rao (2 rn 2) 


Because of the magnitudes of the various terms involved, the use of only the first term on the Taylor expansion of 
Eq. (29) will not result in substantial error. 
Further, since r depends on the square root of ¢, the difference between 2/r,L, and 2/rj/, is small for the ve- 


locity changes under consideration and may be neglected 


Rw = Ray A 2" OC Rw ~~ Ray (2/rioL1) dR,' - 
R; = Rw + = o - (30) 
T,7Ro a 12*Rao (2 roo L2) Para Rai = Rai (2 rol.1) dt 
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Fic. 6. Examination of cosine behavior. 


If C is the total specific heat of the wire in joules per 
°C., then 


c= A»x?acR42 Pa2 (31) 


The reason for the differentiating circuit is now clear. 
Because of the differentiation, the oscilloscope deflec- 
tion is proportional to de’/dt and from a knowledge of 
the constant of proportionality and by use of Eq. (30), 
the instantaneous value of the standard resistance may 
be calculated and, therefore, also the corresponding 
instantaneous mass flow rate. Since the curve of Fig. 
5 is fairly flat, the oscilloscope deflection is approxi- 
mately proportional to the velocity fluctuation, and this 
fact may be used whenever it is necessary to calculate 
approximately a large number of velocity patterns. 

In the previous analysis, it was tacitly assumed that 
the rate of heat loss // of steady flow may be used in 
cases of variable airflow. This is true only if the heat 
gradient field is set up around the wire rapidly in com- 
parison with the velocity fluctuation rates; an order- 
of-magnitude analysis shows this to be probably true 
for the flow behind the axial compressor. However, 
there exists a method of experimentally checking the 
above assumption; such a check will be described 
shortly. 

It may also be noted that the current is not quite 
steady because of the variation in hot-wire resistance; 
however, an easy calculation shows the effect to be 
negligible. 

Because of the relatively minor change in slope of the 
potential-velocity curve as the wire flows plastically 
[see Eq. (12)], a wire may be used indefinitely for the 
measurement of velocity fluctuation without recalibra- 
tion. 

Further refinements in the theory are possible; for 
example, the fact that the thermocouple reads only the 
mean temperature of a steady flow could be considered. 
If it is assumed that the air is adiabatically compressible, 
the temperature fluctuation corresponding to any ve- 
locity fluctuation could be calculated and incorporated 


into the analysis is straightforward and will not be jp, 
cluded here, since it is of little importance for the syp, 
sonic flows of present interest. 

It should be emphasized that, in distinction to th, 
case of turbulence analysis, the preceding analysi 
holds for large velocity fluctuations (of the order of thi 
main-stream velocity). 


(VII) ANGULAR VARIATION 


In the preceding analysis, the axis of the hot-wire 
has been assumed perpendicular to the direction o/ 
the air stream. In practice, of course, this is not the 
case; if the angle between the normal to the wire axis 
and the air-stream direction is denoted by @, then jt 
seems reasonable that the velocity component parallel 
to the wire have little cooling effect, so that the wire re. 
sistance corresponds to a velocity of V’ cos 6. Various 
experimenters have checked this and found good agree. 
ment for @ less than about 70°. 

In the measurement of mean velocities, this approxi 
mation is not sufficiently accurate, and it is necessary 
to calibrate each wire. <A typical test is shown in Fig 
6; the dotted line represents the e.m.f. that the wire 
would indicate on the basis of the cosine hypothesis, 
and the solid line passes through the experimental 
results. It will be noted that, as expected, the wire is 
cooled slightly more than by the amount predicted oa 
the cosine hypothesis. 

Since the hot-wire responds (approximately) only to 
that velocity component perpendicular to its axis, the 
magnitude and direction of a given air stream may be 
determined by merely orienting the hot-wire axis along 
each of three noncoplanar directions. The trigonome- 
try is straightforward—e.g., if the directions are mu- 
tually perpendicular, it is easy to show that the square 
of the magnitude of the velocity is given by one-half 
the sum of the squares of the three velocity components 
perpendicular to the axis. In accurate work, a two- 
step procedure is necessary, as to the approximate 
angles of incidence must first bé calculated on the as- 
sumption that the correction due to the deviation from 
the cosine law is zero; given these angles, the exact cor- 
rection may be obtained from the calibration curve (as 
in Fig. 6). 

If the flow is two-dimensional—as it is approxt- 
mately, behind the rotor of an axial compressor —then 
the hot-wire e.m.f.'s obtained by placing the wire axis 
in each of two mutually perpendicular directions, each 
in the plane of the flow, will give the two-velocity com- 
ponent directly. (The two-dimensionality of the flow 
may be checked by comparing predicted and observed 
e.m.f.’s at other orientations. ) 

The above remarks hold for both the average and the 
instantaneous components of air velocity. If the fluc- 
tuations repeat themselves—as again happens in at 
axial compressor——then the two components of velocity 
fluctuation may be photographed in sequence; other- 
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MEASUREMENT BY 
two simultaneous photographs are necessary, and 


wise, 
may be obtained by the use of two perpendicular 


these 
hot-wires. 


(VIII) MEASUREMENT OF AVERAGE VELOCITY 


Even after compensation for plastic flow and _ air- 
temperature changes, several difficulties remain in the 
measurement of average velocities. The most annoy- 
ing of these are the collection of heat-insulating dirt 
on the wire and changes in the thermal conductivity 
of air (due, for example, to the presence of oil vapor). 
if an absolute measurement of the average velocity at 
any point is required, special precautions must be taken 
for example, the wire may be inserted with its axis 
parallel to the air stream so as to collect no dirt; when 
a velocity reading is required, the wire is then rotated 
into the correct orientation and held there a few sec- 
A typical measurement (carried out with only 
Here the velocities 


onds). 
average care) is Shown in Table 1. 
feet per second) in the right-hand columns are calcu- 
lated from the hot-wire readings at 6 = 0° and 90 
before, 0 is the angle between the air stream direction 
and the normal to the hot-wire axis). These readings 
were taken about '/, in. behind the rotor of a single- 
stage axial compressor. “ The corresponding pressure 
probe reading for 6 = 0° was 169 ft. per sec. 

In some applications, the presence of airflow fluc- 
tuations will affect the average velocity reading. This 
is because of the fact that, in averaging Eq. (23), the 
lat term is given correctly by Eq. (24) only if R»’ is 
extremely small, for, in general, the average of a prod- 
uct is not the product of the averages. This difficulty 
may again be removed by the use of a two-step averag- 


(as 


ing process, in which the average velocities are calcu- 
lated in the ordinary way, the fluctuation velocities 
then being calculated as to be described. Knowing 
these fluctuation velocities, the corrections to average 
velocity may be obtained from Eq. (23). To date, 
this procedure has not been found necessary in com- 
pressor studies. 

The most common use of the hot-wire instrument to 
the present has been the carrying out of traverses of 
axial and tangential velocities across compressor pas- 
sages. A convenient method for doing this is to choose 
ne point in the passage as a reference point and con- 
tnually to refer e.m.f.’s to the e.m.f. experienced at 
this point. Even though the level of the e.m.f. at this 
point will change because of dirt collection, etc., the 
slope of the e.m.f.-velocity curve changes only slowly, 


TABLE | 
Measured Calculated 

6 Velocity Velocity 
0 171 
30 188 191 
45 182 184 
60 158 162 
90 88 
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and relative velocities may therefore be measured ac- 
curately without any special precautions. This “‘slope 
method”’ results in a velocity profile across the pas- 
sage; the height of the profile is then obtained by con- 
tinuity (the total flow being assumed known). 

VELOCITY FLUCTUATIONS 


(IX) MEASUREMENT OF 


As indicated previously, a knowledge of the average 
velocities will permit calculation of the instantaneous 
velocity value from a photograph of the (amplified) 
differentiated voltage the hot-wire. Before 
Eq. (30) can be used, the value of the specific heat C 
must be determined and the amplifications and phase 
shifts of the amplifier and oscilloscope amplifier must 


across 


be known. 

The value of C may be found either by calculation 
from the physical dimensions of the wire or by carrying 
out an experiment in which the wire is placed in a steady 
air stream and an alternating current superposed on the 
Both have been used 
with fair agreement; in the 
method is probably the more accurate. 
used in finding the velocity patterns to be discussed, 
the value of c/aR,4. was 2.71 & 10~° joules per ohm. 

In order to eliminate phase shift in the external am 
plifier and in the oscilloscope amplifier, direct-coupling 
is necessary (i.e., no coupling condensers can be used). 
The differentiation is obtained by use of the well-known 


heating current. methods 


general, calculation 


For the wire 
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capacity-resistance circuit; other circuits may, of 
course, be substituted. The attenuation due to the 
differentiating circuit must be included in the calcula- 
tion of the overall amplification. 

The results of a typical experiment are shown in Fig. 
7; the photographs (tracings being reproduced) were 
taken with the hot-wire behind the rotor of an axial 
compressor. Each minimum in the axial velocity com- 
ponent corresponds to the passage of a blade. 
ing the axial and tangential components, the velocity 
pattern for 6 = 30° may be calculated; 
compared with the experimental pattern for 6 = 


Know- 


the result is 
30°. 
The reasonable agreement lends support to the funda- 
mental assumptions of the analysis. 

As only the alternating component of hot-wire volt- 
age is presented on the oscilloscope screen, hot-wire 
probes without separate potential leads are quite 
satisfactory for use in obtaining velocity fluctua- 
tions. 

There are some cases (e.g., passage of a rotor blade 
wake) in which the instantaneous value of @ falls out- 
side the 70° range; 
to apply and an obvious two-step calculation tech- 


in such cases the cosine law fails 


nique must again be applied. 

As an example of the degree of resolution obtained 
in practice, a typical photograph of axial velocity vari- 
ation is reproduced in Fig. 8. 
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(X) GENERAL REMARKS 


Although considerable experimental technique ang 
experience are necessary for efficient operation of the 
hot-wire, the results obtained from it are quite satis. 
factory. 
theory, the knowledge of the effect of the differentiating 
circuit allows rapid and fruitful qualitative surveys 
of airflows to be made; such surveys are particularly 
useful in the understanding of such phenomena as com. 
pressor surging (pulsation). 


Even without the preceding quantitatiye 


For brevity, many of the experimental checks carried 
out on the theory are not recorded here. For example, 
the effect of a change in hot-wire heating current js 
almost exactly as predicted. The compensation method 
has so far not been found in error even over periods of g 
year. 

The use of hot-wires in high-speed flows (up to 700 ft, 
per sec.) requires no new theory, except for obvious 
effects of compressibility; the details are therefore 
omitted. 

Some further hot-wire possibilities may be suggested, 
If a high-frequency square wave of heating current were 
used and the two portions of the response were sepa- 
rated by a gate circuit, enough evidence would be pres- 
ent for the calculation of both temperature and ve- 
locity; if a more oxidation-resistant material than tung- 
sten could be found, such an instrument might be useful 
in combustion studies. Alternatively, it is possible to 
construct a hot-wire that will give an accurate (but not 
continuous) measurement of the instantaneous velocity 
without the use of our assumption concerning rate of 
heat loss. Suppose a wire to be initially unheated and 
at a certain instant apply to it a fixed heating current. 
Then the rate at which energy accumulates in the wire 
depends only on the instantaneous air velocity (pro- 
therefore the 
initial rate of rise of wire resistance is a measure of the 
instantaneous air velocity. The data for 
use of this instrument would, of course, be obtained by 


vided the air temperature is known); 
necessary 


calibration in steady flow. 
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Air Drag on Cubes at Mach Numbers 
0.5 to 3.5° 


GEORGE E. HANSCHE? anp JOHN S. RINEHARTT 
New Mexico School of Mines 


ABSTRACT 


The air drag on two sizes of steel cubes, '/, and */s in., has been 
carefully determined in the approximate range from Mach 0.5 to 
35. In general, the results are consistent and in substantial 
agreement with similar data obtained by Charters and Thomas 
onspheres. The present data do not yield an explicit value of the 
drag coefficient Cp but rather an average value of the quantity 
(pA for a rotating cube where A is proportional to the area of the 
An “‘average”’ value for A of 1.50d?, where d is the length 
The 


two sizes of cubes give essentially equal values for Cp at cor- 


cube. 
of one edge of the cube, is assumed and Cp is calculated. 


responding Mach Numbers. 


INTRODUCTION 


ennen OF THE AERODYNAMIC BEHAVIOR of a 
missile having a cubical shape is of considerable 
practical importance in connection with certain ord- 
nance problems. The present study was made in order 
to obtain an experimental value for the drag coefficient 
of a cube in flight as a function of Mach Number. The 
air drag on two sizes of steel cubes, !/,; and */s in., has 
been carefully determined in the approximate range 


from Mach 0.5 to 3.5. 


EXPERIMENTAL 


The general experimental procedure used to deter- 
mine the drag coefficients was somewhat similar to that 
used by Charters and Thomas! for spheres with the ex- 
ception that spark photographs of the cubes were not 
taken. 

The '/,- and */s-in. cubes were fired one at a time from 
smooth-bore 30- and 50-cal. guns, respectively. No 
sabots were used in launching them, except at the high- 
est velocities when a small rag wad was placed just be- 
hind them. The */s-in. cubes were fired at several ve- 
locities ranging from about 600 to about 4,000 ft. per 
the '/,-in. cubes were fired over the range from 
The highest velocity in each 


eC. 5 
400 to 3,500 ft. per sec. 
case was the maximum obtainable by this method of 
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firing. Thedatafrom19 */s- and 26 '/4-in. cubes were 
used in the drag calculations. 

A six-spark Aberdeen chronograph, which had six 
grid stations spaced at 25-ft. intervals, was used to ob- 
tain a distance-time trace for each cube. The grids were 
30 in. square and consisted of two sheets of 0.001-in. 
aluminum foil spaced slightly less than !/,; in. apart and 
separated by a sheet of crepe paper. This type of grid 
proved the most satisfactory of a number of types 
tried. The amount that the grids retarded the cubes 
was measured by suspending a grid on a ballistic pendu- 
lum and firing a cube through it. 

The velocity lost by a !/4-in. steel cube in passing 
through the grid was found to be about 3 ft. per sec. 
In calculating the present results, this error was neg- 
lected since it is considerably less than other probable 
errors. 

The masses of all cubes were kept within | per cent 
of one another. One dimension of the '/4-in. cubes 
differed by about 0.003 in. from the other two dimen- 
The uncertainty in area for the steel cubes was 
Recovery of several cubes after 


sions. 
less than 2 per cent. 
firing showed that the loss in mass due to travel through 
the gun barrel was negligible, although the corners of 
the cubes in most cases were rolled over slightly. 


RESULTS 


The calculation of drag, from chronograph readings, 
is subject to certain inaccuracies and approximations 
that have been adequately discussed by Charters and 
Thomas. The manner in which the present data were 
reduced is described below. 

The differential equation of the forward motion of a 
body through air may be written in the form 


m(dv/dt) = —CppAv?/2 (1) 


where m is the mass of the body, p is the density of the 
air, v is the velocity of the body in motion, A is its pres- 


entation area, and C, is the drag coefficient. A solu- 
tion of the above equation is 
—Cp(p/2)(A/m)(x — Xo (2) 


U = Ve 
where v 
locity at Xp. 

It is obvious that the values both of Cp and of A will 
The presen- 


is the velocity at the point x and % is the ve- 


depend upon the orientation of the cube. 
tation area of a cube may vary by a factor of 1.73. The 
shape of this presented area may be flat, wedgelike, or 
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nearly conical. No precautions were taken to prevent 
rotation of the cubes, and examination of the holes 
made by them in the chronograph grids showed that 
they rotated. Hence, the present data were not ex- 
pected to yield explicit values for Cp and A but rather 
an average value of the quantity, C)A, for a rotating 
cube. If it is assumed that A is a sort of average value 
of area which is proportional to the square of one side 
of the cube, then an average value of Cp may be calcu- 
lated by arbitrarily choosing a value for the propor- 
tionality constant. 

In calculating values of CpA, the difference, (x — 
xo), was taken to be 25 ft. Information from three 
consecutive grids was used. 
the cube during its travel from the first to the second 


The average velocities of 


grid and during its travel from the second to the third 
grid were calculated from the chronograph time-dis- 
tance trace. These values were substituted for vw and 2, 
respectively, in Eq. (2), and CpA was calculated. Thus, 
each six-point trace yielded four values of C,,A. 

l 


The values of CpA so obtained for '/4- and */,-in. 


cubes are plotted in Figs. 1 and 2, respectively, against 


formance of Small Spheres from Subsonic 
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the midpoint of the average velocities used in obtaining 
them. 

If it is assumed that the scatter of the points is due ty 
the area and shape variations mentioned above, the 
average values of Cp for the two sizes of cubes should be 
equal. Values of Cp, calculated using one and one-hajf 
times the area of one face of the cube as the “average 
presentation area, are plotted against Mach Number 
Agreement between the two sets of values js 
Use of the data to predict the behavior 


in Fig. 3. 
seen to be good. 
of cubes of sizes other than '/, and */s in. seems justi. 
fied. 

The general form of the C, curve follows closely that 
of the curve for spheres obtained by Charters and 
Thomas.' The transition from subsonic to supersonic 
takes place over a relatively wide range of velocities as 
is the case with spheres. There is no abrupt change as 
the cube reaches the velocity of sound. 
gins to rise at about Mach 0.6 and begins to start down 
at about Mach 1.4. The Charters and Thomas curve 


for a sphere behaves similarly at substantially the same 


The curve be. 


points. 

In passing through the transonic region Cp increases 
from about 0.84 to 
This increase is considerably less than that 


.26 or, percentagewise, about 50 
per cent. 
for spheres, where Cp increases by about 100 per cent. 
The scarcity of points in the subsonic region makes the 
0.84 value little 
weight has been given the lowest '/,-in. cube points 
since their accuracy is somewhat in doubt. It may be 
that the subsonic end of the true curve lies lower than 


somewhat uncertain. Moreover, 


» 


the curve drawn in Fig. 3. 
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On Unsteady Nonlinearized Conical Flow’ 


S. F. BORGT 
U.S. Naval Postgraduate School 


SUMMARY 


It is shown that the equations of unsteady, nonlinear, non- 
yiscous, nonheat-conducting flow may be put in a conical form. 
The problem of a shock striking an infinite wedge is considered. 
An invariance relation is established, and certain properties of 
the reflected shock are examined. These include a discussion 
of a possible form of reciprocal flows; a calculation indicating 
the impossibility of a certain reflected shock shape for air; and, 
based upon an analysis of the characteristics of the flow, a non- 
existence proof for certain combinations of shock strengths and 


wedge angles. 


SYMBOLS 


subscript designating ahead of incident shock 
subscript designating behind incident shock 
subscript designating behind reflected shock 


velocity of sound 


ll 


p = pressure 

x,y, = physica! plane and time coordinates 

u, 1 = flow velocity components in the physical plane 

I = letter designating incident shock 

M = Mach Number 

mM = letter designating Mach reflection shock 

Q = reduced velocity (u — & v — 7) 

R = letter designating reflected shock 

S = entropy 

$s. = slipstream 

7,0,¥ physical coordinates and velocities referred to 

a moving origin 

a,w, w’ = various angular quantities relating to the 
shock reflection pattern 

p = density 

6 = wedge angle 

6. = angle between attached wedge shock and 
horizontal line 

é = x/t | ; . 

conical coordinates 
7 = y/t 5 
g = p/p, strength of incident shock 


INTRODUCTION 


[' WILL BE SHOWN THAT, subject to appropriate 
conditions, the nonviscous nonheat-conducting time- 
dependent equations can be given in terms of new 
variables 

& = x/t, 


The flow equations in this form do not contain the 
length dimension alone, and in the &, 7 field the flow 
is represented by a single map that represents the 
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tion 


flow at all times. This implies that such time-de- 
pendent flow fields grow uniformly with time. If, 
in addition, the boundary does not specifically contain 
a length parameter, then in the sense of Lagerstrom’s 
the field is conical. The more de- 
scriptive phrase for this “unsteady 
flow’’—-is applied to the fields under discussion. 

In order to satisfy the conical conditions for the 
boundary, we will consider the case of a straight shock 
traveling in an infinite field and striking an infinite 


discussion! flow 


case conical 


wedge (Fig. 1). 

The shock reflection phenomena was first considered 
by Mach.’ Since 1943, Von Neumann and his co- 
workers at Princeton have issued various reports*~"! 
describing their work on this problem, both theoretical 
and experimental. 

Additional related experimental work has been done 
by Liepmann'* '* and at the University of Michigan."* 

Von Neumann? considered local conditions only, in 
the neighborhood the shock intersection. By 
applying the Rankine-Hugoniot relations across the 
shocks, with reference to origins at the shock inter- 
sections, it was possible to obtain various angular 


of 
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relations that must hold. Thus he considered the 
regular two-shock configuration (Fig. 2) and the Mach 
reflection, three-shock condition (Fig. 3). 

For the case shown in Fig. 2, with respect to the 
point 0, and for a given shock strength £, the angles 
w and w’ can be computed. 

For the case shown in Fig. 3, again with respect to 
point 0, and for a given shock strength £, it is possible 
to obtain solutions corresponding to the following 
conditions:!® 

(1) The pressure across the slipstream (ss.) must 
be continuous. 

(2) The velocity of the fluid must be tangential to 
ss. for flow through the lower (7) shock and also for 
flow through the upper two (J, R) shocks. 

(3) The kinematic, dynamic, and thermodynamic 
relations across the shocks must be satisfied. 

Shock-tube experiments performed by Smith® and 
the photographs he obtained permit an experimental 
check on the results obtained by Von Neumann and 
that the two-shock 
solution angles check fairly well throughout the range 
of incident shock 


his coworkers. It was found 
three- 
shock solutions for glancing incidence do not check the 
experimental data. 


strengths. However, the 
Bargmann‘ considered the prob- 
lem of the infinitesimal shock hitting an infinitesimal 
concave corner, and by showing that, to the approxima- 
tion considered, the flow is irrotational, he solved the 
resulting Laplace equation subject to the given bound- 
ary conditions. 

Lighthill'® ! 
theoretical 


considered the from the 


point of 


problem 


view and obtained linearized 


solutions to the following: 
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(1) A shock of any strength striking a convey 
corner of angle 7 + 6 with 6 small, Fig. 4. 

(2) A shock of any strength striking head-on , 
very nearly plane corner of angle t + 6 with 6 small, 
Fig. 5. 

In both cases Lighthill shows that, by use of the 
Busemann conical flow transformation, the pressure 
A solution by 
means of complex variable is then obtained. 


field satisfies the Laplace equation. 

More recently, Ting and Ludloff—in a short note jy 
the February, 
AERONAUTICAL 


1951, issue of the JOURNAL OF THE 


ScIENCES—show a_ solution for a 
finite shock hitting a concave wedge of very small 
angle. The wedge face need not necessarily be a 
straight line in this linearized analysis. Their solution 
does not utilize the ~, 7 coordinates and is based on a4 
Possio integral. 

The present investigation considers the problem 
from a different point of view than has heretofore been 
done. The complete nonlinearized form of the equa- 


tions is considered (although viscosity and_ heat 
conduction effects are neglected). 
A brief outline of the present investigation follows: 
(1) The mathematical equations are obtained in 
conical form and an invariance relation is established. 
(2) It will be shown that a two-shock, attached 
shock pattern cannot exist for air. 
(3) The possibility of obtaining reciprocal flows is 
discussed. 
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UNSTEADY NONLINEA 


existence of will be 


The 
shown which leads to (a) a form of discontinuity 


(4) patching curves 
in the solutions and (b) a nonexistence for certain 
flows. 


DEVELOPMENT OF THE EQUATIONS 


We consider the two-dimensional, unsteady non- 
viscous, nonheat-conducting flow, with 


u = U(x, y, t), v = v(x, y, t) 


Physically, the problem under consideration is a 
shock striking a (for the present) symmetrical wedge, 
whose point is at the origin in the x, y plane and which 
extends to infinity in the positive x, y directions, Fig. 1. 
If we assume that u and v can be given in terms of the 
new variables, &, y, defined by 

u = U(x/t, y/t) = 


pie 
v = v(£, n) 


u(&, n) 


then it will be found that our equations are independent 
of x and y, and because our boundary is also inde- 
pendent of any length parameter, it follows that the 
flow pattern can be represented in the &, 7 plane. 
That is, the unsteady pattern will be represented by a 
» plane. The 
initially deduced by Von Neumann 
amply verified experimentally. 


single map in the &, above fact was 


and has been 
The equations governing the flow are conservation 


of mass 


pr + Up, + Upy + pu, + pry = Of 


pi + (pu), + (pv)y = 0\ 0 

also conservation of momentum 
(pu), + (pu?), + (pun), + py = 0 (2) 
(pu), + (puv), + (pv?), + py = 0 (3) 

In addition, an equation of state 
p = p(p, a (4) 


and if we assume that changes of state are adiabatic 
(i.e., do not consider a shock or other dissipation inside 
the field under consideration), then 


S,; + uS, + vS, = 0 (5) 


so that, corresponding to the five unknowns, 4, 2, p, 
p, S, we have five equations 


Now, from Eqs. (4) and (5), it follows that 


1 6) 1 1S 
dp _ ( ) dp 4 (22) r (6) 
dt Op/ s dt OS/ », dt 
and 
1 oO ra) re) re) 1 
as i. ef) 2 
dt Of Ox Oy \Op/ 5 dt 
or 
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o o o 1 /o fe) o 
e+e e+e = P u Pts 3 (7) 
Ot Ox oy a* \Oot Ox Ov 


Therefore, Eq. (1) becomes 


l E 1 Op 4 “4 4 E n ov — 
u v = (8) 
a* Lot Ox Oy ? Lox oy 
Now, if we assume f = /(£, 7) only, then 
oO Od: O On 
ot O£dt dO 
oO Od: O On 
Ox O§& Ox On Ox 
oO Od: O On 
Oy O& Oy On Oy 
and, for — = x/t, 7 = y/t, these become 
re) i. —&O n O 
of toE £On 
o 10 Oo 10 
ox = t OF Oy On 
so that Eq. (8) becomes 
op op dp P) 4 
—£ ~e + a: v + 
| "ne Sy ee me 
| 4 
p + = 0 (8a) 
O£& On 
and Eqs. (2) and (3) become 
, Ou Ou 4 Ou 4 Ou 1 Op 0 (2 
-§ - u v = 2a) 
o£ i On 3 On pose 
_ Ov Ov 4 Ov 4 Oz 1 Op 0 ¢: 
of -U v = (3a) 
df "dn DE dn pdr 


and we see that these equations are independent of 
x, y, and ¢ as assumed. 

We may rewrite Eqs. (8a), (2a), and (3a) in the fol- 
lowing form: 


u— ~Op v—nOp Ou Ov 
a? O€ sg a? On 4 (5. 7 4 = 
(u — &) ~ + (v al = + — = 0 (2b) 
dé On poet 
(u — &) = + e~-9q > + +. 0 (3b) 
og On = pOn 


Substituting 0p/d0& from Eq. (2b) and 0f/0n from 
Eq. (3b) into Eq. (Sb) we obtain, after collecting terms, 


Ou fa? 2) + Ov la? 2) 
a” — (u — &)? a* — (v— n)*] — 
O€é 7 On . 


(9) 


' ey E + | . 
u— t)(v — n) = 
’ On rey: 


Note, the two-dimensional steady-state irrotational 


equation is 








88 JOURNAL OF THE AERONAUTICAL SCIENCES—FEBRUARY, 1952 





Ou la? 2) 4 Ov ts ss es ~) also 
a> — nu” a? — v*| — uv + = as 
Ox oy. Oy Ov . 


Eq. (9) holds everywhere in the flow in the field ABCDR 
of Fig. 6 and will be used later in connection with a 


study of the characteristics of the field. and 
As a first problem, we consider the form the equations shoc 
take in a moving X, Y system of coordinates, where sail 


the origin of the X, Y system is moving along a straight 
line making an angle 6 with the x axis at a constant 





velocity A, Fig. 7. Let u, v be the velocity components 


in x, y, and let U’, V be the components in X, Y. Now 


Then it may be shown that our equations of flow are 


invariant with respect to the new X, Y system. This 





can be proved in several ways. The simplest proof 





is that which utilizes a well-known invariance property 














of a vector equation —namely, the equation is invariant 7 
with respect to a translated origin. We can show easily 
that in the &, 7 plane our equation of continuity becomes 
E-Vp = V-(pu) ther 
and conservation of momentum becomes 
. My? 
(4 — &)-Vu + (1/p)Vp = O 
where the bars represent the corresponding vector 
quantities. J 
This invariance property of the flow equations will met 
be considered shortly in connection with the possibility may 
a al of having reciprocal flow fields. For the present, we sum 
may note that the invariance property is the essential Tab 
- requirement behind the two-shock and_ three-shock lore 


angle solutions obtained by the Princeton group. 
That is, if we take as an origin the triple shock point, 





for example, then the equation in the X, Y plane of 
this point indicates the existence of conical flow with 
- respect to this point. This in turn permits the as- 
sumption that the field is growing uniformly about this 
point as an origin, which is the assumption made in 
& ra computing the three-shock angular relations. 
T 


The following theorem will now be proved: 





First Theorem.—Given a two-shock pattern, with will 

air, Fig. 8, then the reflected shock cannot be attached S 

° to the wedge point. pat 
Fic. 8. Proof: We prove this theorem by assuming an cuss 


attached reflected shock and then proving that a two- 
shock pattern is impossible under this condition. The 











U various regions and quantities referred to are indicated 
es. ara on Fig. 9, in which we use the notation of Lighthill I. I 
R g, v7 In addition, some of the formulas that follow concern- es : 
fG/ON4 iy Haan © ing relations between velocities in regions @) and (0 
M-£ are from Lighthill I. — 
R YY With respect to point 0, we have I 
er i U2 U D 
0: V2 = qa? + — — 3 cos B anc 
cos- 6 cos B floy 
ang 


or 





Fic. 9 V = V qi + (U?/cos? B) — 2q,U ide 
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UNSTEADY 


also 


V u 
WM, = = M;? + 


ay ay 


l 
— 2M, 


2 cos” B ay 


and we will show that 3% < 1, which indicates that a 


shock R is not possible as shown on Fig. 8. That is, 


we will show that 
‘ i U 
NM? + ~~ + 38, — <1 
a," cos B ay 


Now we may easily show that 
q\ = S rf ; (1 ') 
ay 6 ay l 72 


l (6/5) (qi/ a1) 6 AM, 


ay | (ao7/ U?) 51 — (1/4?) 
therefore, we must show that 
j 56 1) ) 
M41 4+ - 
7. 25 cos? B[1 — (1/A7*)] dll — (i /M?)19 


< | 


Mi, M, and cos 6 are interdependent, and the simplest 
method for proving the above inequality is to plot the 
maximum values of the left-hand side for various as- 
sumed values of, M. This has been done in 
Table |, from which it is apparent that 9% < 1; there- 
fore, the attached two-shock pattern is not possible. 


say, 


TABLE | 
Min. Val. 
V V7, cos? B mM, 
2.1 1.0 1.0 0.35 
i) 1.66 0.92 0.55 
10 1.83 O89 0.77 
1.89 0.87 0.86 


The extension of this theorem to gases other than air 
will be considered in another paper. 

Should it be found that the two-shock attached 
pattern exists, then reciprocal flows of the type dis- 


cussed in the next section may be possible. 


RECIPROCAL FLows 


Because of the previously derived invariance property 
of the flow field, we are led to examine the possibility 
of the existance of reciprocal fields. To explain this, 
consider Fig. 10. 

Because the equations of flow are identical for both 
D and 1D’ as origins, it follows that, (1) if shock angle 
and flow deflection at C and A equal shock angle and 
flow deflection at C’ and A’, (2) if angle CDA equals 
angle C’D’A’, and (3) if conditions along AD are the 
Same as along A’D’, then the two flow fields are 


identical. 
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While the reciprocal property seems possible from 
both the physical and mathematical points of view, 
the difficulty in obtaining reciprocal pairs is magnified 
by our ignorance concerning the slipstream shape and 
position. The only known solutions available at this 
time (Lighthill, Ludloff, and Bargmann) can be shown 
to be without reciprocals. Whether this is due to the 
linearized form of the equations used or to some other 
reason is unknown. 

The reciprocal relations, if they exist, would seem 
to be an important property of the shock reflection 
phenomena. It is therefore suggested that this matter 
might bear further investigation, particularly from the 
experimental point of view—especially in the direction 
of checking photographs for possible evidence of this 
phenomena. 
CHARACTERISTICS AND THEIR APPLICATION TO 
THE GIVEN PROBLEM 


THE 


In this section the flow field will be considered for 
the special case in which the shock is attached to the 
wedge point (or equivalently for the regular two-shock 
patterns). Under these conditions, a necessary condi 
tion for the existence of characteristic curves will be 
determined. This necessary condition, when used in 
conjunction with an inequality to be derived in this 
section, will then be examined further, and it will be 
found that a type of transonic singularity or discon- 
tinuity occurs. In addition, it will be shown that for 
certain shocks an attached straight shock solution is 


not possible. 


The Characteristics 


When the Mach Number behind the incident shock 
is greater than unity, for small enough wedge angles, 
straight shocks attached to the wedge point can exist, 
since with respect to this point as an origin the steady- 
state solution is possible (see Fig. 11). The statements 
that follow can also be applied to the case of Fig. 12, 
with the two-shock point considered as the origin. 

Referring to Fig. 11 and confining the discussion to 
the region ABFDA, we see that there is a region behind 
the straight shock which is in uniform irrotational 


motion. 
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Because there are particles in the disturbed flow 
field which have come through the curved shock, it 
follows that there is also a region of rotational flow 
behind the shock." 

Thus, in general, we have two separate pairs of 
regions distinguished by discontinuities, both originat- 
ing at the point B, Fig. 11: 

(1) A region of uniform flow separated from a 
region of nonuniform flow; the nonuniform flow region 
is not necessarily rotational. 

(2) A region of irrotational flow separated from a 
region of rotational flow. 

In mathematical terms, the first pair is characterized 
by the requirement that there be a curve in the (é, 7) 
plane along which uw and v are continuous, but 0u/0é, 


a* — (u — &)? —(u — &)(v — n) 
dé dn 
da da 
1] 0 
0 | 


and this will have a solution if 


(u — &)?+ (wv — yn)? > a? (5) 


In this equation, a is the velocity of sound in the 
uniform flow irrotational region just behind the straight 
portion of the reflected shock. 
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Ou/On, Ov/OE, 
That 1s, there is a patching curve along which th 


v/On are not necessarily continuoys 


derivatives in the velocities may be discontinuoys. 

Similarly, the second pair is characterized by the 
requirement that the vorticity is continuous along 
the corresponding patching curve, but the rate of change 
of vorticity is discontinuous. Physically, this patching 
curve must be a particle path (BE, Fig. 11) through 
the intersection of the straight and curved parts of the 
reflected shock. This will satisfy the known physical 
fact that particles passing through a curved shock 
become rotational and, in addition, satisfies the known 
fact that a particle path may be a characteristic or 
curve of permissible discontinuity of a flow function.” 

It is the first patching curve, however, which is oj 
interest at this time. In the mathematical treatment 
of this curve we use essentially the same method used 
by Murnaghan*! in his analysis of steady-state com 
pressible flows. 

We have our fundamental equation 


ja? — (uw — &)*|] + la? — (v — n)?| 


ry EK 4 =| 
(u — &)(v — n) = () 
On O& 


and along the patching curve 
(Ov/OE) — (Ou/On) = O 


Now, we consider an arbitrary curve C in the (&, 9 
plane given by —& = &a), 7 = n(a), a being any con- 
venient parameter or independent variable, and we 
suppose that the value of u, v and vorticity are pre- 


scribed along this curve. Then, since 


du Ou dé Ou dn . 
da dé da On da a 
dv _ ov dé Ov dn (4 
da 0§ da On da 


A necessary condition for the existence of the velocity 


characteristic is that 


= — 20 — 9) => oe = er | 
0 0 
= 0 
dé dn 
da da 
—] 0 


We may note, too, that the Eq. (5) corresponds to 
the usual steady-state characteristic theory requirement 
that 


g=u+0 >a’ 


or that the flow be supersonic. 
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If the determinant of the coefficients of Eqs. (1) 
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posed on w and v in order that the four equations con 


4) is zero, then the further condition that must be im- necting the derivatives may have a solution is 
dé du 
. 0 0 
da da 
dv dé dy 
0 : = () 
da da da 
0 0 l 0 
=i — 2)" 0 -(s — 2 — 9) a-=— yy — 9) 


from which a sufficient condition for the existence of solutions of the four equations for the four derivatives 


of wu and v along C (when the determinant of these four equations is zero) is 


dé dy 


a? — (v — »)’ 
| ad ee 


From Eq. (6) we can determine the slope d¢y/d& of the patching curve; 


and solving the resulting quadratic equation, we have 


or, if we call Q, the reduced velocity, having components 


u— &v — n), then 
dn (u — é)(v—n) t+ aVQ* — a’ 
= = = ace (S) 
dé a> — (u — &)° 
For use in later applications of the patching curve 


relations, the following Theorem II is proved, Fig. 13: 


“The attached shock at the wedge point will 
always intersect the circle of radius a (velocity 
of sound behind the shock) which is centered at 
(u, v), these being the uniform velocity components 
behind the shock.’ 


Note that the point (#, v) is on the wedge angle line. 
To prove the above theorem, we must show (see Fig. 
13) that un» + 0 > e«, and, since « = 6B + 86, we must 


show that » > 8. This is equivalent to showing that 


sin uw > sin B 


= 1/M., where MJ. is the Mach Number of 
We must therefore show 


But sin yu 
the flow behind the shock. 
that 


1 > M. sin B 


This is done by noting (see Eq. 125, page 45, ref 


that 


tan / 2 l y¥-—1 
= ay, sear 
tan 6 y + 1 \M.? sin? B 2 


which obviously >1 for any finite wedge angle 0. 


erence |S) 


The 


equality holds for 6 = 0. 


dn -2(u — t)(v — yn) + VA(u — &)*( 
| 


: ,, dn du 
— (%@ — &) = (0) (6) 


dada 


for, dividing through by (dé)* 
— yn)? — 4[a* — (v — n)*|[a® — (u — &) 
(u £) 
From the above, 
2 jt Ys = 
Y + l 2a ~ 
where a = J/,* sin® 6 or 
(ay + 2 — a)/(ay + a) 2 1 
and we see that for a = | (that is, @ = 0) the equality 


holds, but only for a < 1 will the inequality be true 
This proves Theorem II. 
Along the patching curve we have 
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that is, if a circle of radius a is drawn centered at 
(u, v), the patching curve must be on or outside this 
circle. Also, Theorem II states that the straight 
attached shock intersects this circle. It follows, 
therefore, that the patching curve must lie in either 
Region A or Region B of Fig. 13. 

One other fact may be noted: For M/, > 1, the 
circle is disposed as shown in Fig. I4a; for IW. < 1, 
it is as shown in Fig. 14b. 

For an infinitesimal shock, the patching curve 
corresponds to the are of circle @) @) of Fig. 4a. 
As the wedge angle increases, the patching curve lies 
in Region A of Fig. 13, except as noted later. 

But, it is a well-known fact (see, for example, ref 
erence 18) that, for stationary attached shocks, as 
the wedge angle increases, the Mach Number J/, 
behind the shock decreases from supersonic values to 
subsonic values. The subsonic values occur only for 
a small range of wedge angles just before separation of 
the shock from the wedge point. The subsonic values 
do, however, occur, and, in addition, they occur for all 
strengths of incident shocks. 

This phenomena leads to the following curious 
behavior of the straight shock in the present problem: 

We consider a shock of given strength (strong enough 
to form attached shocks for a given range of wedge 
angles), and we examine the sequence of events as the 
wedge angle @ is increased. 

For the smaller angles 0, the patching curve lies in 
the Region A, Fig. 13. As @ increases, Region A 
becomes smaller and smaller, but the patching curve 
still lies in this region; hence, the straight attached 
shock becomes smaller and smaller. Finally, a wedge 
angle @,, 1s reached such that \/, = |. For this value, 
Region A becomes a point, and the patching curve 
also degenerates to a point. For @ > 6,,, if the shock 
is still straight and still attached, then because VW. < 1, 
the patching curve must lie in the Region B. 

Now, we see from Fig. 14b that the velocity of the 
curve AB is greater than g, the velocity of the fluid 
behind the straight shock. Therefore, AB must 
sweep over fluid particles that passed through the 
curved portion of the shock AC. This means that 








Curveo 








Fic. 15. 


AB sweeps over a nonisentropic fluid, and, if AB js, 
patching curve, this fluid must then be of unifory 
velocity and isentropic. But the energy equation 


DS/Dt = 0 


states that in the disturbed field the entropy does jo 
change. We are thus led to a contradiction of require 
ments and must conclude that the straight attached 
shock solution for 4/2. < | shown in Fig. 14b cannoj 
exist. 

However, an attached curved shock solution as 
shown in Fig. 15 can exist for those cases in which 
Ms <1; the shock slope at the wedge point is the slope 
of the steady-state straight shock for the wedge angle 
6 and Mach Number .V/,, but there is no region of unj 
form flow behind the shock, and the flow is rotational 
throughout the entire field. 

This effect may be thought of as representing a 
discontinuity in the solution as follows: For J, > |, 
we have the attached straight shock, with a uniform 
region corresponding to the steady-state solution 
However, for \/, < 1, this quasi-steady state attached 
straight shock solution cannot exist; hence, at M/ = |, 
we have a discontinuity in this type of solution. 


CONCLUSION 


The following has been done in this paper: (1) The 
mathematical nonlinearized equations (neglecting vis 
cosity and heat conduction) have been obtained in 
conical form; (2) an invariance relation for these 
nonlinear equations has been proved; (3) it was shown 
that an attached shock, two-shock pattern is not 
possible for air; (4) the possibility of the existence of 
reciprocal flows has been discussed; (5) the patching 
curve relations were obtained, and with the aid of these, 
(6) the existence of a form of transonic singularity in 
the straight attached wedge shock was _ obtained; 
(7) 1t was shown that a certain type of attached shock 
flow cannot exist. 
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A Study of an Analogous Model Giving the 
Nonlinear Characteristics in the Buckling 
Theory of Sandwich Cylinders’ 


CHI-TEH WANG? anv G. V. R. RAO$T 
New York Unwersity 


SUMMARY 


In previous papers,'? the problem of buckling of sandwich 
cylinders under axial compression has been investigated both 
theoretically and experimentally. For sandwich cylinders with 
weak core, it was found that the buckling load given by the small 
deflection theory checks well with the experimental one. Whereas 
for homogeneous cylinders, it 1s well known that this is not the 
case, and the buckling phenomenon must be explained by the 
nonlinear large deflection theory. The limiting case of a sand 
wich cylinder with a strong core is an homogeneous cylinder. 
Thus the weakness of the shear rigidity of the sandwich cylinder 
must have compensated the importance of the nonlinear terms in 
the large deflection theory. To show this effect, nonlinear theory 
must be applied to sandwich cylinders with various degrees of 
weakness in the shear rigidity of the core. This, however, is a 
difficult problem. With a view to getting a qualitative picture in 
this respect, a model originally proposed by von Karman, Dunn, 
and Tsien* which will reproduce the nonliner behavior in the 
buckling phenomena is studied, and a satisfactory explanation is 


obtained. 


INTRODUCTION 


I PREVIOUS PAPERS, ! ? it was found that for sandwich 
cylinders with weak core the theoretical buckling 
loads based on the small-deflection theory agreed well 
with the experimental values. The limiting case of a 
sandwich cylinder with strong core is a homogeneous 
thin cylinder, but in that case it is well known that the 
classical linear theory fails to predict the buckling stress. 
Many investigators have attempted to explain the dis- 
crepancies between theoretical and experimental results 
in this latter case. Although the question of just how 
to define the ‘‘buckling load’”’ is still a subject of con- 
troversy,*~® it is clear, however, that the reason for this 
discrepancy is due to the fact that the nonlinear terms 
in the governing differential equations which have been 
left out in the classical theory become important here. 
By introducing second-order terms in the shell theory, 
von Karman and Tsien found that there exists a buck- 
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ling phenomenon entirely different from that given by 
the linear theory. Based on these results, Tsien* ® ‘ 
proposed new definitions of the buckling loads which de- 
pend on whether the load was applied by a rigid testing 
machine or dead-weight loading. On the other hand, 
Donnell and Wan’ have suggested that the buckling 
load should be the maximum load that a cylinder with 
initial imperfections can sustain. 

The investigation of the buckling of sandwich cylin- 
ders introduces another factor into this already com- 
plicated phenomenon—namely, the influence of shear. 
A study of the effect of shear on the buckling load is not 
only important in the buckling theory of the sandwich 
cylinder itself but may also throw some light on the 
buckling behavior of curved structure in general. A 
rigorous treatment of the buckling of a sandwich cylin- 
der by the second-order approximate theory should be 


based on the large deflection equations. These equa- 
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tions have been formulated in reference 2, but in that 
case five nonlinear differential equations must be solved 
simultaneously. 
these equations will be exceedingly difficult and in- 
volved. 


However, the essential characteristics of buckling 
phenomenon retaining these nonlinear terms in the 
buckling equations can be studied by an analogous sim- 
ple model that does not involve any mathematical dif- 
ficulty in its treatment. In the paper by von Karman, 
Dunn, and Tsien,* it was shown that the typical fea- 
tures of the buckling of curved shells can be reproduced 
by the buckling of a column with nonlinear lateral sup- 
ports. This follows from the fact that a longitudinal 
strip of a cylindrical shell may be considered as a column 
supported on rings (Fig. 1). With finite deflections, 
these rings behave as nonlinear springs. Following the 
same considerations, it appears that the nonlinear be- 
havior of the buckling of sandwich cylinders may also 
be seen from a study of sandwich columns supported on 
sandwich rings. Based on this analogy, a study is 
carried out in this paper to show the effect of shear 
rigidity of the core on the nonlinear characteristics of 
the buckling phenomenon. 


THE CHARACTERISTICS OF THE NONLINEAR SUPPORT 


In the original paper by von Karman, Dunn, and 
Tsien,* the nonlinear support was considered to be 
derived from a semicircular ring clamped at both ends. 
A theoretical treatment of this problem has been carried 
out by Biezeno and Koch.’ In the actual buckling of 
thin circular cylinders, tests have shown that there exist 
somewhere from 16 to 20 buckles in the circumferential 
In the case of sandwich cylinders, a similar 
Therefore, 


directions. 
number of buckles has also been observed. 
it appears to be a better approximation to consider the 
nonlinear support to be derived from a segment of 
circular are with simply supported ends rather than the 
semicircular ring with clamped ends. At the same time, 
the treatment of a circular are with simply supported 
ends is a much simpler problem. 


It is evident that the actual solution of 
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The load-deflection relation of an arch with smalj 
curvature and with hinged ends has been studied in de. 
tail by Timoshenko.® A small segment of circular are 
can be treated with good approximation as a flat 
arch, the initial centerline of which is given by the equa- 


tion 
vy. = a sin (rx/I) (1 


where a is the rise of the arch and / is the projected 
length (Fig. 2a). Under loading, deflections y» will be 
produced so that the final ordinates of the deflection 
curve are y = y,; — ye. Assume that y2 can be repre- 


sented by the sine series 


to 


yo = >> a, sin (nrx/l) (: 
n=1 

where the undetermined parameters a, can be deter- 
mined by the variational method, which is more usually 
referred to as the energy method. 

Denote the shear at any section of the arch by Q. 
The change of curvature produced by the bending 
moment and the shear force is 


d*y, dy d*y._ 1 ( 
dx? dx? dx? EIN’ 


ak] ) . 
(3 
AG dx 


where F is the Young’s modulus, / is the moment of 
inertia, G is the modulus in shear, A is cross-sectional 
area, and a is a factor with which the average shearing 
stress must be multiplied in order to obtain the shearing 


stress at the centroid of the cross sections. For homo- 


geneous beams with rectangular cross sections, a = 1.5. 
For sandwich sections with thin faces and weak cores, 
1. In the following 
calculations, for simplicity a will be taken as 1. 
Eq. (3), 


a iS approximately equal to 
From 


M = EI [(d*y2/dx?) — (1/AG)(dQ/dx) | (4 


In order to include the effect of shearing force on the de- 
flection of arch, the concentrated load F should be con- 
sidered as the limiting case of a uniformly distributed 
load g over an extremely short portion ¢ of the beam 


(Fig. 2b). The shear at any section x is, therefore, 
Q = (F/2) — H(dy/dx) 
when 0 < x & (1/2)(/ — €) 
Q= 4 — H = ~ a(x = Z - “) 
re dx 2 2 (5) 
when (1/2)(1 — €) < x (1/2)(1 + e) 
Q = —[(F/2) + H(dy/dx)], 


when (1/2)(1 + €—§ <x <1. 


The bending strain energy is 


1 f' i EI ['(d*y.2—-1«dQ\? : 
dx = — - — —— dx (6) 
2/0 EI 2 0 dx? AG dx 
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The shearing strain energy is and the work done by the lateral load F is 


iz Cag Wr = F(y2)e=1/2 (9) 


- AX (7) 
2/0 AG 


Vv, = 
The ‘principle of minimum potential energy requires 
that the potential energy LU’ 
The work done by the axial forces /7 due to end shorten- , ; ; : : 
a U=V,+ V, — Wu — Wr (10) 
ing 1S 

to be a minimum. The coefficients a, may be deter- 


>) ‘ ° ¥ ; es pie rere ae 
Wy = : H | (%2) a () | dx (8) mined by setting the first variation of l’ with respect to 


dx dx ad, to zero. Carrying out the calculation, we obtain 


a Eee ‘a =) (2m 4 Fl )] fi HP ( P ey] 
me Lee PAG)” \wEI AG El PAG | 


(= an) _ ar | ‘| ‘( “) | 
a, =— — ge ae | n*in*tl — -}— : 
wiKI tT°*AG 2 j AG rkI 


In obtaining the above relationship, we have used the condition that as « approaches zero, ge > F. 


(11) 


Taking only 


one term in the series, the rise of the arch under load becomes 


7 ( aL )] | ~“, rely) (12a) 
oo = US ~ \48.7e1 * 4.934G/ } / rE! PAG ” 


In the subsequent calculation, we shall only retain one term in the series; in such a case, a better approximation 
may be obtained by replacing (a — a,) in Eq. (12a) by the exact value of maximum rise, which is 


( FR ¥ FI )] [; Hr (: ‘ =) | (12b) 
(a — ae io | SS - 
1 1) exact a 18K] 4AG rk l?AG 


The centerline of the arch after loading then becomes 


7 fi FIs (: ‘ sd) | _* ‘) if aE (, m a ‘ci 
.~ 48aEI PAG) Io) | El PAG 


The deflection curve, Eq. (13), of the arch, however, is still undetermined unless HW is known. The equation for 


calculating the axial load H is obtained by putting the change in length of the span due to deflection equal to the 


compression of the arch due to axial load. Assuming for flat curves that the compressive force along the length 


of arch is constant and equal to //, we obtain 


7 


Hi 1 [' : #*. 
‘ER sf (dy,/dx)* dx — sf (dy/dx)* dx (14) 


Substituting Eqs. (1) and (13) into Eq. (14) and integrating, we obtain 


Hl _ way [ FIs (: - =} If i, PM oa ae 
AE 41 \ 48aEI rAG/ | / rE PAG) | § 


A sandwich-type structure consists of two external layers of thin high strength material and a thick internal 
layer of lightweight weak material. The former are usually called the ‘‘face layer’’ and the latter “the core.”’ 
For such type of structure, AG in the above formulas may be replaced by A.G., where the subscript c refers to the 
We shall not neglect the normal stress carrying ability of the core as usually assumed, in 


properties of the core. 
From the analysis of composite 


order that our final equations will reduce to those of homogeneous structure. 


beams, it can be easily shown that the flexural formula 1/p = M/EI for homogeneous beam becomes 1/p = 
The 


M/(E,l, + E.I.), where p is the radius of curvature and the subscript f refers to properties of the face layers. 
axial compressive strain H//AE should be written as H//(A,E,; + A.E,). E, may be expressed in terms of G, by 
using the relation E, = 2(1 + v,)G,, where v, is the Poisson’s ratio for the core material. 
Writing EJ, + EJ, = EI, we obtain, for sandwich beams, 


am _7 2"), [ FI* (: ‘. 12k] )f [fh {1 . )\$ - 
AE, + 21+ JAG. 4 | iSaEI PAG.) | | EI PAG.) | § Sa) 
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Eq. (15a) reduces to the equation for the homogeneous beam with finite shear rigidity if we let E. = E,, v, = y, 
[,+I,=I1,and A = Ay+A,=A,. Let ¢ be the thickness of the face layer, 2h be the thickness of the core. 
and b be the width of the beam. Let t < hk. Thus, 

A, = 20t, A. = Doh 

IT, = 2bth?, I. = 2bh?/3 


To simplify the writing, let us introduce the following symbols: 


u Fl*/48aEI, 
Ki = £,/G,, 
Kz <= a’ [2, 


II 


B = Hl?/z*EI 
Ke = AVA; = t/h 
K, = h*/a?* 


m = 4EI/E,Aa*® = 4(E,l, + Eel,)/EjA.a? = 4K4[Ke + 2(1 + v,)/3Ki] 


In terms of these parameters, Eq. (15a) becomes 


mB 
Ky + 2(1 + »,)/Ky 


For given values of Ay, Ay, Ay, Ay, and »,, 8 can be solved in terms of u from Eq. (16). 


= 1 — [1 — u(1 + 3mK,4A3)|?/[1 — BU + 0.25r2mK1K3) |? 6 


With 8 known, the de- 


flection 6 of the arch measured from the unloaded position may be found from Eq. (13) to be 


6/a = 1 — [1 — u(1 + 3mK,K3)]/[1 — BO + 0.259?mK1 Ks) | 17 


For numerical calculation, the circumference of the cylinder is assumed to buckle into 20 buckles, which gives 


a/l = ~/ Ks; = 0.04. 


For K, = 2.6, we take », = » = 0.3. 


For A; = 500, 1,000, and 5,000, we take v, = 0. 


Taking As = 0.1 and Ay = 0.2, the nondimensional load deflection curves of the arch are computed from Eqs. (16 


and (17) for these values of A. 


BUCKLING OF A SANDWICH COLUMN WITH NONLINEAR 
LATERAL SUPPORTS 


Having determined the characteristics of the non- 
linear support, we shall now proceed to treat the prob- 
lem of buckling of sandwich columns with such lateral 
Let the column be initially straight and its 
Assume that 


supports. 
centerline be taken as the x-axis (Fig. 4). 
the column is supported laterally by g nonlinear springs 

Let 6;, 6s, . . ., 6, be the deflections at 
Depending upon the load-deflection 


GE Hig Bay «. 6 sy Bas 


these points. 
characteristics of the springs, the spring forces at these 
points are F\(6,), F2(d2),..., F,(6,). The elastic energy 
Vp stored in these springs at the deflections 6), 62, . . ., 6, 
can be easily calculated as 

*dm 


F,, (é)dé (18) 


Vo = 2, 


m 1 0 


Assuming the column as simply supported at the two 
ends, the deflection of the column may be represented 
by the following Fourier series: 


y= > asin (nrx/L) (19) 


n l 


where a, are again undetermined parameters that may 


be ascertained by the variational method. The de- 
flection of spring m at x = x,, is, therefore, 
im = V(Xm) = > a, sin (mrx,,/L) (20) 


n=1 


e 


The curves showing w vs. 6/a are plotted in Fig. 3. 


Similar to Eq. (+), the bending moment at any section 
of the column may be written in terms of the curvature 


/2y ) 
M = El, ( a. -) (21 
dx? A.G.dx 


where /, is the moment of inertia of the cross section of 


as 


the column. 
we shall consider the concentrated load F,, as the limit- 


Again in the computation of the deflection 
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x ing case of a uniformly distributed load g,, over a short 
portion ¢€, of the column, and 


lim gmém = Fr. 


em—>0 


The shear at any section x is, therefore, 




















p 
nh 
t =. dy 
Q= r,, ~ P— 
} 2 az dx 
when 0 < x & [x — (4 /2)] 
1 ] > (22) 
: ay € 
Q = > F,, = P = = @ (« — Fi + ) 
2m=1 dx 2 
when [x, — (6/2)] < x & [mi + (6/2) 
8, l <= dy 
J=-) F, —-P—-F 
Ff o— OW 3 .* 
2 
| L when |x; + (€ 2)] SxS [Xe — (e/2)],..., ete. 
| The bending strain energy is 
“ El, [{" (d*y 1 dO\? 
| Xp 2l/s Vy = Z (; 7.2) @ 
z 0 dx? A.G,dx 
| 8, and the shearing strain energy is 
F —__ | | *] 
Vy. = / )? dx (24) 
2A.G, Jo ( 
| The work done by the axial compression P applied at 
Y | X=L/3 the ends of the column due to the lateral deflection is 
t 
: ag “ (dy\? an 
W, = dx (25) 
2 JO dx 
Y= The potential energy L’ in this case is, then, 
U = Vg + Vs; + Vsp — Wa (26) 
P and the coefficients a, may again be calculated by 
setting the first variation of lL’ with respect to a, to 
Fic. 4 ZeTO 


Now, consider the case of sandwich column supported on two equally spaced springs with the same characteris- 
tics. Carrying out the calculation and noting that Fj = F, = F, we see that the condition 0U’/0a, = 0 leads 
to 


EI,L P \2/nr\! PL p nr \? : r _ nn _ 2Qnr P 
(1 — ) ( )a, ( - )( ) a,+ Fli-— : (sin + sin + El,{ 1 — —— i 
2 AG, 8 2 A.G, L AG, 3 3 AG, 


F (nr\?/. nr _ 2nr FL (OF, . nx . OF, . 2nr P /OF,. . me OF. . 2nr 
) (sin + sin ) + ( sin + sin ~ 4 6; sin + Ds dy sin = 0 
3 3 5 5 y. bo ‘ 


A.G. ¥, 0 3A.G, 06; 5 O06» ” .G, 06; 3 
where the relation 06,,/0a, = sin (nx,,/L) has been used. The coefficient a, can therefore be written as 
ae a NTX yp, ) j RO(F,,/P-) Pp | b./L OF. /P,) 1 
sin (1 + ?R) + — R 1+nR+— 
An vA > P, * ( 4 { t 3 O( dm L) P. F,, P, O(b», L) f 


Lr ( RP\? , P/P. ] 
1 — ) n- — n* 
P, 1 — (RP/P.) 











where, for simplicity of writing, the notations P, = r°EI,/L? and R = P,/A.G, have been introduced. Eq. (27 


gives the coefficients a, if the Fn, O0(Fin/P.)/O(6m/L), and (6m/L)/(Fn/P,) are known. The quantities, however 


again depend ona,. To determine a, uniquely, we may proceed as follows: Recalling Eq. (20), 


6,/L = > (a,/L) sin (nrx,/L) (28) 


n=1 


Substitution of Eq. (27) into Eq. (28) gives 


6; : Fy R re) F,, é P, P bm L re) Fe P, ) 
Me Be oes 4 1 Rei eS | 
# m=1 Ho \ 3 O(b ny L) P, Fy P, OS» L) 


( , , Sin (n1Xm/L) sin (nx; £)) ‘ F 7 
——.. > of PrP. R— ; 
[x°[1 — (RP P,) |? =1 | | \ ~, P. [1 — (RP P.)\~* 


> sin (nrx,/L) sin (nrx;/L) 


= P/P, (29 
F — (RP/P,) "| 
If we define 
K. = 2 - - sin (mmx,/L) sin (nrx,/L) 
™ {1 — (RP/P,)]? n=1 n2[(P/P,.)/1— (RP/P,)] — n4 
(30 
eT oe 2 sin (mmrx,/L) sin (nrx,/L) | 
™ {1 — (RP/P,)] n= ((P/P.)/1 — (RP/P,)] — n? 
Eq. (29) can then be written as 
bs : . Faf RO(Fn/P.) =I bn/L OCF m <I} ai Fn 
= Kim —41 + = ~ 1+ + DL Kim’ R s=1,2) (31 
L » P, \ 3 0(6,/L) P, Fn/P. O(8m/L) 4§ * &, P. si | 


With a given value of P/P., the coefficients Kym, Ksm’ can. be calculated. With A,,, and K,,,’ calculated and 
the spring characteristics given, Eq. (31) then consists of two simultaneous equations for 6; and 62. Solving, 6,, bs, 
then F, and F, can be calculated. The Fourier coefficients a, /L can now be determined by means of Eq. (27). The 
problem is thus completely solved. 


The coefficients Kym and Ks,’ may be seen to be symmetrical from their definition—i.e., Kym = Kin, Ksms’ = 
Kyns’. Also, they are independent of the spring characteristics. In the case of the two equally spaced springs, 
Ku, Ko, Ky, Ki’, Ke’, and Ky’ can be easily summed up by well-known relations® ' to be 


’ : 2 l/, 0 r r 
Ku = Ke = 9 — 40 scoté — — / st es Ro 
ee 0\] /P P 
Ky = Ka = | ~ (3 ese 8 — csc )| ] (: — R :) 
9 40 3 P, P, 
Siti, > alll a\ / P 
Kn’ = Ke’ = 3 cot 8 — cot (1 — R— 
40 3 re 
ate 0\ / P 
Ky’ = Ky’ = — —(3 csc 0 — ese (1 — R 
46 3 P, 


i] / ive 
where 0 = x i{(P/P.)/(1 — (RP/P.)]} “. 
Now assume that the deflection of the column is symmetrical. Since the two springs have the same character- 


(32) 


sitics, we have 6; = & = dand F, = Ff; = F. Eq. (31) reduces to 
6 ; i 4 RO(F/P,) P 5/L O(F/P.) ]\ F ai 
= (f Ky) =41 — R l - Ky’ + Ky’)R (33) 
i “Su + Au) by! T 3 OG6/L) al F/P, a6 rol + An + Salk 


From Eq. (32), we have 


: : | l/, 6 60 /P P 
KAyn+ Kp» = E = ul? tan 3 — tan ‘)] / a — R ) | 


Ky’ + Ky’ = : 3 t : t =) /(; R-) 
Au A220 = AG an 2 al” / — P. 
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(34) 
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WevVer, 


(28) 


(29 


(31) 
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BUCKLING THEORY OF 


The load-deflection of nonlinear spring as worked out 
in the previous section is too complicated for the solu- 
tion of 6/L when it is substituted into Eq. (33). It was 
found, however, that this relation can be represented 
with good approximation by the following power series: 


F i(*) m (:) ‘ y 
= 0— - 2 > 95 
P, L a "ae (eo) 


The constants 6, c, and e may be determined by the 
initial slope and the positions of the maximum and 


minimum of the F vs. 6 curve. 


as defined before, 


Now, 
a FB - FL? nr] (+) a F "(=) K;? 
“ 48aEI wEI48a\L) ~ P.I \48/VK; 


where A; = //L. In terms of u and 6/a, Eq. (35) be- 


comes 


A)G) Cr) «= Na) 


or 


48/1, nn... 5\° =" : 
( Ju = »( )ke +e( ) V Ei Kot e( ) K3K;° 

r°\l a a a 
(36) 


Combining Eq. (35) with Eq. (33), we obtain 


() = n(2) + 9() +9(2) 4 
nen) on 


B, = }({1—(2RP/P.)|(KutKw) + R(Ku’ + 
Ky’)} b’ + [R(Knu + Ky)b’?/3] 


By _ } [1 am (83RP P.) (Aun + Ky ) + 


where 


R(Kiu’ + Ky’ bie’ + R( Ky + Ky)b’c’ 
B; = }f1 — (4RP/P,)|\(Ku + Ky) + 
R(Ky)’ + Ky’) te’ + 2PR(Ky + Ky) x 
(2b’e’ + c’*)/3] 
By = 5(Ku + Ky)R c’e’/3 
By = (Ku + Ky)Re” 
where b’ = b//I,, c’ = cI/Ih, ande’ = el/I,. A trivial 
solution of the equation is 
6/L, = @ (38) 


which represents the unbuckled straight column. Other 
solutions may be obtained by dividing Eq. (37) through 
by the common factor 6/Z and then solving the resulting 


fourth-order equation in 6/L. For any given value of 
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P/P., (Ky + Ky) and (Ay,’ + Ay’) can be calculated 


from Eq. (34) and (6/L) can then be solved from Eq. 
(37). 

In terms of the constants defined previously, we find 
that 


R = P./AG, = 7EI/L°A.G 


= w(Eql, + Ecle)/L? 


5 "(ey l Nee 2(1 + "2 
= x? : -— 4 
(| l (; G.A, 3 

= 9°K;KiKs?{ KiKe + [2(1 + »)/3]} 
we shall take K. = 0.1, K3; = 0.0016, Ky = 0.2, Ks = 
0.1, and J/J, = 2.77 X 10-4. For K, = 2.6, 500, 1,000, 
5,000, we have R = 3.56 XK 10~°, 1.60 & 107%, 3.16 X 
10-2, and 1.58 X 10~*, respectively. From Fig. 3, the 
constants b, c, e are then computed. To compare the 
characteristics of the P vs. 6 curves, we plot P/P, 
against 6/L in Fig. 5, where P, = r°E//L’, EI being the 
EI of the homogeneous column. 


A Ge 


CONCLUSIONS 


From the P/P, vs. 6/L curves, it is seen that, as the 
shear rigidity of the core of the sandwich structure is 
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reduced, influence of nonlinearity of the equations on 
the characteristics of the buckling decreases. Because 
of the decrease of shear rigidity of the column, the 
classical buckling load is greatly reduced, while the 
“valley” in the P vs. 6 curve due to retaining the non- 
linear terms does not lower as much. This explains 
that, for sandwich cylinders with weak cores, the 
buckling loads as calculated by the small deflection 
theory check well with those found by experiments; 
whereas for sandwich cylinders with strong core, this 
agreement cannot be expected. 
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Theoretical Investigation of Several Types of 
Single Degree of Freedom Flutter 


HARRY L. RUNYAN,* HERBERT J. CONNINGHAM,* ann CHARLES E. WATKINS* 
Langley Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


The possibility of an undamped oscillation of an airfoil (or 
aerodynamic body) with a single degree of freedom in potential 
flow was first noted by Glauert in 1929 but has received relatively 
little attention until recently. These undamped oscillations 
might have a direct bearing on problems of both ‘‘dynamic sta 
bility’ and ‘‘flutter’”’ of airplanes. Because of the many con- 
figurations of airplane design now being considered and because 
of the wide operational ranges of speed and altitude, knowledge 
of these single degree phenomena becomes more significant. 

The present paper gives theoretically predicted effects of 
several primary parameters associated with single degree of free 
dom flutter. Specific problems dealt with are: (2) pitching os 
cillation of an unswept wing, oscillation of an aileron about its 
hinge axis, and bending oscillation of a swept wing in subsonic 
streams; (b) pitching of an unswept rectangular wing, pitching 
of a delta or triangular wing, and oscillations of an aileron about 
its hinge line, in supersonic streams 

Results for one series of experiments which confirm the ex 
istence of single degree of freedom pitching flutter of an unswept 


wing at low speeds are included. 


INTRODUCTION 


[' HAS BEEN RECOGNIZED for some time that the 
linearized aerodynamic theory for an airfoil, os- 
cillating in a two-dimensional potential subsonic or 
supersonic flow, indicates that under certain conditions 
a single degree of freedom pitching instability is possi- 
ble. The possibility of occurrence of this instability, 
which may be referred to as a type of single degree of 
freedom flutter, was first noted by Glauert' about 1929. 
Glauert, following Wagner,” made theoretical studies of 
the air forces acting on an airfoil oscillating in two- 
dimensional incompressible flow. Although Glauert 
discussed briefly the conditions under which the theory 
for incompressible flow indicated that such an insta- 
bility might occur and went so far as to propose an ex- 
perimental investigation to test the theory, it appears 
the problem was completely overlooked for several 
years. 

In 1937, Possio* published the results of a theoretical 
study of air forces acting on an airfoil oscillating in 
steady two-dimensional supersonic flow. In his publi- 
cation it was pointed out that, for certain ranges of 
Mach Number and location of axis of rotation, a pitch- 
ing instability is indicated by linearized theory for two- 


dimensional supersonic flow. There followed consider- 


Presented at the Aeroelasticity Session, Nineteenth Annual 
Meeting, I.A.S., New York, January 29-February 1, 1951. 


* Aeronautical Research Scientist. 


able interest in the possibility of single degree of free- 
dom flutter in supersonic flow, and a number of in- 
vestigators'~? both in this country and in England 
studied the problem. 

In 1946, Garrick and Rubinow® presented the results 
of a general investigation of oscillating air forces and 
flutter in two-dimensional supersonic flow. Their in- 
vestigation was partly devoted to the problem of single 
degree of freedom flutter. Ranges of Mach Number 
and of location of axis cf rotation were delineated for 
which single degree flutter was predicted to be possible. 
Furthermore, these authors pointed out again that a 
similar type of instability is indicated by the theory for 
incompressible flow. After this, Smilg* made a calcula- 
tion study to determine the conditions, with regard to 
moment of inertia and location of airfoil pitch axis, 
under which the instability might occur in incompressi- 
ble flow. 

The results of Garrick and Rubinow® for supersonic 
flow and of Smilg* for incompressible flow may be con- 
sidered as starting points for some recent®~'! and current 
work on the problem of single degree of freedom flutter. 
In some of these investigations the possibility of un- 
damped oscillations other than pitching of the wing has 
been studied. 

The purpose of the present paper is to discuss some 
results and salient features of these investigations. 
Theoretically predicted effects are given of several pri- 
mary parameters associated with single degree of free- 
dom flutter. Specific problems dealt with are: (a) at 
subsonic speeds, pitching oscillation of an unswept 
wing, oscillation of an aileron about its hinge axis, and 
bending oscillations of a swept wing; (b) at supersonic 
speeds, pitching of an unswept rectangular wing, pitch- 
ing of a delta or triangular wing, and oscillations of an 
aileron about its hinge line. 

Results are included for one series of experiments 
which confirm the existence of single degree of freedom 
pitching flutter of an unswept wing at low speeds. 

Until recently, whatever interest was shown in the 
phenomena of undamped single degree oscillations was 
mainly academic. However, with present-day con- 
figurations of airplane design for high speeds at high al- 
titudes, undamped single degree oscillations even of ex- 
tremely small amplitude may become important. 
Furthermore, calculations of single-degree flutter may 
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be a useful easily obtained limit for cases of coupled 
flutter. 


SYMBOLS 


a = nondimensional distance of axis of rotation 
from midchord, measured in half chords, 
positive aft 


Ao, A1 aerodynamic flutter coefficients 

b = airfoil half chord 

Bo, By = mass, inertia, and elastic force coefficients 

c = nondimensional location of control surface 
hinge axis, measured in half chords from 
the midchord, positive aft 

e = cotA 

Cc. = structural stiffness per unit of span length 
(z = a, B, h) 

d = coefficient of viscous damping 

dy = 1/EZ’ 

d; = defined in section following Eq. (5-2) 

E’ = complete elliptic integral of second kind 

af = complete elliptic integral of first kind 

F,G = functions of & for oscillating plane flow 

g, = coefficient of structural damping (z = a, 

; B, h) 

h = vertical bending displacement 

is = moment of inertia about axis of rotation 
per unit of span length (z = a, 8) 

k = reduced frequency parameter (bw/v) 

7 = wing semispan 

M = Mach Number 

m = mass per unit length 

M, = complex aerodynamic moment about axis 
of rotation (s = a, B) 

M,, Mo, M3, My = aerodynamic flutter coefficients 

Fe = aerodynamic lift force 

v = velocity, ft. per sec. 

Xe = location of axis of rotation, based on 
chord, measured from leading edge, 
positive aft 

y’ = spanwise coordinate 

£ = a single degree of freedom (z = a, 8, h) 

be te a &e = aerodynamic flutter coefficients 

a = angular displacement of airfoil about axis 
of rotation a 

B, = VM -1 

B = angular displacement of control surface 
about axis of rotation c 

w = circular frequency 

o, = natural circular frequency (z = a, 8, h) 

A = sweep angle 

n nondimensional spanwise coordinate, based 
on semispan /’ 

p = air density 

Subscripts 

r = reference 

R = root 

h = bending degree of freedom 

z = a single degree of freedom 

a = pitching degree of freedom 

B = control surface degree of freedom 


Basic CONSIDERATIONS 


Before discussing specific examples of single degree of 
freedom flutter, it may be advantageous first to review 
the concept of a single degree of freedom vibrating sys- 
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tem and then to show the relation of this example to an 
aerodynamic system. The linear differential equation 
for a free system consisting of a mass m, a spring having 
a spring constant C., and a viscous damper having a co. 
efficient d is 


mk + dx + Cx = 0 (j 


The motion represented by this equation is damped if 
d is a positive quantity, as is ordinarily the case. If { 
should be negative, the motion is undamped and a con. 
dition of dynamic instability exists, and, if d is zero, 
harmonic oscillations may exist, corresponding to , 
borderline condition between damped and undamped 
motion. 

For a system such as an aircraft wing oscillating in g 
steady air stream, the same type of equation would 
apply as for the mass-spring-damper system previous] 
mentioned. However, the coefficients m, d, and C, of 
Eq. (i) will now have added components associated 
with the aerodynamics. A typical equation of an os- 
cillating wing is 


Ie + (1 + 2,)C.z = M,(2, 2, 8, ...) (ii 


where \/, represents the complex aerodynamic moment 
that is a function in part of amplitude z, velocity 3, 
acceleration 2, a reduced frequency k = wb/v, location of 
axis of rotation, Mach Number /, and sweep angle A. 
Eq. (ii) can be separated into two components, one 
associated with the damping of the system (sometimes 
referred to as the imaginary part) and the other asso 
ciated with the flutter frequency and velocity as follows 
(sometimes called the real part): damping (out-of- 
phase component) 


—C.g, + Im(M,.) = 0 (iti 
and frequency equation (in-phase component) 
wt, — C, = —Re(M,) (iv 
where 
M, = Re(M,) + tuIm(M,) 


Eq. (iti) is equivalent to the vanishing of the damping 
coefficient d of Eq. (i) and thus represents a borderline 
condition between damped and undamped oscillations. 
This is the condition sought in flutter work. The 
flutter frequency and air speed can then be determined 
from Eq. (iv). 


(1) SuBsoNIc FLow 


(1) Pitching Oscillations of an Airfoil 


The first example to be considered is that of an air- 
foil oscillating in pitch. As was pointed out in the in- 
troduction, Smilg* has reported some calculation studies 
for the Mach Number zero case in which it was shown 
that the axis of rotation must be ahead of the one- 


quarter chord. This work has recently been extended 
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by Runyan to include the effect of Mach Number and 
structural damping. 
Eq. (ii) in terms of the pitching (or yawing) degree of 


freedom a results in 


Tate + (1 + tga) Cea = M, (1-1) 


where \/, is the complex aerodynamic moment. M. 
for two-dimensional flow was given analytically by 
Theodorsen”? for the J = 0 case, and some numerical 


values were given by Garrick!’ for the compressible-flow 


case. 
If Eq. (1-1) is separated into two components, there 
results: (out-of-phase component) 


4 iL #3 ; i 
>. My; —_ ( + a)(a + me = Z:) = 
ke 2\2 2 


= M: + (“*) 2 ne 
Rk? 4¥i2 La me xpb' ei 


(in-phase component) 





Plot of reduced velocity coefficient v/bwg against inertia parameter J,/apb‘ for four axes of rotation (gq 
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(ce) « = -1.0. 
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f l a 
re M; — +a}{ Mi + Zs; — 5 21) 1+ 
l 
2 VW, = (0) 1-3) 
k* 
Eq. (1-2) represents the condition that the total 
damping, structural plus aerodynamic, vanishes. Be- 


cause of the transcendental nature of this equation, the 
solution must be obtained by step-by-step or graphical 
methods. The calculation procedure used here is to 
assume values of Mach Number .1/, structural damping 
14, location of axis of rotation a, and the inertia parame- 
ter /,,/rpb’ and to solve Eq. (1-2) graphically for k. The 
value of reduced frequency & thus obtained, taken with 
the corresponding assumed parameters, then constitutes 
a set of parameters under which a borderline condition 
between damped and undamped oscillation exists. This 


set of conditions can thus be applied to Eq. (1-3) in 











104 JOURNAL OF THE 
order to determine the frequency of oscillation w and 
the flutter speed v associated with these conditions. 
Some results of calculations based on the above pro- 
cedure for gg = O are presented in Fig. 1-1 where a 
flutter speed coefficient v/bw, is plotted against an 
inertia parameter /,/rpb‘! for three Mach Numbers, 
M = O, 0.5, and 0.7 and for four locations of axis of 
The effect of alti- 
tude may be observed since an increase in the inertia 


rotation a = —1, —2, —3, and —4. 


parameter /,/mpb‘ may represent an increase in altitude. 
The region of instability for each curve is to the right or 
above the curve, and the vertical asymptotes to the re- 
spective curves indicate the smallest value of the inertia 
parameter at which the oscillation could occur. The 
results of Smilg* correspond to the curves 1/ = 0. One 
of the more important trends to be noted by an examina- 
tion of curves is the large effect that Mach Number has 
on the range of the inertia parameter in which the insta- 
bility can occur. For instance, for the axis of rotation 
located ata = —4 (Fig. 1-ld), the minimum (asymp- 
totes) value of the inertia parameter for 1/J = 0 is 242,- 
000 whereas for J = 0.7 it is reduced to 3,223. If the 
aircraft center of gravity is regarded as the axis of rota- 
tion and the vertical tail is considered as the airfoil, 
values of the inertia parameter for modern aircraft 
vary from 2,000 to 20,000 at sea level, which are in- 
creased by a factor of 10 for aircraft operating at 60,000 
ft. It can be seen, then, that the parameters of modern 
aircraft may fall within the ranges of the parameter in- 
dicated here. However, it should be remembered that 
the calculations presented here are based on two-dimen- 
sional flow considerations and that such effects as as- 
pect ratio and other degrees of freedom have not been 
considered. It is probable that these additional factors 
could have a great effect on the single degree of freedom 
instability, and the results presented in this paper must 
be considered as possible limit conditions for actual con- 
figurations. 

The results of calculations in which structural damp- 
ing coefficient g. has been considered are shown on Fig. 
1-2 for one axis of rotation location a = —1.24. The 
abscissa is the inertia parameter, and the ordinate is the 
flutter speed coefficient v/bw,. The type of structural 
damping used in flutter analysis is in phase with the 
velocity but is proportional to amplitude. Damping has 
a great effect on the flutter speed, as can be seen by an 
examination of Fig. 1-2. As an example, the flutter 
speed for zero structural damping is raised by a factor 
of 5, with /,/mpb4 = 18,000 for ga = 0.02. 
the use of structural damping should be an effective 
method of eliminating this type of flutter. 

A preliminary experimental investigation at low air 


If possible, 


speeds of single degree of freedom oscillation has been 
made by Runyan in which a rigid wing was hinged in 
front of the leading edge (a = —1.24) and was re- 
The results of the investiga- 
The experimental data agree 


strained by a coil spring. 
tion are shown on Fig. 1-2. 
with theory for = O and g, = 0.008 at high values of 
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the inertia parameter which correspond to high alti. 
tudes. However, at low values of inertia parameter the 
agreement between theory and experiment is poor and 


requires further investigation. 


(2) Control-Surface Oscillations 

By a similar analysis, a single degree of freedom os. 
cillation can theoretically be shown to exist if the con- 
trol surface degree of freedom is examined. The equi- 
librium of moments on a section of a control surface 
about its axis of rotation may be written as 


18 + (1 oa igg)CgB — Ms ‘ey 


Theodorsen'? has given the expression for .\/g for the 
JJ = O condition. 
been given by Dietze!'! for one control-surface chord 


Tables for Wg for subsonic flow have 


wing-chord ratio, and these tables were extended by 
Turner and Rabinowitz'® to include several control 
Tables of Mg for 
supersonic flow have been published by Huckel and 


surface chord/wing-chord ratios. 


Durling.'® 


The results of calculations based on the values of \, 
from reports mentioned above are presented in Fig. 2-1, 
where the abscissa is the inertia parameter /,/mpb‘ and 
the ordinate is the flutter speed coefficient v/bwg. One 
ratio of control-surface chord to wing chord of 0.15 
(c = 0.7) is given for three Mach Numbers. These 
curves exhibit the same characteristics as previously 
noted for the pitching oscillations. A limiting (asymp- 
totic) value of the inertia parameter for each curve is 
indicated below which the configuration would be stable. 
Therefore, a control surface that is stable at low alti- 
tudes could conceivably become unstable at high alti 
tudes. Moreover, the effect of Mach Number is ex- 
tremely great and increases the region of the inertia 


parameter in which the oscillation could occur. For in- 
stance, at \J = O, the limiting value of the inertia 
parameter is 15 and is reduced to 5.57 at WW = 0.7 and 


is further reduced to 0.0079 for \J = 10/9. 

Recently single degree of freedom flutter of a control 
surface has been experienced in high-speed flight, and 
the explanation of the cause has not been entirely clear, 
although most investigators have been of the opinion 
that separated flow is involved. The mechanism of in- 
stability presented here, however, should provide a 
framework for a more complete and logical study, even 
though it is not necessarily a complete explanation of 
the cause. However, as for the pitching wing pre 
viously discussed, similar qualifications with regard to 
the effect of aspect ratios and other degrees oi freedom 
apply to the control-surface oscillations 


(3) Bending Oscillation of a Swept Wing 


This section deals with a cantilever swept wing per- 
forming only harmonic bending oscillations. The wing 
is shown in Fig. 3-1, and the equation of equilibrium ol 


forces for a wing section dy’, similar to Eq. (ii), is 
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Number, and reduced frequency k,, which is wb/v cos A. 
The semichord } is for a section normal or nearly normal 
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Fic. 1-2. Curves of experimental and theoretical single degree pitching oscillations for various values of structural damping (a = 
— 1.24) 
—~hm — C11 + ig,)h = —Pylh,h,...) (3-1) the swept wing, f,(y’) is the shape of the bending mode, 
; ; : , s 7 and 
where /i is vertical bending displacement of the section 
dy’, m is mass of a unit-length section, C, is effective 100 ; ae 
elastic restraint, and g, is structural damping coefficient | 
for bending vibration. P, is the aerodynamic lifting —— -. - oa 
- . ' 
force for harmonic bending oscillation and is a function so | th | 
of several parameters including angle of sweep, Mach | 
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to the leading edge. 
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The case of cantilever wing bending inherently treats 





a wing whose different spanwise sections have different 
amplitudes——that is, the wing has a mode shape that is 
not a constant over the span. An equation of equilib- 
rium of forces for harmonic motion, therefore, is not, in 














general, that of a representative section but is an equa- 
tion whose elements are integrations over the wing. 














The bending displacement of a section as a function of ad ” 





time and of span position is represented by : 


rer , 9« ” . — . 
= POL Tay 2) 3-2) Fic. 2-1. Plot of reduced velocity coefficient v/bwg against 
h Al f(y’) ] ( 3 ag 
‘ : _ inertia parameter /g/mpb' for one location of axis of rotation ¢ = 
where y’ is the coordinate measured along the length of 0.7 and for several Mach Numbers 
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Fic. 3-1. 





Swept wing. 


h = Ie (3-3) 
of which fp is a constant. 


In consideration of the work done by all forces in the 
bending direction, the equation of equilibrium of bend- 
ing forces for the entire wing is 


oy 
-i f m[ frly’) |?dy’ — 
0 
l’ bd AA 
h 7 Ci(1 + tgn)[fa(y’) ?dy’ = — | P»{ fr(y’) dy’ 
0 0 


The limits 0 and /’ are values of y’ for wing root and tip, 
respectively. 


Through a series of steps, Eq. (3-4) can be converted 
toa more useful form. The steps are: (1) Wing bending 
frequency w, is introduced to represent elastic restoring 
force; (2) substitute for P, the aerodynamic forces of 
Eq. (10) of reference 17 which apply for bending vibra- 
tion; (3) convert the lengthwise variable y’ to the non- 
dimensional yn, which is zero at the root and 1.0 at the 
tip; (4) replace h by —w°h; and (5) collect terms. The 
resulting equation of equilibrium of complex forces can 
then, as in preceding sections, be separated into two 
equations of component forces, expressing the simul- 
taneous equilibrium of forces in phase with bending dis- 
placement and of forces 90° out of phase wth displace- 
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ment. The two equations in symbolic form are: 


equilibrium of in-phase forces 


(By + Ao)rpb,*w?l’ = 0 (a: 


equilibrium of forces 90° out of phase 


(B, + A;)rpb,*wl’ = 0 


a 
9-0 


where By, is the integrated coefficient of mass inertia and 
elastic forces in phase with bending displacement, B, js 
the integrated coefficient of structural damping forces, 
A, is the integrated coefficient of aerodynamic forces jn 
phase with displacement, and A, is the integrated ¢o. 
efficient of aerodynamic damping force. The coefficients 
in expanded form as developed are 


Wh 2 sd b\? 1 . 

ho it« [Fa(n) }2dn 
w 0 b,/ k 
wn \? oa b\? 1 : 

B, = =—25 [ F,(n) ]*dn 
w 0 b, K 


"1.0 b 2 a 
Ag = = ( Ren{ Fn(n) }*dn — 
0 b, 
=) tan A | Lal (F dF, j 
onl (F,(n)] 
f F ii” oe 
s: (7) tan * l Fat a’F, _ 
7 (n)] j 
0 b ¥ b, k,,? wis dn? ”" 
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Fig. 3-2. Minimum mass parameter as a function of the sweep 
parameter to enable single degree bending flutter of a uniform 
swept wing 
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1.0 b 2 ' ; 
Ay iad _ A 5 Tn{ Fu(n) |? dy + 
"0 /b\3 tan A 1 dF, 
(-—1 Ren) [Fr(n)] 1 
J () '/b, Ry + Mal eg 


where x is rpb?/m mass parameter; 0, is the semichord 


of a section chosen as a reference, 
Fi(n) = frll’n) = fry’) 


R., and J,, are the real and imaginary parts, respec- 
tively, of a complex function, which accounts for part of 
the aerodynamic effects. The reduced frequency k, is 
wb/v cos A. 

Since a dynamic instability (w # 0) is being treated, 
the quantity within the parentheses in each equation, 
Eqs. (3-5) and (3-6), is zero as follows: 


By + Ao = 0 
B, + A, 0 (3-S) 


For incompressible flow, the aerodynamic functions 
R., and J., used are those of two-dimensional incom- 
That is, 


Il 


pressible flow. 


Re _ =e (2G h,): Las a (2F k,) (3-9) 


where the transcendental functions F and G of Theodor- 
sen account for the effects of the wing wake. 

An attempt to approximate some effects of compressi- 
bility at subsonic speeds can be made by substituting, 
for R» and J, values for two-dimensional subsonic 
compressible flow according to the normal-component 
Mach Number M cos A. For this purpose, the func- 
tions R,, and J,, of Eq. (3-9) are replaced by R,,’ and 
(1/R»)Icn’, respectively, in the notation of Garrick.'* 

The simultaneous equations can be solved for any two 
quantities or parameters when the other parameters are 
specified. For the present problem, it is convenient to 
solve for the sweep parameter tan A/(/’/b,) and an 
integrated mean or effective value of mass density ratio 


[- 4) l [F,(n) ]2d 
Ny ) « 
1 /0 (= K — ued 


= (3-10) 


hi a tot {> 
[Fn(n) |*dn 
0 b, 


The sweep parameter increases as the sweep angle A in- 
creases and as the length-to-semichord ratio /’/b, de- 
creases. The effective mass parameter 1/x, signifies 
also effect of altitude and increases with increased alti- 
tude for a given wing. 

Fig. 3-2 for the case of a uniform cantilever swept 
wing shows, as ordinate, the minimum value of mass 
parameter 1/x, as a function of the sweep parameter 
tan A/(l'/b,), for which single-degree bending flutter is 
theoretically possible in incompressible flow and in the 
compressible flow J cos A = 0.7. Conversely, for any 
mass parameter, the curves show the minimum value of 
the sweep parameter for which single degree flutter is 
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Fic. 3-3. Flutter-speed coefficient as a function of sweep angle 
for single degree bending flutter of a uniform swept wing of 
length-semichord ratio 6 and mass parameter 1/« of 1,000 


possible. For values above each curve single degree 
bending flutter of a uniform wing is possible; for values 
below the curves, it is not possible. 

It may be noted from Fig. 3-2 that, for a uniform 
wing with A = 45° and length-to-semichord ratio of 5 
(sweep parameter of 0.2), a wing mass parameter of at 
least 1,000 is necessary to enable single degree flutter in 
incompressible air. Such a large mass parameter is far 
above any contemplated value for actual aircraft. 
Present-day jet fighters have a wing section 1/«x of the 
order of 25 for sea-level air density and of 100 at high 
altitudes. The factor may be increased by heavier 
wings, higher altitudes, and weights on the wing. 

It is to be noted further from Fig. 3-2 that the nor- 
mal-component Mach Number M cos A = 0.7 causes a 
marked decrease of the minimum values of mass 
parameter and sweep parameter to enable single degree 
bending flutter. A wing whose sweep parameter is 0.2 
is predicted to need a mass parameter greater than only 
about 90, a practical combination. 

Fig. 3-3 shows curves of flutter-speed coefficient as 
functions of sweepback angle for a uniform wing with 
the parameters: 1,000; length-to-semichord 
ratio 6; and three different structural damping co- 
efficients 0, 0.01, and 0.03 in incompressible flow and in 
compressible flow with M/ cos A = 0.7. 

As was noted in Fig. 3-2 and in previous sections for 


l/xk = 


single degree wing pitching and aileron oscillation, com- 
pressibility is seen in Fig. 3-3 to have a marked effect 
on flutter characteristics. For example, for a sweep 
angie of 60° the flutter-speed coefficient for M cos A = 
0.7 flow is only about 31 per cent of that for incom- 
pressible flow. Consequently, for the two calculated 
flutter speeds to be the same, the wing stiffness in terms 
of bw, must in the cos A = 0.7 flow be about 327 per 
cent of that required for incompressible flow. 

Also in Fig. 3-3 structural damping is predicted to 
have a large effect on single degree flutter speed for the 
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particular wing parameters given. At a sweep angle of 
60° for an incompressible fluid, an increase of structural 
damping coefficient from 0 to 0.03 is predicted to in- 
crease flutter-speed coefficient by about S85 per cent. 
Such a large effect of structural damping is associated 
with the large mass parameter of 1,000. 


In addition to effects of mass parameter, sweep 
angle, and structural damping, there is an important 
effect of ratio of wing length-to-semichord which is not 
treated here. Also of interest is the fact that single 
degree of freedom bending flutter is the limiting case if 
bending-torsion flutter as the ratio of wing bending fre- 
quency to torsion frequency approaches zero under 
certain conditions on mass parameter and sweep pa- 
rameter. Those conditions are given for a uniform wing 
in Fig. 3-2. Some effects of permitting the frequency 
ratio w,/w, to increase from zero to finite values when 
the sweep parameter and mass parameter are such as to 
permit single degree bending flutter have been analyzed 
by Cunningham, but these effects are not within the 
scope of the present paper. 


In conclusion, it should be stated of this section that 
the three-dimensional flow effects of finiteness of span, 
of mode shape, and of compressibility obviously need 
more exact analysis than is presented here, particularly 
for wings of small length-to-chord ratio. Structural 
effects of the root of such a wing in their relation to 
bending vibration also need more rigorous treatment. 


(II) SUPERSONIC FLOW 


The final two sections deal with wings of rectangular 
and triangular or delta plan forms in supersonic flow. 
As mentioned in the introduction, the problem of un- 
damped pitching oscillations of wings in two-dimen- 
sional supersonic flow has been studied by a number of 
investigators.‘~’ In three-dimensional flow, the rec- 
tangular wing with Mach lines behind the wing di- 
agonals has been treated by Watkins® and by Stewart 
and Li.'"° The wide delta wing—that is, one with super- 
sonic leading edges—has been treated by Miles,'' Wat- 
kins, Nelson, and others. The narrow delta wing —that 
is, with subsonic leading edges—has recently been 
treated by Watkins. 


It is convenient to treat these cases by expanding the 
aerodynamic moment coefficient in terms of the re- 
duced frequency parameter k. Considerable informa- 
tion regarding undamped pitching oscillations can be 
obtained by retaining only terms of the moment co- 
efficient which depend on the first order of k. By such 
an approach one can indeed determine, with satisfac- 
tory accuracy, the range of values of such parameters as 
Mach Number, location of pitch axis, and aspect ratio 
in which the oscillations can occur. In order to study 
such factors as critical speed and frequency associated 
with the oscillations, more complete expressions for the 


moment coefficient are required. It is expedient here, 
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however, to discuss only the first-order effects of the 
reduced frequency. 
The equations to be considered are 

(1../mpb*) [1 — (wa/w)?] — Ms = 0 (II-] 
and 

([q/apb*) (w/w)? ga + My = 0 (II-2 
where .\/; and \/; are the in-phase and out of-phase 
components, respectively, of the moment coefficient, 
(4) Pitching Oscillation of Rectangular Wings 


For wings of rectangular plan form with Mach lines 
behind the wing diagonals, the components of moment 
coefficient, to the first order of the reduced frequency 


parameter, are 


I 2 — 3x0 
M, = —— | 38\(1 — 2%) — (4-1 
387k? A.R. 
2 
M, = 384k 28:3(1 — 3x9 + 3x0?) — Bi X 
2 B17X0(3X0 — 2) + lyn — 3 
(2 — 3x0) — (4-2 
2A.R. 


where x» is the distance, referred to wing chord, of the 
pitch axis from the leading edge; A.R. is aspect ratio; 
and B, is V M2? — 1. 

The range of values of V/, x), and A.R. in which un- 
damped oscillations can occur in this case is delineated 
easily by consideration of Eq. (II-2) with the damping 
coefficient ga zero. The condition that 1/; vanish, as 
may be seen from Eq. (4-2), leads to an algebraic equa- 
tion in 6), %», and A.R. which may be solved readily for 
any one of these quantities when particular values are 
assigned to the other two. 

Fig. 4-1 shows values of x» and / which satisfy this 
borderline condition for some selected values of aspect 
ratio. The curve for infinite aspect ratio corresponds 
to the two-dimensional flow case treated in other pub- 
lications.‘~** The curves for finite aspect ratio end in 
the lower right portion of the figure for combinations of 
aspect ratio and Mach Number for which the Mach 
lines from the wing tips coincide with the wing diago- 
nals. 

The regions inside the different curves of Fig. 4-1 
give the ranges of value of x» and .V/ in which single de- 
gree pitching instability might occur—that is, in which 
M, has a negative value. It may be noted that, as the 
aspect ratio is decreased, the ranges of values of both \/ 
and x», in which the instability could occur, are de- 
creased. For example, the curve for infinite aspect 
ratio indicates that the instability might occur at some 
Mach Number between | and a maximum of about 1.58 
and for any pitch axis located in an unlimited range 
ahead of the two-thirds chord position. The curve for 
aspect ratio 2, however, indicates that the Mach Num- 





of the 


(IT-| 


I[-2 


hase 
it. 


lines 
ment 
ency 


the 
tio: 


un- 
ited 
ing 
, as 
ua- 
for 
are 


his 
ect 
ids 
1b- 

in 
of 
ch 


r()- 
4 


}- | 


ch 








SINGLE DEGREE OF FREEDOM FLUTTER 109 


1.6 





Stable region 





1.5 











1.4 





















M13 
1.2 
Oe 
1.0 
-.5 ° -.3 ok ey o3 é 5 6 of 
XO 
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for delta wings with selected values of leading-edge sweep angle. 
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ber range is approximately only from 1.12 to 1.23 and 
that the location of pitch axis is limited from between 
0.275 chord (behind) and —0.25 chord (ahead) relative 
to the leading edge. Examination of Eq. (4-2) (or of 
more complete expressions for 7/4) will show that, as 1/ 
and xo are held constant, the value of J/, decreases 
(with x») not greater than about 0.67) as 1/2 A.R. as A.R. 
decreases. Hence, decreasing the aspect ratio has a 
highly stabilizing effect. It may be concluded there- 
fore that, with wings of finite aspect ratio, small values 
of the damping term ga(ga # 0) could exclude the possi- 
bility of the undamped pitching oscillations. 


(5) Pitching Oscillation of Delta Wings 

For delta wings, the components of the coefficient of 
total moment acting about the pitch axis x», to the first 
order of the reduced frequency k = bpw/v (bp is the root 
chord), are 


M3; = (mc’?/3k?)do(2 — 3x0) (5-1) 
My = 7 {3¢,(2M? — 1) — do(2M? + 1) + 
2B,*k 
Axod,[do —_ (2M? = 1) ] + txy?d, (AM? — 1)} (5-2) 


where c’ is the cotangent of the angle of sweepback A, 


I 1 — Bi2c’? 


do +) dy -_ 9 19m fe A 
E Bi2c"?F’ + (1 — 2Bc")E 


F’ and E’ are complete elliptic integrals of the first and 
second kind, respectively, with moduh Ft By2c"?. 
The location of pitch axis xo is referred to the wing root 
chord and is measured from the wing apex. 

Fig. 5-1 shows values of Mach Number and location 
of pitch axis for some particular values of angle of 
sweepback A for which the condition that M/, vanish is 
satisfied. Asin Fig. 4-1, the regions inside the different 
curves give the ranges of M and x» for which M, is 
negative (regions of instability). 

For all angles of sweepback, the range of location of 
pitch axis extends from the three-quarter chord to an 
infinite distance ahead of the wing apex (— ~ < x) < 
0.75). The pitch-axis location at the half chord (x) = 
0.5) shows an instability over the greatest range of 
Mach Number, from 1.0 to 1/2, and this range applies 
to a delta wing whose angle of sweepback is equal to, or 
less than, 45°. 

As the angle of sweepback increases, the range of 
Mach Number of the region of instability decreases. 
As the sweep angle increases to about 60°, the range of 
Mach Number in which instability might occur, as well 
as the magnitude of the negative damping moment, is 
so small that, here again, a small amount of structural 
damping could eliminate the possibility of instability. 
This behavior is in accord with that of rectangular 
plan-form wings in that reducing aspect ratio has a 
highly stabilizing effect on single degree of freedom pitch- 
ing oscillations. 
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It may be seen further in Fig. 5-1 that the upper-limi 
Mach Number for a zero damped oscillation (\/, = 9 
increases as the angle of sweepback decreases until the 
wing leading edge becomes sonic (Mach line coincides 
with the leading edge); but the Mach number for , = 
0 does not change with a further decrease of A. Thus, 
when the flow at the leading edges of a delta wing js 
sonic or supersonic (f;c’ > 1), the conditions under 
which M, vanish no longer depend on the angle of 
sweepback. Hence, the dashed curve represents the 
upper limit of Mach Number for which J, vanish for 
all wide delta wings. 


CONCLUDING REMARKS 


The main purpose of the paper has been to direct 
attention to the possible loss of damping, theoretically 
indicated, for various single degree of freedom systems 
rather than to give a complete discussion of the prob- 
lem. Each case considered is a limiting one of coupled 
flutter involving two or more degrees of freedom, and no 
attempt has been made to examine the effects of re- 
laxing the condition of infinite structural stiffness in 
other degrees of freedom. 

A study of the examples treated indicates that Mach 
Number has a marked influence and can cause a reduc- 
tion of the stable region as defined by a mass or inertia 
parameter. Structural damping for all cases considered 
raised the flutter speed and could be an effective means 
of reducing or eliminating the oscillation. For the ex- 
amples treated in supersonic flow, the range of variables 
in which the oscillation could occur was greatly reduced 
by a decrease in aspect ratio. The effect of aspect ratio 
has not been introduced in the examples calculated for 
subsonic flow, and it may be expected that large effects 
exist. This simplified treatment, however, should prove 
useful in connection with further experimental studies 
and in providing a logical grouping of variables. 
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Supersonic Flow Around Cones at Large Yaw 
I £ 
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The RAND Corporation 


SUMMARY 


Equations have been developed for the flow properties in the 
region bounded by a supersonic cone body at large angle of attack 
and its attached shock wave for use with the tabulations of 
Kopal.-* A comparison of the theoretical results is made with 
the experimental data of Cronvich.? A slight discrepancy 
that increases with increasing Mach Number has been found 
in the value of Apyy, the yawing head drag as tabulated by 


Kopal. 


SYMBOLS 


a = local velocity of sound 
( = speed of sound in a vacuum 
Cz = axial pressure force coefficient 
Cx = normal pressure force coefficient 
Cp = force pressure drag coefficient 
C1 = lift coefficient 
p = static pressure 
u,v, w = orthogonal velocity components 
V = resultant velocity 
B = angle of velocity vector with respect to the wind axis 
€ = yaw angle (equal to negative angle of attack) 
p = density 
barred values = nonyaw properties 
primes = derivatives with respect to 6 
Subscripts 
1 = free stream 
r = rotated conical field régime 
s = solid cone surface 
w = shock wave surface 
B = body axis coordinates 


INTRODUCTION 


De PRESENTATION IS INTENDED as an aid to aero- 
dynamicists and missile designers who will utilize 
Kopal’s tabulation of Stone’s solution for supersonic 
flow around cones at large yaw. Relations for the flow 
characteristics in the conical flow field and for the nor- 
mal and axial pressure force coefficients are presented 
in terms of the parameters used by Kopal. Since the 
general method of derivation is similar to that given in 
reference 6, many of the details are omitted. 

Under the direction of Kopal,'-* the Computing Sec- 
tion of the Center of Analysis at M.I.T. has tabulated 
results for supersonic flow around cones at zero angle of 
attack, at small angles of attack, and at large angles of 
attack. The presentation of the results for nonyaw 
cones is complete, and the analysis of the problem by 
Taylor and Maccoll‘ is generally available. The tabu- 
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lations for cones at small and large yaws contain all the 
essential results but are not in the form whereby the 
flow characteristics in the conical flow field are readily 
obtainable. Stone’s analysis® for slightly yawing cones, 
which is employed by Kopal, has been published and 
has been elaborated upon by Young.* The analysis for 
cones at large yaw is only briefly reviewed in the intro- 
duction in reference 3. A brief résumé of the present 
paper has been published in this journal.’ The purpose 
of the résumé was to draw early attention to the proper 
use of the M.I.T. tables. 


(GENERAL DISCUSSION 


Stone* obtained a second-order solution to the prob- 
lem of a supersonic cone at large yaw by assuming the 
following expressions for the flow properties, as the 
nonyaw variables (#, v, p, and p) are perturbed by the 
effects of angle of attack. See Fig. la for the wind coor- 


dinates as used in references 2, 3, and 5. 


ud) = ate dx, cosng te YS Um COS ME 


n=0 m=0 


>> Um COS MEY 
0 


v6) =vte Dy, cosng t+ é& 
n=0 m 
wd) = e z, sinng + & >> Wm sin me (1) 
n=0 m=0 


be Ei mcosng + & DY pm cosme 
n=0 m=0 


pO) = pte > &, COS Ng + € b Pm COS MY 


n=0 m=O 


Kopal® * evaluated the coefficients in the following cor- 


responding expressions 









rd 
i Conical surface 


(8 = const) 


Fig. la — Perturbed velocity components, Eq.(2), in the 
non-yawed conical flow field with wind 
axis coordinate system (r, 6, d) 
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v,(8,) 






u(O,) 
p(r,8,,) 





=e oreo Pica surface 


(0, =8+A8; 8=const) 


Fig. |b — Velocity components, Eq.(3), in yawed conical 
flow field with wind axis coordinate system (r, 6, ) 


u(d) = a + ex cos g + e(uy + U2 Cos 2¢) 
v0) = v+ ey cos g + e(% + v cos 2¢) 
w(@) = es sin g + e*w. sin a (2) 


p(d) = p + encos "(po + p2 cos 2¢) 
9 


p(@) = p + eEcos¢g ie €*(po + p2 cos 2¢) 


The coefficients x, up, M2, etc., are tabulated for various 
6 with cone angles and Mach Numbers as parameters. 
The determination of these coefficients essentially com- 
pletes the problem. However, upon examination of the 
manner in which the boundary conditions were formu- 
lated, it becomes evident that the properties 1(@), v(@), 
etc., apply only in the conical field of the nonyaw cone® 
In order to obtain the flow properties 
Taylor series 


(see reference 6). 
in the conical field of the yawed cone, a 
expansion is required, this being equivalent to rotating 


the nonyaw conical field to a yawed position. The ex- 
pansion is as follows: 
u,(0,) = u(6 + Ad) = 

u(O) + u’(0)A0 + u"(8)(Ad?/2) +... 
v,(0), = v(6 + Ad) = 

v(0) + v'(0)A0 + v"(A)(A6?/2) +... 

w,(0,) = w(6 + Ad) = ‘ded 

w(6) + w’(O)Ad + w"(6)(AG2/2) +... | ? 
p-(6,) = p(@ + Ad) = 

p(O) + p’(@)Ad + p”(6)(A6?/2) +... 
p(8,) = p(6+ Ab) = 


p(O) + p’(6)Ad + p”(6)(A?/2) +... 


where u,, v,, etc., are the flow properties for the yawed 
field and are functions of 0,, the yawed conical surface 
variable (see Fig. 1b). 

The expressions @, and A@ for the solid cone surface 
and the shock wave surface are well defined by Kopal.* 
They are 
6,, = 6, + Ab, = 6, + € cos g —(1/2)e cot 6, sin? %) 

(4a) 


and 


= 6, + Ad, = 8, + ea cos g + (By + Bo cos 2¢) 


(4b) 


(aL SCIBENCES—FEBRUARY, 1952 


Now the disturbing element is how to evaluate Aé fo, 
the intermediate conical surfaces. For small angles of 
attack, good results can be obtained by considering the 
intermediate surfaces to be rotated an amount € [see 
Eq. (5)]; these surfaces remain circular cones, as does 
the solid cone surface. 
tack, in the region of the bottom plane (g¢ = 180°), the 


However, at high angles of at. 


properties directly behind the shock wave assume a dis. 
continuity in slope. Also, the properties in a consider. 
able range of 6,, in the region of the top plane (¢ = 0) 

cannot be computed. To circumvent this situation an 
alternate method is proposed. The intermediate conical] 
surfaces are considered to be distorted linearly, from 
the circular cone body to the distorted shock wave. 
Both methods are illustrated in the next section, and 4 
comparison of the flow properties are shown in Figs. 6, 


7, and 8. 


FLOW CHARACTERISTICS IN THE YAWED CONICAL FLoy 
FIELD WITH WIND AxIS COORDINATES (7, 6,, ¢) 


(1) Method A 


In this method it is assumed that all intermediate 
surfaces remain circular cones and are rotated an 
amount «. The equation for these surfaces therefore 
is 
6. =0+ ecos¢g — 1/2€ cot Osin? ¢ (for@, <0<6 

(5 
Carrying out the operation indicated by Eq. (3) and 
using Eq. (5) and the first line of Eq. (2), for an example, 
the radial velocity component in the yawed conical flow 


field is given by 
u,(0,) = u(@) + u’(O)Ad + u”"(6)(Ab2/2) + 


= n+ ex cos g + e7(uM% + Us COS 29g) + 
[a’ + ex’ cos g + €(u’ + m2’ cos 2¢)] X 
[e cos g — hi 2)¢* cot 6 sin? 6| + 
[a” + ex” cos g + &(uy” + ue” cos 2¢)] X 
[(1/2)e? cos? go +... 


Upon expansion and collection of terms it becomes, to 
the order e’, 


u(0,) =u + (x + a’ )ecos g + 
[uo — (1/4)a’ cot 6 + (1/2)x’ + (1/4)a"Je + 
[uw + (1/4)a’ cot 6 + (1/2)x’ + (1/4)a”)e cos 2¢ 


This equation for u, would hold for the flow field from 
the surface up to the shock wave. 

For the radial velocity immediately behind the shock 
wave, Eq. (4b) is employed, yielding 


Ury (Or) = Uu(B,) + u’(0,,)A6,, + u”(6,,)(A8,2/2) +... 


= u, + eX, COS g + €7(M, + U2, Cos 2¢) + 

[u%,’ + ex,’ cos g + €2(u’y,, + U2,’ cos 2¢)] X 

[ea cos g + €°(By + Bo cos 2y)|] + 

UH,” + ex,” cos g + €(Uo,,” + U2,"Cos 2¢)] X 
[(1/2)€?a? cos? g + 
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Again, after expanding and collecting terms, 


(Ore) = Wy + (Xu + ati,,’)e cos g + 
Fru 
Uy + (1 2)ax,,’ + Bothy’ + (1 4)an,,”)e + 
[u2, + CU 2)ax,’ + Beit,’ + (1/4)a*H,,” Je cos 2¢ 


In a similar manner all the necessary flow relations in 
the conical field at angle of attack can be obtained. The 
primed quantities in the equations are derivatives with 
respect to the variable @ and can be found in explicit 
form in the latter part of this section. Since the equa- 
tions of the flow properties at the solid cone surface 
simplify when the boundary conditions are inserted, 
this special set of equations is included with those of the 
general flow field and at the shock wave as follows: 


At the yawed solid cone surface 6,,, 


+ ex, cos ¢ + €7?[u, + (1/2)a,] + 
e* cos 2g|[u, + (1/2)a 


u = Us 


= ¢€[m%, — x, — (1/2)z, ese 6,] + 


€ cos 2¢[v2, — a, cot 0, — x, — (1/2)z; ese 6, | 
Ww, = €,sin g + ée sin 2g 
[w, — ty esc 6, — (1/2)z, cot 8, 
n 1h. ¥r 
b, = ps|1 + €cos ¢ + 2 (— +—] + ; (6) 
a P/ s p pad 
Dy n* 
€’ cos 2 (2 + be ) 
p 237) . 
e 
7 , [p0 ul 
i, = B I + ecos ¢ (*) +e( 4 = + 
Pp’ s p 28°) 4 
» 2 u* 
e’ cos 2¢ (° + 5) 
p 26°) 
At the yawed shock wave surface 6,,,, 
Nin ti, + € cos g(x + av), + 
[uo + (1/2)ay + Bo + (1/4)a°*H" |, + 
€? cos 2g[u2 + (1/2)av + Bw + (1/4) aa"), 
Uy = vV, + €cos ¢(y + at”), + 
€ [vy + (1/2)ay’ + Bit” + (1/4a?n’”’|, + 
e cos 2g[v2 + (1/2)ay’ + Bou” + (1/4)a°*n’”’], 
Wry = sin g + e& sin 2¢[we + (1/2)az’]y 
/ hb! 
Pro = Pu [ + €cos ¢ (” + at ) 4. 
p Plu 
/ y (4) 
(Po an p a? p' 
( + fh — + = 
p 2p p : pe. 


p w 
(ow , a8 B’ a? p” 
¢ ar + Bo a 4 ee t 
p 26 p 1 plu 
, a Dp’ a? p” 
€ cos 2¢ (° e- > eg ESE | 
p 2p p +p 
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At the yawed intermediate conical surfaces 6,, < 6, < @,, 


uy = u + €cos g(x + 2) + 

e [uy — (1/4)8 cot 0+ (1/2)y + (1/4)a"] 4 
e& cos 2¢/u. + (1/4)3 cot 6 + (1 2)y + (1/4) nu” | 
v, = t+ ecos g(y + a’’) + 
e*[v — (1/4)a’’ cot 6 + (1/2)y’ + (1/4)a’"’] 4 


1) a’’ cot 6 4 


yee 


€e*> cos 2gjv2 + (1 
(1/2)y + (1/4)a 


es sin g + e& sin 29 [we +(1/2)s’ 


= pl + coos e ("4") 4 
Pp Pp 


Do , ” 
2 (é _ lp cot 6 4 ln 4 l p ) , 
p 4p 2p ap 


e’ cos 2¢ ( + Aa cot 6 + iy + LP ) 
p t p 2 p 1 p 


p 1 p 2 Lp 
ee Ta, wee: Pe 
€- CoS 2¢ os cot 6 + + 
pb 49 2p 49 


It must be remembered that the barred quantities and 
first-and second-order coefficients are tabulated asa func- 
* Thus @ is 
employed as a parameter relating the flow properties 
given by Eqs. (6), (7), and (S) to the spatial positions, 


tion of 6, the nonyaw conical variable.! 


in the yawed conical flow field as given by Eqs. (4a), 
(4b), or (5). A typical calculation scheme is given in 


Table 1. 


(2) Method B 


This method assumes that the intermediate conical 
surfaces are distorted linearly in shape from that of the 
solid cone surface to that of the shock wave at angle of 
attack. It is evident that Eqs. (4a) and (4b) have the 
same form, and therefore a linear variation between the 
two can be achieved by linearly varying the coefficients 


of the g terms. That is, 


6, = 6 + ek, cos ¢ + e*(ky + k; cos 2¢) (for 6,, < 0, < B,,,) 
(9a) 


where 


§ — @. 
ki =1-(1- a) | | 
6, — 6 
l l ‘ 6—06 
ki = —- cot @, + ( cot 6, + as) (9b) 
| | 16 6 
l ‘ 6— 06 
l cot #@. — ( cot #, — s:)) | 
| \4 


k; = 
i-~—§ 


The results of carrying out the operation indicated in 
Eq. (3) and using Eq. (9a) are 
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TABLE | i 
Calculation of Static Pressure Ratio p,/p, in ¢ = 0° Plane (Method A 
6. = 15 
M, = 1.9541 
e = -12 
54 
1 2 3 4 5 6 7 
@ (deg. ) Ag (rad. ) A9 (deg. ) 6, (deg. ) a? u a 
? : . 180 : . . ss 
See Eq. (5) 2 1+ 323 Ref. 1 Ref. 1 Ref. | 
7 
15 —(). 20944 —12 3 0. 12960 0. 60000 0 
17 —(). 20944 —12 5 0.12945 0.59930 —(). 03954 
l 
4 
8 9 10 11 12 13 14 9 
v? cot @ v cot 6 uv + vcot 9 a? — 7? u+ nu” ese @ ' 
{ 
11 
7? cot 1 7x9 6+ 10 5-8 : 12 ese 1 - 
0 3.¢ae. 0 0.6000 0. 12960 0.6000 3.8637 
0.00156 3.2709 —(0. 1293 0.4700 0.12789 —0.4757 3.4203 
P 
15 16 17 18 19 20 21 P 
esc? 6 esc? 6 —27 ii” —ii" cot 0 16+ 17+ 19 5 X 20 
14? xis —2 x 7 13 — 6 —9 x 18 
14.928 0 0 —1.2000 1.479 4.479 0.5805 0 
11.698 —(). 4625 —(). 07908 —1.0750 3.516 $. 138 0.4056 U 
22 23 24 25 26 27 28 
p ae v3 , 21 + 22 + 23 K 
20u"(u + u") (y — 1 ) 52 (i + u 24 X 25 26 u M, 
. 27 
—13 X 17 X 18 ee 0.405 13? 12 
0 0 0 0.3600 0 0.5805 $.479 
—(0 04044 —(0, 0000617 —0. 000193 0. 2263 0 0.3651 2.855 
29 30 31 32 33 34 35 
= P = = (a+ nu") es 32 25 ul “= Lb ul 
@° p a* a* @* v, 
05 29 x 8 
ee 7X 13X29 -—13 x 18 Xx 29 ; 7 + 28 -7 X 29 x 34 
10.84 0 —7 .805 0 0 $479 0 
10.85 —().2041 —5.548 0.1308 0.0296 2.815 1.208 
36 37 38 39 40 41 42 
P y —vcot é 2 = 2 y u+ u % : 4 
p . a- a- a- 
‘ —17 X 37 7 8 w, 
31 + 33 + 35 Ref. 2 —) X 37 5 13 39 Ref. 2 = 
—7.805 1.2000 —4.479 0 0 —(). 2503 0 
—4.310 0. 98383 —3.216 —(). 6007 0.2858 —(). 2124 0.01205 
43 44 45 46 47 48 49 
2--> —x (2 _ =) —scse 0 ux vy (ux + vy 
a* a- sd 
2 — 42 —41 X 43 Ref. 2 —14 x 45 6 X 41 7 X 37 47 + 48 fc 
2.000 0.5006 0.9803 —3.788 —(). 1502 0 —(0. 1502 : 
1.988 0.4223 0.7352 —2.515 —(. 1273 —(). 03888 —(). 1662 u 
T 
50 51 52 53 54 55 56 
v 38 + 40 + 44 + ; n y — 1) x + y’ 
om 9 x x ; ial —_ 
(¥y 1) as 13 49 x 50 “6-4. 53 y > a2 49 a4 a 
42 5 x 52 49 41 + 53 
0.405 Ref. 2 —(0.405 
5 12 . 5 13 
0 0 —7.766 —7.766 —1.6055 0.4694 —13.361 
0.03770 0.0030 —5.020 —5.081 —1.7813 0.5200 —11.128 
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i TABLE 1 
Calculation of Static Pressure Ratio p,/p; in ¢ = 0° Plane (Method A) (Continued from page 114 
. 57 58 59 690 61 62 
54 + 55 + 56 —30 X 57 -Wa+a) ! (@ + r) ¢cos ¢ 61 
a°* P P P 
—13 X 29 X 37 58 + 59 30 — 54 —0. 20944 X 61 
—14.497 0 7.805 7.805 1.6055 —0.3633 
—12.389 —2.529 5.075 2.546 1.5772 —(0.3303 
64 65 66 67 68 69 
p’ 1 7’ 1 p” 0 
LP vot 8 aS 63 + 65 — 64 «66 a. 64 + 65 + 68 
4p =~ /p + p Pp 
4 ] 
° * = 5 60 + | 36 0.04386 X 66 Ref. 3 
0 1.9512 —3.119 —0. 1368 3.370 5.321 
—(). 1669 0.1955 —2.317 —().1016 4.419 1.448 
71 72 73 74 75 76 
pr Ds - 937/(y — 1) p p 
Dp pi a,* pi Pi 
+ : 5 
 Sallaal Ref. | o ine 739.169 2X7 71X75 
0.7603 1.5488 1.0000 1.0000 1.5488 1.1776 
0.7632 1.5488 0.9988 0.9959 1.5424 1.1772 
up = ut + €cos g(x + kw) + P bh’ 
¥ ; br = pll+ecos¢ 3 + ky t + 
P | ae » p 
e( um + kd +—kiy + ea") + hg ; 
2 s » ( Po a OS Be 
\ | i + ky — + 5 : +7 he! 
“ 7 D y) Zz ) ‘ 
€’ cos 2¢( u. + ky +—-kyt bi") / / : t 
2 4 a De b 1 n’ ] 
e’ cos 2¢ (? + hs! + 7 +—k 
p p 2p 1 


The derivatives with respect to § are obtained in explicit form from Eqs. (5)—(12) in reference 
follows :* 

u = Vv 

u" = [a2(2a4 + v cot 6) — av*]/(a? — v?) 

[ Eq. (11) continued on page 116) 


* Eqs. (11) can be reduced at the body surface by means of the boundary conditions [reference 2, Eq. (13)]. 
v, = 0, y, = 2Uu, 
In the evaluation of sz’ at the body surface, use must be made of L’Hospital’s rule, giving 
’ 


s, = —2u, csc 0, — z, cos 6, 





63 
Po 
p 
Ref, 3 
—5.070 


—2.679 

70 
€? cos 2¢ 69 
0.04386 * 69 


0.2334 


0.1951 


) 
3 be 
p 


: A) 
b/1\ a 


Q) 


Cont. 


| “x ee ee 
%? = 7 + e cos ¢(y 7 kyii'’) + { (10) yy = @ k + € cos ¢ (: a. ky r) fe 
l l =/ / = J 
2 a oti’! Ly ,/ L 2 —sis ) , ] E l 
e(m that hy + hea”) + | Cee me) + 
7 Bp pb 2 p 4 p 
2 9 rida I k Pdi 1 act p2 p’ l z¢ l p” 
€° COS 29 | v2 + kgit’’ + 9° 4 kyu €’ cos 2¢ (2 + ks— + 2 — + -— ki? — | 
p p 26 4 Bp 
When 6, = 6,, and 6, = 6,,,, Eq. (10) reduces to Eqs. (6) 
w, = essing + e? sin 2g [we + (1/2) ki2’] | and (7), respectively. 


2, and are as 


(11) 
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—ftp l =9 2D 9 -" A Ld Pr) - 3 t= iw —\9 
u = — - | a*(v csc* @ — 2v — "cot 6) + Qu" (ua + un") + v? + (u4 + un")? 
a? — vy? a? 
y =x’ 
P a ss ow OP to = a y? a 
y = — -| — ycootd + 2— (4 + #”) ¥ (2 — —)— zgesc6 4+ (yy — 1) = (4 + 4”) (4x + vy) 
Go =: 7° a~ a’ a} : 
: ‘ . u ese 6 7d 
2 = —yecesc@—2 (cot or ) — E . : 
v v y(y — 1) 
p'/p= — (yv/a?)(a + i”) 
. n’’ _ Px , v ve 
/ =— oe (4 + u’") + (4+ u'')* Y (v + n’’’) (11) 
p a? a‘ a? Cont. 
P = — (v/a2)(a + &”") 
p 
p’’ x pa v epee rh a 
Fin i ae —~(4@+n' )—-—WUt+a )+@-y/ ) ita )2 
Dp o* a* a‘ 
/ - ~ , ~W 
u(u u) (y — 1) ,- % y v(u “) 
ps ~ = = | - Sa = (ue = DY) — T : a T 
p a? pb a? “+ nu” a? 
¢’ (a + 1”) FE (y—1).- x + y’ via + 1”) 
= = — = — (ux -- vy) — - : =2 
D a’ Dp a’ n+ a” a? 


CONE SURFACE CHARACTERISTICS WITH Bopy AXIS 
COORDINATES (7, (), W) 


For the purpose of obtaining and integrating surface 
distributions, it is more convenient to obtain the cone 
surface characteristics in body axis coordinates. Using 
Eq. (4a) for the body surface, the unit surface vectors 


are, to the order e’, 


n = [0,1 — (1/2)? csc? 0, sin? g, € csc 8, sin ¢] 

t; = (1, 0, 0) 

(12) 
to = (ty X n) = 10, — ecsce @, sin g, 1 — 


(e7/2) csc? 0, sin? ¢] 


And so, the velocity components normal and tangential 
to the body surface become 


Ung tye Vig (teres Urey Ws) | 
Veg = N° Veq( trey Urey Wrs) (13) 
We, = toe Vr.( tres Urss gs 


The directions of “x,, p,, and wp, are those given by an 
increasing spherical coordinate system with reference to 
the body axis, as shown in Fig. 2. Substituting Eqs. 
(6) and (12) into Eqs. (13), one obtains 


Up, = Us + ex, CoS g + €7[U0, + 
(1/2)#,] + & cos 2¢[uo, + (1/2)H,] 


Up, = €7(%0, — X,) + (14) 
e? cos 2¢(vo, — u, cot 0, — x, — 2, ese 8.) 
Wa, = €2, Sin g + e? sin 2g X 


[we, — a, csc 6, — (1/2)z, cot 6,] 


Now, one of the relations between the coordinate svs- 


tems (7, 6,, ¢) and (r, @, W) is 


cot ¢ = cos e cot y + sin e cot @ csc wy (15) 


From this, the following relations can be obtained: 


ecos ¢ = ecosy + (1/2)e? cot ® — 
(1/2)? cot @ cos 2y 
esing = esin y — (1/2)e? cot @ sin 2y (16) 
e’ sin 2g = e€’ cos 2y 
e sin 2g = e€ sin 2y 








Body oxis 


Wind oxis 


"Wigtees surface 





Fig. 2— Velocity components, Eq. (17), in body 
axis coordinate system (r, @, ,w) 
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(®, = 


Substituting Eqs. (16) into Eqs. (14) and since 
the relations for the surface velocity components 


are 
tn = ts t+ x, cosy + 
e?[uo, + (1/2)x, cot 8, + 1/2)u,] + 
e? cos 2W[u2, — (1/2)x, cot 6, + (1/2)%,] 
‘i, = €7(Uo, — Xs) + €? cos 2y X (17) 
(vo, — u, cot 6, —X; — Zs CSC 6,) = 0 
we, = €, sin py + 
e? sin 2Y(w., — ts csc 0, — 2, cot 6s) 


In Eq. (17), since vg, is the normal velocity to the sur- 
that this is satisfied can be 


face, Up, Should equal zero; 
(21) in reference 3. 


seen from Eqs. (31), (22), and 

In the case of the pressure and density on the solid 
cone surface, only a transformation of variables is 
Therefore, using Eqs. (6) and (16) and re- 
membering that @, = @,, the surface pressure and den- 
sity relations in a body axis coordinate system are 


be, = Ps [ + ecos yp (3). 
p 


necessary. 


p 2a? 2b P 
) vu l 
€’ cos 2 hs — ” cot i) | 
p 24? 2p 
> (18) 
pa; = Ds [ + ecosy (‘) a 
a’, 
) The . z 
2 “U) + e + cot i) aa 
p 2a" 2p : 


NORMAL AND AXIAL PRESSURE FORCE 
COEFFICIENTS 
Referring to Fig. 3, the axial and normal forces on the 
surface element dS are 


dFy = (pr, — pi)x tan? ®, dx dy 


and 
dFy = —(pp, — pi)x tan ®, cos p dx dy 
Introducing the pressure relation from Eq. (18) and 


integrating gives 


Fy = xi? tan? ®, X 


J ' Po . = yu" | 
et + < (2+ 5 root b+ 2) | — Pre 


and 
2)/? tan 0,p,€(n/p), 


Fy = (tT 


Thus, the pressure coefficients are given by 
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Fig. 3 — Definition of cone body geometry 


- {Fy Ps Pi 
Cx = - 
rD*q, 7/2M,? 


° Po yu 2 l 
pita (M4! 7 cot 6 + =) |- (19) 
p 2p 28") 7/2M;? 


and 


’ 1Fy (p./ pi) (n/p), cot 8 
Cy = SEX = — £ Pes PnP) CHE — (ag) 
rD*q 2 7/2M;? 
The drag and lift coefficient to the order e? can be ob- 
tained as follows (where positive angle of attack is 
=F: 


> 


Cp = Cxll — (2 2)] — Cy 
CL a= Cy[1 —_ (e? 2)] + eCx 


Employing Eqs. (19) and (20) gives 


Cp = P/Pr | | + 4« (Poy t+" cot 6 — ") — 
7/2M,? >p 2a p 2/, 
[1 — (&?/2)] (21) 
2M? 
and 
C, _ €/2(p/ pr)s(n/ p), cot 8, e[(p,/pr) — 1] (29) 


7 2M,? ¥Y 2M;,? 


The above aerodynamic coefficients are related to 
Kopal’s parameters (Ay, Ap, and Apyy) in the follow- 


ing manner: 


— (Se w)Ky 


yu.xd l 
—- + 
24,2(y — 1) 2 


Cy = 


Cy = AKp + 


Ilar - ae 
e (Ako + = )| Korn 
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Fig.4—Surface static pressure variation around cone at M=200 
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Fig.6 — Flow field velocity -methods "A" and "B" (second order) 


Ci, = —e« [(8 T)Ky — AKp| (2: 


Cp = A Kp +e 


: : yux,d : z | 

e) Kpyx — — - (4%, + :) + | 
24,7(y — 1) 7M, yM4 

where A = 2 if Kp is taken from reference | and 4 = 


8/x if Kp is taken from references 2 and 3. The term 
yuxd/2a,*(y — 1) appears in the above equation be. 
cause Eqs. (41) and (42) in references 3 are incom. 
patible. However, the tabulated values correspond to 
Eq. (42). Thus, in order to use the tabulated values of 


Kpyu, Eq. (41) should be corrected from 


Kp = Kp,|{1 + e’Apyy a 


: . c 1X 
Kp = Kp, | 1+ € (Kova - YU ,X 6 ) 4 
2a,*(y — 1) 


The effect of this term on the drag coefficient is small at 
moderate Mach Number but increases with increasing 
Mach Number as shown in Table 2. 


TABLE 2 
Calculation of Drag Coefficients, Eq. (23) 


Cot Co*) , 

6s MM, € Co* Cot ( Cp* 

roe 5 0.03722 0.03724 0.05 

£6602 10 0.07903 0.07911 0.10 

5 13.196 5 0.03116 0.03797 21.85 

eee 10 0.07147 0.09871 38.11 

- {5 0.09476 0.09542 0.70 

sa sha la 10 0.14217 0.14479 184 

ox . nie 5 0.47519 0.47682 0.34 

= 2.0665 10 0.49929 0.50585 1.3 
* Eq. (23) without the term yiisxsd/2452(4 ] 

T Eq. (23) 


CONCLUDING REMARKS 


The theoretical results for the pressure distribution 
around the cone and for the angle of attack variation in 
the axial force coefficient are compared with the ex- 
perimental data of Cronvich’ in Figs. 4 and 5. The in- 
dication is that the method presented for use with 
Kopal’s tabulations is satisfactory both in regard to 
magnitudes and trends of variation. In these compari- 
sons the theoretical results were linearly interpolated 
with Mach Number to correspond to the experimental 
Mach Number. 
mental values given in Fig. 5 are integrated surface 


It must be remarked that the experi- 


pressures from experimental values’ averaged for a 
series of test runs. In Figs. 6, 7, and 8 are presented 
the comparisons of flow field properties as described by 
methods A and B. It is believed that method B is su- 
perior for establishing the initial conditions for a charac- 
teristic solution of the flow field for a pointed body of 
revolution. The comparisons of the first- and second- 
order solutions for some flow properties are shown in 
Figs. 9 and 10. 

It should be mentioned that Stone’s solutions yield 
the following information on the entropy distribution in 
the yawed conical flow field: In the first-order analysis 
In the 


second-order analysis the entropy is constant only in the 


the entropy is constant in each meridian plane. 
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Fig.9 — Flow field velocity for first and second order theory 


top and bottom meridian planes. The true entropy dis- 
tribution should be such that the entropy along each 
stream line is constant. This dictates the method of 
distorting the intermediate conical surfaces (see method 
B above). Thus with the proper distortion of these 
surfaces, the second-order solution should approach the 
physical conditions. Method B is only a first attempt 
in this direction. When properly interpreted, Stone’s 
results for the second-order flow characteristics on the 
solid cone surface are not affected by the method of dis- 
tortion and should be exact to the order of analysis. In 


AROUND 
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view of the first-order results of constant entropy in 
each meridian plane, Ferri’ has introduced a vortical 
layer through which the entropy can vary and has de- 
veloped a method of using the first-order results of 
Stone-Kopal.?. A comparison of the experimental and 
analytical results given by Ferri’ for the pressure co- 
efficients over a 10° cone at M = 6.86 are made with 
Stone-Kopal’s second-order values in Figs. 11 and 12. 

The authors are indebted to Dr. M. D. Van Dyke for 
the many constructive discussions and to Miss B. 
Dodge for computational assistance. 
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Heat Transfer by Laminar Flow from a 
Rotating Plate’ 


KNOX MILLSAPS? anp KARL POHLHAUSEN? 
Wright Air Development Center 


ABSTRACI 


An exact solution of the heat-transfer problem for the von 
Karman example of the laminar flow of a viscous fluid over a 
rotating plate is given in dimensionless form and physically dis- 
The solution is explicitly given for a constant tempera- 
The nu- 


cussed. 
ture on the plate with viscous dissipation included. 
merical results are given for Prandtl Numbers from (0.5 to 10. 


INTRODUCTION 


6 Bees ARE ONLY A FEW CASES KNOWN where the 
mathematical structure of the flow of a viscous 
fluid can be rigorously determined without simplify- 
ing the basic hydrodynamic equations by neglecting 
the effects of certain terms in the differential system 
describing the motion. One of these few solutions for 
the case of the laminar flow of a fluid over a plate rotat- 
ing at a constant velocity is due to von Karman.' It 
is the purpose of this paper to show that the associated 
heat-transfer problem for the von Karman example can 
also be solved exactly.**: ” 

It should be noted that the physical significance of 
these exact solutions for the rotating plate lies not in 
their engineering importance, since in practical cases 
the motion rapidly becomes turbulent,t} * but rather 
in the possibilities of making critical comparisons with 
approximate methods. 


(I) Tue Basic EQuaTIONS 


The equation of continuity, the Navier-Stokes equa- 
tions, and the energy equation for the steady motion of 
an incompressible fluid (density p is constant) when ex- 
pressed in cylindrical coordinates (7, 3, z) with velocity 


Received August 13, 1951. 

* The authors wish to express their thanks to Frank W. Bubb, 
Chief Scientist, Flight Research Laboratory, for making this in- 
vestigation possible and to Mrs. Olga Kay Francis for her help 
in performing the tedious numerical calculations. 

Tt Now, Professor of Physics, Alabama Polytechnic Institute 

t Research Scientist, Applied Mechanics Research Group. 

** Two further cases of fluid flow, together with the associated 
heat-transfer problem that can be solved exactly, are the Hamel 
Jeffery examples of flow in convergent and divergent channels 
The explicit solution of the hydrodynamic problem in terms of 
elliptic functions was given by Jeffery? and Hamel.? The solu- 
tion of the energy equation in terms of elliptic functions arising 
from a generalized inhomogeneous Lamé equation will be given 
in a forthcoming paper by the present authors. 

tt A summary of the results from experimental investigations 
is given in Goldstein.* 


components (“,, Uy, uz) and kinematic viscosity vy are, 
respectively, ff 4 
Ou ‘ u“, . 10, Ou, 0 
+ = (1.1 
or r r Ov Oz 
Ou, 4 Uy OU; Un” i Ou, 1 Op \ 
Uu; - = Uu; = — + 
or r Ov r Oz por 
O°u, és 1 Ou, u, 4 1 O7u, 2 Ou, | rad 
v = a i 
or? r Or y? y* o- r? Ov 02? 
(1.2 
Ou» fi Uy OUy 4 U,Uy 4 Ou y 1 Op 
Uu, - u- = — 
or r Ov r Oz pr Ov 
O75 1 Ou, Uy l O°u, 2 Ou, O"u, 
v > oF se leer as a + 
or’ r Or ‘he r* oOv- r? OW 02? 
(1.3 
Ou; ‘ Uy OU; 1 Ou. 1 Op 4 
Uu, Uu; = — 
or r Ov ~ OZ p Oz 


= 4 | Ou; 
. or r oO 


and, denoting the temperature of the fluid by /, 


Ol Uy Ot Ot 
oe (1 + + Uu, = 


1 O7u- O7u ) 1“ 
r? Ov? Oz? 


Or r Ov " OZ 
k (= 4 1 Of 4 1 OF i: ~*) rn 
p \Or? r or r> Ov? 02? 


v 26,5" + 2e 59° = 2622” ae é, ” + C52" + tes") (1.5 


where c, is the specific heat expressed in mechanical 
mass units, k is the thermal conductivity, and 


Ou, 1 OUs u, 
e;, = ; ess = 
or r Ov r 
Ou. | Ou. Ou » 
ézz = ’ Cs: = + (1.6) 
Oz r Ov Oz 
Ou, 4 Ou, 1 Ou, 4 Ou, Uy 
e.. = P e+9 = ates 
Oz or r ow or r 


As shown in the schematic diagram in Fig. 1, the 
problem of an infinite plate being rotated about a point 
in its plane at a constant angular velocity w and being 
heated to a constant temperature /,, with the fluid hav- 


tt These forms of the continuity and dynamical equations 
. ‘4 -r . e ¢. 
are given by Hopf.‘ This form of energy equation may be in- 


ferred from Goldstein.‘ 
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Fic. 1. Schematic diagram showing coordinate system and 
velocity components. 


ing an ambient temperature /,, is characterized by the 


boundary conditions 


t(r, 3, O) = ft, tir, 3, o) = t, 
u(r, 3,0) =O0 u(r, 3, o) =O | (1.7) 
u(r, 0, ©) =O u(r, 3, O) = O 

U(r, 3, O) = rw 


Since the angular velocity is assumed to be constant, 
the physical problem is axially symmetric, and, con- 
sequently, the partial derivatives with respect to J 
vanish. Accordingly, the partial differential system 
Eqs. (1.1)-(1.5)] may be reduced for laminar flow to a 
set of ordinary nonlinear differential equations by as- 


suming 
u, = rf(z), My = rg(s), uz = h(z) 
for the velocity components, 
p = p(s) (1.8) 


for the pressure, and 


t = r’s(z) + q(z) +1 


for the thermal distribution. 
The reduced equations are 


27 +h’ =0 (1.9) 


f—-g thf’ = of" (1.10) 
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2fg + hg’ = vg" (1.11) 

hh’ = —(p’/p) + vh" (1-32) 

2c,fs + Cohs’ = (Rk/p)s” + v(f’? + g’?) (1.13) 
cyhq’ = (k/p) (4s + q”) + 12vf? (1.14) 


where the primes denote differentiation with respect 


to 2. 
The equations may be made dimensionless by intro- 


ducing the following variables: 


f=wF, g= 0G, h = (w)'*H 


Ww Vw 


b = pwP, s=-S, q= QO 


Se 
r=(’) R, 


where o is a modified Prandtl Number in which the spe 
cific heat at constant volume has been used instead of 
the specific heat at constant pressure. (In the case of 
incompressible fluids the distinction vanishes.) The 
dimensionless equations resulting from the application 


Lal 
|| 
a 
eis 
—” 
Str 
Q 
| 
.-) 
xr s 


t = (pw/c,)T 


of these transformations are 


H’ +2F=0 (1.16) 

F” — HF’ — F? = -—G@ (1.17) 

G” — HG’ + 2FG = 0 (1.18) 

P’ — HH" + Hi’ =0 (1.19) 

S” — oHS’ + oH'S = —o(F’? + G") (1.20) 
Q” — oHQ’ = —(4S + 12¢F?) (1.21) 


where the primes now denote differentiation with 
respect to £& The boundary conditions for the ve- 
locity components have been transformed into 


F(0) = 0, F(a) = 0 (1.22) 
G(O) = 1, G(~) = 0 (1.23) 
H(0) = 0 (1.24) 


One obtains the thermal distribution by combining 
Sand Q into 


T=RS+Q04+T (1.25) 
with the boundary conditions becoming 


x0) = 0, S(o) = 0 (1.26) 


Q(0) = T, — T Q(o) = 0 20 


(II) SOLUTION OF THE EQUATIONS 


The first four equations [Eqs. (1.16)—(1.19)] with 


the boundary conditions |Eqs. (1.22)—(1.24)| consti- 
tute the hydrodynamic part of the problem and have 
been integrated numerically by Cochran,® who used a 


process of continuous joining of power expansions from 
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ae I 7 a ‘ 40 
[| | | layer hypothesis is made, one finds a boundary laye, 
‘ |_| SgISaee ah _| of approximately constant thickness over the entire ‘ 
| | -H(eo)= 0.886 | plate, provided the small influence of flow toward the 
s Lo, : 4 { plate is neglected. This interesting conclusion, firg 
ea 2 ae 30 
noted by von Karman, may be drawn from the folloy. 
7 e . el rir . ° 
ing elementary considerations: The magnitude of the 
° e ry — lifes ° 25+ 
. fluid velocity V(R, £) divided by (8w) “* is given by 
ig | 
, (R?F? + R°G? + H?)” (2.] 20 
and for large R this expression becomes 
“ 15 
x 1971/9 
R[F? + G?]”* (2.2 
3 
——" 10 
Eq. (2.1) clearly shows that a boundary layer consistent 


with the Prandtl hypothesis does exist in the von Kar. 
, " a< “— . 5 
man example, and Eq. (2.2) shows that it is one of ap- . 


proximately constant thickness, since in almost any 





F ° a ° 
particular physical example (w/v) ’* is quite large, o- 


1 2 3 4 & For a physically significant definition of the thickness 


Fic. 2. Numerical values of the velocity components calculated of the boundary layer 6, an equalization of areas in the 
by Cochran. ; 


velocity profile is perhaps to be recommended, and, in 
symbolism, this definition is 

> > | ‘ NW0 
/ V(R, &)| — |V(R, ~)|| dé Fig. 
JS0 | func 


> a: = oul 
max. V(R, &)| — | V(R, ©) _ 
O<s¢:« | | Sx 


tem 
plat 
Sat 
7 fyin 
Se en ee ee eee diffe 


two 


Fig. 4 is an application of the definition to the case of 
the rotating plate with flow toward the plate neglected 
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Fic. 4. Approximate velocity profile showing thickness of the 
boundary layer. by 





Fic. 3. Three-dimensional model of the flow over a rotating 
plate. 


the origin with asymptotic expansions from infinity. 


P= 1400) 
The numerical results of Cochran are given in Fig. 2. 





The three-dimensional composition of the velocity 
components shows that a spiral fanlike motion of the 
fluid occurs in the von Karman example. <A _photo- 5 
graph of a three-dimensional model of this fluid motion 
is given in Fig. 3. 
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One should note that the spiral motion shown in ee F- e e & Se S UF i 

sa ats 5 RZ) 
Fig. 3 is an exact solution of the hydrodynamic equa- en ) 
. oe ° ‘ Ww 
tions, and, if a comparison with the Prandtl boundary- Fic. 5. Thickness of the boundary layer as function of R. 
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Numerical values of S with Prandtl Number as a param- 
eter. 


Fig. 5 shows the thickness of the boundary layer 6 as a 
function of R with flow toward the plate taken into ac- 
count. 

Somewhat similar to the hydrodynamic problem, the 
temperature field in the fluid due to a heated rotating 
plate is characterized by the dimensionless functions 
Sand Q which are solutions of Eqs. (1.20)—(1.21) satis- 
conditions (1.26)—(1.27). The 
.20)—(1.21)] is composed of 


fying the boundary 
differential system [Eqs. ( 
two inhomogeneous linear equations of the 
order, with variable coefficients that are functions of 
the hydrodynamic variables and the Prandtl Number 


second 


6. 

An analytical solution of the equation for S(&) is un- 
known; hence, it appears necessary to obtain solu- 
tions for parametric values of o by numerical methods. 
The straightforward method of obtaining the required 
numerical solution would be to assume initial values of 
S’(é) and to integrate the differential equation for S(é) 
by some numerical method and by subsequent inter- 
polations to obtain a correct value for S’(0) such that 
the boundary condition at infinity is satisfied. How- 


ever, since the equation is linear and one of the second 


B= 


Therefore, the required numerical solution is 


"'* 3 a 
S(~) = 91(€) y | M(é)y, dé dé — 
| o M1 JO 


res 


where \/(£) is defined in E 
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"oY we ee ee 
_ = M(&)y, dé dé / - 
0 V1 JO / 0 v1 
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order, the following simple method may be applied 
without the necessity of trial integrations :° 

The first-order term, —o//S’, in Eq. (1.20) may be 
eliminated by the substitution 


y = S(t) exp [—(0/2) fH dé] (2.4) 
The resulting normal equation is 
2 30H’ oT” - 
ee cee Ss iin 
a il 
—o(F’? + G’*) exp (- = / i] it) = M(t) (2.5) 


If y,(€) is a numerical solution of the homogeneous 
normal equation for arbitrary nonvanishing initial 
values of the function and its first derivative, which is 
obtained by the Runge-Kutta method (for instance), a 
second fundamental integral is given by 


(2.6) 


yo() = 


and the particular integral is given by 


— ] “s Z 
vp(é) = VM f - M(&) y, dé dé (2.7) 
Ji Vl Ski 


Since superposition is a valid procedure, the general 
solution is given by 
y = Ay + Bro + yy (2.8) 


where A and B are arbitrary constants to be deter- 
mined by the boundary conditions. If the boundary 


conditions are 


y(é1) = a, v(é2) = b (2.9) 
then 
A = da ‘wi (1) (2.10) 
n(t —_ n(fo) — vy,(£o)v,(& 
B= byi(&1) ayi(&) Vp (Eo) 1 (1) (2.11) 


yi(&i ) Vo( Ee) 


For the case under consideration, it is seen from Eqs. 
(1.26) and (2.4) that Eqs. (2.9) reduce to 


y(O) = 0, y(o) = 0 (2.12) 


and the particular values of the constants are 


(2.13) 


9 


pm a l "a ! 
_ -f M(é)y, dé dé . 
bis eG fire) eas 


j. (2.5) and y,(£) is an arbitrary nonvanishing solution of the homogeneous normal 
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equation. Since it was necessary to assign a definite value to the upper limit of integration in the infinite jp. 
tegrals, this numerical method was checked for the indicated Prandtl Numbers by the conventional Runge-Kutty 
method using interpolation to determine the initial slope satisfying the boundary condition at infinity. Fig, ¢ 
shows the numerical results for parametric values of o up to 10. 

The equation for Q, 


Q” — oHQ’ = —(4S + 12cF?) = N(é) (2.15 


does not contain the dependent variable explicitly; hence, the equation may be reduced to a first-order linear 
equation for which the solution is well known. The resulting Q may then be obtained by an additional integra. 
tion. 

In terms of quadratures the solution of Eq. (2.15) is 


Q= y| [exp (o J Hadé) | | f N(é) exp (—o fH dé) dé + G| dé+ Cy (2.16 
0 0 
where C, and C; are integration constants to be determined from conditions (1.27) and are given by 
CG, = T,—T., (2.17 
C _ —[l2+ (Ty, - Fatt I, (2.18 
where 
Ii, = | exp (o JH dé) dé | 
0 
x "E ) 
2 [ exp (o fH dé) | N(&) exp (—o fH dé) dé dé 
0 


J/J0 


(2.19 


~ 
I 


It is then seen that the solution for Q may be written as 


| : exp (o fH dé) dé I “ 
Q=[T.—T.a]J}1-— ] - A J exp (o fH dé) dé + 
1 1 0 


I exp (o fH dé) I N(£) exp (—o J'H dé) dé dé (2.20 
0 0 


Now, it is clear that Q may be divided into two parts, Combining Eqs. (1.25) and (2.21) one obtains the 
Q; and Q2, one being a factor times the temperature dimensionless thermal distribution 
difference 7, — 7. and the other being only a function - -_ - - = 

tha Ah tae ,  . T=RS+O(T,s —To) +Q+T. (2.22 
of the viscous dissipation and independent of the tem- 
perature difference. Accordingly, one has Replacing T by (c,/vw)t and R by (w/v) ‘’*r, one has the 


final result 


Q = Q(T, — Tu.) + Qe (2:21) 


Figs. 7 and 8 give the numerical results for Q; and Q, @ 
° ° 14 
for parametric values of o up to 10. 
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HEAT TRANSFER 


vo | wr? _ 
+ QOilto — to) + - | S + 0.| (2.23) 


Cy v 


t= le 


It is interesting to note that in the cases where vis 
cous dissipation can be neglected —i.e., N(~) = 0, 
Eq. (2.23) reduces to 


t= te + Olt — ta) 


where / is given in degrees centigrade and, of course, 
wr?/v is the Reynolds Number for a given value of r. 

The foregoing determination of the thermal distribu 
tion enables one to calculate the amount of heat that is 
transferred to the fluid. This calculation may be made 
by noting that the amount of heat transferred by con 
duction per unit area and unit time is given by 
—k(0t/Oz)--» and that the total heat flow per unit time 


from one side of a plate of radius 7 is 


° 2x "ro Ol 
¢ = k / | ( ) r dr dvd (2.24) 
< 0 Oz 2=—0 


By using transformations (1.15), Eq. (2.24) may be put 


into the dimensionless form 


a [ {~ (7) 
: pv wo = a JO /0 O& : 


‘ 


RdR dv 


0 


(2.25) 


and, after substituting from Eq. (2.22) and perform- 
ing the indicated integrations, one deduces that 


rR” 


oO 


Ro? . al _ , 
& S’(0) + (Ty — To)Qi’(0) + Qe | (2.26) 
L 2 
In Eq. (2.26), S’(0) may be calculated from Eq. (2.14), 
and Q,'(0) and Q2’(0) may be easily obtained from Eq. 
(2.20). Fig. 9 gives their numerical values as func- 
tions of o, and one sees that the dimensionless heat 
flow is inversely proportional to the Prandtl Number 
and approximately proportional to the Reynolds Num- 


ber. 
It will be shown later that in most practical cases the 














term containing 7, — 7°,, as a factor is extremely large 
when compared to the remaining two terms. Con- 
sequently, since Eq. (2.20) shows that 
3 | —~——_ 

| ao 
2 | | = 

Lo Sw 
0 } o 
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| Aw 
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Fic. 9, Initial slopes of .S, Q:, and Q» as functions of the Prandtl 


Number. 
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O,/(0) = —(1/h) (2:27) 
one has the approximations 
E = (#Ro?/al,) (Ty fo) (2.28) 
and, in dimensional form, 
e wka *ry"(ty ._/ 9" es (2.29) 
where /, is given by Eq. (2.19). 
(III) EXAMPLES 
A. Water 
A plate heated so that ¢,, = 35°C. rotates at 30 revolu 
tions per min. (w = 2(30)/60 = 3.14 rad. per sec.) in 
water in which ¢, = 30°C. It is assumed that for the 


temperature range under consideration the necessary 


physical constants have the following values: 


= 0.0075 sq.cm. per sec. 

0.00149 cal. per cm. sec. 
ae 

= 0.997 cal. per Gm. °C. 


v (kinematic viscosity) 
k (thermal conductivity) = 


c, (specific heat) 


p (density) = 0.996 Gm. per cu.cm. 


With these basic physical constants and the value of 
the mechanical equivalent of heat (4.184 xX 10’ ergs 
per cal.) it may be shown that 


0 = 5.0 
w/v = 42 & 10 per sq.cm. 
va (0.0075) (3.14) - 
= - — = 5.6 X 10-" °C, 
C (4.184 & 10°) (0.997) 


Also, since the critical Reynolds Number (Re), for 
transition from laminar to turbulent flow is approxi- 
mately 4.7 X 10** the critical radius r, is found by 


r- = [p(Re)./w] “” = 10 cm. 


The thermal distribution is given by Eq. (2.23) and, 
with the particular values of the present example, takes 


the form 


t= 30° + , 

(5.6 & 10~™)°[(42.X& 10)r2S + Qe] + 5°Q, 
in which S, Qi, and Q» for the parametric value ¢ = 5 
are given in Figs. 6, 7, and 8, respectively, and the value 


of ry ranges from 0 cm. to 10 cm. (critical radius). Since 
max.|(42 & 10)r°S] = 11 X 10%, the radial dependence 
of the thermal distributions is negligible. Hence, it is 
clearly seen’ |that 5°Q; is the decidingly predomi 
nant term in determining the temperature profiles and 


that 
t = 30° + 5°Q, 


For the heat flow from one side of the rotating plate of 


10 cm. radius, one obtains from Eq. (2.29) 
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who "rot — t,) 


e€ = = 


py], 
(3.14) (0.00149) (3.14)'*(10)? (5) 
(0.0075) '”* (1.172) 
cal. 


e 1] 


sec, 


B. Argon 


If one desires to observe the radial variations in the 
thermal distributions, the observable magnitudes of 
(te — ¢,,)Q: and (vw/c,) (r’w/v) S must be comparable 
and physically significant. A heated plate rotating in 
argon at low pressures appears to offer definite possi- 
bilities for making experimental investigations on the 
radial dependence of the temperature profiles. These 
possibilities can be attributed to the distinctive physi- 
cal properties of argon. 

If a plate heated so that ¢, = 
287 X 10 revolutions per min. (w = 27(287 X 10)/60 = 
300 rad. per sec.) in argon at 5.00 X 10~? atmospheres 
so that ¢, = —2°C., the pertinent physical constants 


have the following values :* 


2°C. is rotating at 


v = 2.57 sq.cm. per sec. 

k = 3.89 X 10~*° cal. per cm. sec. °C. 
C, = 7.45 X 10~* cal. per Gm. °C. 

p = 8.18 XK 10~° Gm. per cu.cm. 


and it may be shown that 


0.400 
117 per sq.cm. 
1 OM 


o 


W/V 


vw/Cy = 2.47 XK 10 


If the critical Reynolds Number is assumed to re 
main approximately the same for the proposed experi- 
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ment as in air at normal pressures, then 


v(Re).] (2.57) (4.7 * 104)]'” ; 
r= = = 20 . 
w 300 _ 
Therefore, since vw/c, max.|(wr?/v) S|] ~ 1.2°C. is clearl 


measurable and roughly of the same order of magni 
tude as 4°Q, it seems that measurements determining 
the radial dependence might be taken on plates rotating 
in argon at low pressures. Furthermore, the suggested 
experiment has additional features in that the turby 
lence level in the argon flowing toward the plate may be 
kept extremely low and in that, as has been previously 
recognized, the determinations of the critical Reynolds 
Number and the curve of neutral stability would ap- 
pear to have particular significance, since the velocity 
and temperature profiles of the laminar flow can be 


determined exactly. 
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An Experimental Investigation of the 
Isothermal Laminar Boundary Laver 
on a Porous Flat Plate’ 


PAUL A. LIBBY,t LAWRENCE KAUFMAN,}| ano R. PAUL HARRINGTON** 
Polytechnic Institute of Brooklyn 


SUMMARY 


The transition Reynolds Numbers and velocity profiles of the 
laminar isothermal boundary layer on a porous flat plate with 
suction or injection have been measured and compared with 
the laminar boundary-layer analyses that have been carried out 
in the past several The indicate that porous 
plates closely approximating the mathematical or ideal porous 
plates can be realized and that the theoretical velocity profiles 


years. results 


are in good agreement with experiment. 


SYMBOLS 


surface area of the porous plate in turbulence 
channel wall 

constants depending on the physical properties 
of the fluid medium and of the wire and on 


the distance from the wall 


A, A’, 8,38" 


E = voltage drop across the hot wire 

Ep = voltage drop across the hot wire when u 0 

E, = voltage drop across the hot wire when u = l’ 

Q = total volume of injected or withdrawn fluid 
flowing per unit time 

R, = Reynolds Number based on x 

R = Reynolds Number based on x; 

R.,, = transition Reynolds Number based on x 

u = local velocity in the x direction 

l = free-stream velocity in the x direction 


= local velocity in the y direction 
= velocity in the y direction at the wall, positive 
for injection, negative for suction 
4 = distance along the plate, positive downstream 
measured from leading edge of imaginary 
flat plate (see Fig. 8) 
x; = distance along the plate from 
edge of the imaginary plate to the leading 


the leading 


edge of the porous section (see Fig. 8) 
distance from and normal to plane of the 


plate 
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n = Blasius variable = (y/x)V R, 

m = coefficient of viscosity 

v = coefficient of kinematic viscosity 

£ = suction or injection parameter = (%/U’)?R, 

é; = suction or injection parameter for the partially 
porous plate = (v/U’)?R,, 

p = mass density of the fluid 


INTRODUCTION 


She PURPOSE OF THIS INVESTIGATION is (1) to obtain 
experimental data on the characteristics of the 
isothermal laminar boundary layer on a porous flat 
plate subjected to homogeneous suction or injection 
of the same fluid as that of the main stream and (2) 
to compare these characteristics with the predictions 
of the laminar boundary-layer theory. This represents 
the first of a series of experimental investigations of 
boundary-layer characteristics on porous walls, with 
and without heat transfer. 

The characteristics of the boundary layer on walls 
made of porous material with suction or injection are 
of considerable interest. Suction applied through such 
walls inhibits boundary-layer separation in the presence 
of adverse pressure gradients' and laminar boundary 
layer transition.? By injecting a coolant through a 
porous wall, a thick heat-insulating boundary layer 
can be established to protect the wall from high- 
temperature gases such as exist in combustion chambers 
and around turbine blades.’ 

In addition to the cited references, which present the 
results of experimental investigation on these uses of 
porous materials, several theoretical investigations of 
the characteristics of the laminar boundary layer 
over a flat plate with suction or injection have been 
carried out. Schlichting,‘ Iglisch,> and Schlichting 
and Bussmann,® have determined the characteristics 
of the incompressible boundary layer over a flat plate 
subjected to uniform suction and to suction or in- 
jection varying as 1/+/x, where x is the space variable 
Yuan’ has calculated the laminar 
boundary-layer velocity temperature charac- 
teristics for incompressible flow over a flat plate subject 
Yuan* and Lew et al.’~'' have 
the compressible 


along the plate. 
and 


to uniform injection. 
calculated the characteristics of 
laminar boundary layer with uniform suction and 
injection, with and without heat transfer. References 
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Fic. 1. Turbulence channel showing adjustable side wall and 
hot-wire probe 
12-15 give the characteristics of the compressible 


laminar boundary layer with partial injection and 
suction. 
the laminar 


theory that has been developed by many investigators 


Moreover, boundary-layer stability 
for the flow over impermeable walls has been applied 
without modification to the theoretical laminar bound- 
ary-layer velocity and temperature profiles presented 
Although this 
application of the stability theory has not been justified 


in some of the references cited above. 


rigorously, it has been assumed that the stability 
theory will give at least the qualitative influence of 
suction and injection on laminar boundary-layer 
stability. * Ulrich'® that 
greatly increases and injection decreases the stability 


Thus, has shown suction 
of the incompressible laminar boundary layer. Lees” 
Mach (M < 1) 


laminar boundary-layer stability with injection varying 


has shown that for small Numbers 
as 1/7 x depends on the opposing effects of the injec- 
tion, which is destabilizing, and on the cooling, which 


is stabilizing, and that for & < 0.25 the cooling pre 


dominates. Libby, Lew, and Romano! applied the 
stability theory to the profiles obtained by the integral 

Their results thus indicated, 
the 


methods in reference 11. 


at least qualitatively, simultaneous effects of 
* The authors have recently been informed by Prof. Lester 
Lees that this application can be rigorously justified for & on 


the order of unity. 
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compressibility, uniform suction or injection, and heat 
transfer. 


Despite the interest in the characteristics of the 
boundary layer on porous walls as indicated by the 
number of experimental and theoretical investigations 
in the cited references, there appears to have been no 
experimental investigation of the details of the actual 
boundary layer subjected to homogeneous suction and 
injection. There is need for such research; the 
theoretical analyses of the laminar flow should be sub 
stantiated by careful tests as the impermeable wall 
boundary-layer theory has been. In addition, the 
question arises as to whether extensive regions of 
laminar flow can exist with even small rates of injection 
over actual porous plates with their inherent roughness 
and with their large but finite number of pores. There 
exists the possibility that the injected fluid may issue 
from these pores in the form of discrete jets that would 
not permit laminar flow to exist. (That suction can 
maintain an indefinitely stable laminar layer has been 


clearly indicated by reference 2.) 


In addition, a detailed examination of turbulent 


boundary layers on porous walls, especially with 


injection, is necessary to establish a theory for such 
boundary layers. With the present state of knowledge 
regarding turbulence such a theory must be semi 
empirical. Measurement of the mean velocity profiles 
and of the skin friction can supply the required data. 


Turbulence channel showing injection bell and connec 
tion to flowmeter and supercharger 
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Flow diagram of turbulence channel 


This report presents the results of an experimental 
investigation of the characteristics of these boundary 
namely, iso- 
Investigations 
with 


layers under the simplest conditions 
thermal laminar flow over a flat plate. 
of the isothermal turbulent boundary 
injection and of the laminar and turbulent boundary 
layer with heat transfer and injection are under way 
and will be reported upon subsequently. 


layer 


EXPERIMENTAL EQUIPMENT AND MEASUREMENT 
TECHNIQUE 


The investigation was conducted in the 30- X 2-in. 
turbulence channel at PIBAL. A_ porous plate 

« X 14in. X 5 ft. was mounted flush with the inner 
surface along the centerline of one 30-in. wall imme- 
diately downstream of the inlet bell. Behind the 
porous plate was an injection air bell that was con- 
nected through a flowmeter to a supercharger. By 
changing the connections to the supercharger inlet 
and exhaust, either suction or injection through the 
porous plate could be applied by the supercharger. 
A 5-ft. section of the other 30-in. wall was hinged at 
the inlet bell so that the pressure gradient caused by 
skin friction and mass exchange through the porous 
plate could be effectively eliminated. Special care was 
exercised so that the pressure gradient from the leading 
edge of the porous plate to the measuring station 
was negligibly small, the deviations of the static pres- 
sure being made less than +1 per cent of the dynamic 
pressure. Although this care resulted in larger devia- 
tions downstream of the measuring station, the pressure 
drop properties of the porous plate were such that 
deviations of the injection or suction velocity from the 
average were negligibly small throughout the length of 
the porous plate. Thus, the assumptions of a flat 
plate with uniform suction or injection were satis- 
factorily achieved. Figs. 1, 2, and 3 show the two 
sides and a flow diagram of the channel. 


The turbulence factor o'® in the free stream at the 
inlet of the channel has been assumed sufficiently small 
so that the free-stream turbulence would have no effect 
on transition. This is believed justified because the 
maximum velocity at the inlet to the channel during 
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these tests was small, approximately 3 ft. per sec., 
and because the air flow through the channel was 
discharged outside the building and therefore not 
continuously circulated within the test room.* 
Because the results of these tests must be associated 
with the particular porous plate used or with porous 
plates closely related in porosity, permeability, and 
grain size, the characteristic of the plate will be dis- 
cussed in some detail. The plate was made from 
sintered bronze by the Wel-Met Company of Kent, 
Ohio. Fig. 4 shows a photomicrograph of the porous 
material in a section normal to the plane of the plate. 
The pores have been filled with plastic under pressure, 
and the surface has been heavily etched to make the 
pores and grains distinct. The pores are black and the 
bronze is white. 
facts of interest may be established by 
The microstructure of the porous 
No simple flow passages 


Several 
inspection of Fig. 4. 
plate used here is complex. 
are evident. Moreover, the 
such that no discrete jets could be expected from the 
The airflow through 


surface conditions are 
low pressure side of the plate. 
this plate would appear to be more of a percolation 
than a flow through small discrete channels. This 
characteristic implies that the porous plate used here 
provides a satisfactory approximation to the ideal or 
mathematically porous plate assumed in the theory 
of the laminar boundary layer with injection—namely, 
one for which the injection velocity v is either a 
constant or a smooth function of distance along the 
plate. This implication has been verified experi- 
mentally by placing a hot-wire at several wire diam- 
eters from the wall with only injection air flowing 
through the plates (U = 0, % > 0). The velocity 
indicated by the hot-wire was steady and independent 
of small changes in the position of the wire in a plane 
parallel to the porous wall. Moreover, during measure- 
ments of the velocity profiles, the fluctuation of the 
wire with the laminar boundary-layer 
velocity oscillations decreased as the wall was ap- 
Near the wall the indications were com- 


associated 
proached. 
paratively steady. 


the time of these tests, a suitable instrument for the 
However, 


* 
qualitative determination of o 4vas not available. 
one is under construction and will be used to check this assump- 


tion in the near future. 





Photomicrograph of a section of the porous plate 
The injection or suction flow is parallel to the short side of the 


Fic. 4. 


photograph. Magnification is 20 times. 
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A further verification of the closeness of the approxi. 
mation to the ideally porous plate may be seen from 
Fig. 5, which shows the variation of velocity v9 with 
the pressure drop across the porous plate as determined 
by two different methods: by a calibrated hot-wire 
located 0.020 in. from the porous plate with U = 9 
and by the ratio Q/A, where Q is the volume of injection 
air flowing per unit time and A is the total poroys 
plate area. It will be seen that the agreement between 
the two velocities is good even at this small distance 
from the wall. Also indicated by Fig. 5 is the clearly 
viscous pressure drop across the plate for low rates of 
flow. 

From Fig. 4 it may also be seen that the surface of 
the plate is comparatively rough. This roughness 
is greater than that associated with the usual boundary- 
layer investigations. However, with porous plates 
made of sintered metals, such surface roughness must 
be expected, since the surface cannot be smoothed 
without damaging the pores and, therefore, the perme- 
ability of the material. The aerodynamic effect of 
this roughness can be established only experimentally, 

In order to eliminate the effect of contamination 
from the sides of the porous plate parallel to the free 
stream, all the measurements with injection were 
made within the most upstream 24-in. section of the 
plate and along its centerline. Thus, even if turbulence 
is assumed, because of the effect of the injection, to 
spread from these edges toward the centerline of the 
plate at a greater angle than that shown by Charters,” 
no transverse contamination could be expected. 








Hot-wire probe and associated electric equipment 
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A hot-wire anemometer was used to make the velocity 
measurements. The wires, which were made of plati- 
num either 0.001 or 0.002 in. in diameter, were mounted 
in the probe shown in Fig. 6. A constant resistance 
type of circuit was used in order that the effects of the 
wall on the hot-wire response could be calibrated at 
zero flow velocity and could be applied independently 
of the flow velocity. The circuit used is shown in 
Fig. 7. 

In carrying out the tests, essentially four quantities 
had to be measured or chosen. These variables are 
the distance from the plate y, the velocity parallel to 
the plate , the injection velocity v, and the distance 
The techniques used to measure 


along the plate x. 
The distance 


or to determine each may be of interest. 
of the wire from the wall was determined by moving 
the wire mounted in the probe slowly toward the 
plate until a simple d.c. circuit connecting the microm- 
eter barrel and the porous plate was closed by the 
contact of the wire on the plate. The micrometer dial 
was then set to zero and the wire backed off the desired 
distance, which then was considered to be the value of y. 

In determining the values of the velocity u, the indi- 
cations of the hot-wire, which are really proportional 
to the resultant velocity, were interpreted as being 
due to the u velocity component only. This is justified 
by the smallness of the v velocity components even at 
the highest rates of injection considered here. The 
technique used to obtain the u data is suggested by a 
consideration of King’s formula for heat transfer from 
a cylinder at low Mach Numbers. The heat balance 
between the electrical energy input and the heat 
convected to the fluid and inducted to the wall from 
the wire can be written as 


FE? = R(T, — T)(AVu + B) (1) 
where 
R = resistance of the wire 
E = voltage drop across the wire 
| je = temperature of the wire 
7; = temperature of the fluid stream and of the 
wall 
1,B = constants depending on the physical proper- 
ties of the fluid medium and of the wire 
and on the distance from the wall 
For a constant resistance (constant temperature) 
hot-wire circuit, Eq. (1) becomes 
E? = A'vV/u + B’ (2) 
or 
u = [(E? — B’)/A’}? (3) 


To determine the constants A’ and B’ for a par- 
ticular wire, the following technique was used: With 
the turbulence channel sealed at both ends and the air 
therein allowed to become quiescent, a heating current 


giving a value of 7, 7, of approximately 50° to 
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H.W. 

















Re 


Hot-wire anemometer, constant resistance circuit 


Fic. 7. 


150°F. above ambient was passed through the wire. 
At various values of y, the voltage drop required to 
temperature was re- 

Fkyo(y). Then from 


wire at constant 


ky = 


maintain the 


corded and denoted by 


Eq. (2) since u = 0, 
B’ = EF,’ (4) 


If the actual magnitude of « was of interest in a par- 
ticular test, the value A’ was determined by recording 
the values of E in a stream of known velocity. By 
plotting the results in the form of E? versus ~/u, A’ 
was easily determined from the slope of the resulting 
However, for most tests, only the ratio 
therefore, 


straight line. 
u/U had to be determined by the hot-wire; 
the constant A’ did not have to be explicitly deter- 
mined. By applying Eqs. 3 and 4 at two points in the 
flow with velocities of u and U’, one finds 


u/U = [(E? — E,*)/(E;? — E?)]? (5) 


where £; is the voltage drop across the wire in the free 
stream of known velocity U. Because of the simplicity 
of applying Eq. (5), it was possible to calibrate the 
wire immediately prior to each test so that daily 
effects, such as wire aging and fluid temperature 
changes, were minimized. 

In the determination of the injection velocity v, 
the pressure drop across the porous plate was measured. 
Fig. 5 was used to determine and to fix the value of 7%. 

Because of the boundary layer that develops on the 
inner wall of the inlet nozzle, the effective distance 
along the plate could not be measured from the be- 
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X= x, ea cca Reynolds Number corresponding to this sudden 
: change was called the transition Reynolds Number 

Yo . . . - - . . . 

Ae are rang ar are oo a ar a oe a e A clear indication of the onset of transition was ob. 
f- =. ae ae tained in this way. The extent of the transition 
Pi eT Late =—s region downstream of the point corresponding to thes 
a — values of R,,, can be determined only by an invest; 


“POROUS PLATE \ 
INJECTION BELL 
INLET BELL 


Fic. 8. Inlet conditions 


The effect of the inlet 
bell with its accelerating pressure gradient was con- 


ginning of the porous section. 


sidered to be equivalent to an impermeable flat plate 
section of length x, as shown in Fig. 8. The length 
With the smallest in- 
jection rate that could be realized, without reversing 


x; was determined as follows: 


the supercharger and using it as a suction pump, the 
impermeable plate conditions were closely approxi- 
mated ( < 1). the 
leading edge of the porous plate was measured under 


A velocity profile close to 


these conditions. The value of x. in the Blasius variable 
n = (y x)VR, = yV U/xv was chosen so that the 
slope of the velocity profile at the wall (Ow/Oy)) was 
equal to the theoretical value of Blasius. The value 
of x; was then simply determined by subtracting the 
location of the wire with respect to the leading edge of 
the porous plate from the value of x so determined. 
This same value of x; (4.57 in.) was then used for all 
subsequent tests. 

The validity of using this effective inlet length is 
established by Fig. 9, where, for several Reynolds 
Numbers and several values of x;/x (including the one 
used to determine x,), the velocity profiles are com- 
pared with the classical Blasius profile. Good agree- 
ment is seen even for large values of » and for value 
of x;/x close to unity. Thus, during the tests with 
injection and suction, the boundary layer was con- 
sidered to be that over a partially porous plate with an 
impermeable section of length x, and a permeable 
section of length (x — x;,). All values of x used here 
are therefore measured from the leading edge of an 
imaginary plate x; ahead of the leading edge of the 
porous plate. 


Two types of data were collected in these tests. 


The effect of injection or suction on the transition’ 


Reynolds Number R,,, for two values of x;/x was 
established, and suction and injection velocity profiles 
were determined. 
with the theoretical laminar boundary-layer stability 


The transition data were compared 


limits, while the velocity profiles were compared with 
The Reynolds 
Number was determined by placing the wire close to 
the wall (y = 0.005 in. — 0.010 in.). With a fixed 
injection velocity, the free-stream velocity was gradu- 


theoretical calculations. transition 


ally increased. A plot of u/U for the given value of 


y increased smoothly until some value of L’ was reached 
The 


at which «/U’ would change disproportionately. 


gation of the turbulent boundary layer, and it is not 
shown by the tests reported here. 
TEST RESULTS AND COMPARISON WITH THEORY 


Fig. 10 shows the effect of suction and injection on 
the transition Reynolds Number R,,, for two values of 


v,/x. Although most of the transition data has been 
collected for x;/x 0.381, three values of R,,, corre 
sponding to x;/x = 0.190 have been measured. Since 
these fall on the same curve as those for x;/x 0.381, 


it has been concluded that the transition curve can be 
considered to be the transition curve for the completely 
0). 

Several comments may be made with respect to the 


porous plate (x;/x = 


comparison between the transition Reynolds Numbers 
shown on Fig. 10 and the results of laminar boundary 
that 
Moreover, the effect 


layer stability theory. It is seen with suction 


the boundary layer is stabilized. 
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of injection is to reduce the stability and to cause 
earlier transition. It will be noted that the stability 
theory appears to predict accurately the suction ratio 
for which the laminar boundary layer is indefinitely 
stable, since the results of reference 16 give —11.8 X 
10-* for this value and Fig. 10 indicates a rapid rise 
in R,,, in the neighborhood of #/U = —10 X 10~°. 
Moreover, it indicates qualitatively the effect of suction 
and injection on boundary-layer stability. 

Figs. lla and Ilb show the agreement for two 
different values of — between the measured velocity 
profiles and those of Yuan.’ The theoretical results 
have been corrected for the partial injection in ac- 
cordance with Eq. 17 of reference 7. It will be noted 
that the agreement is good and, therefore, that Yuan's 
analysis of the effect of injection on skin friction is 
experimentally substantiated. Thus, it can be ex- 
pected that heat-transfer data based on the analysis 
of reference 7 will be accurate for at least low rates of 
heat transfer where the assumption of constant gas 
properties is satisfactory. Moreover, the range of 
validity of these data can probably be increased if the 
fluid properties (p, u) are evaluated at the wall tempera- 


ture. The effect of heat transfer on the stability of the 
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laminar layer can be estimated by application of the 
method given in reference 18. 


Fig. 12 shows the agreement for one value of & 
corresponding to suction between the measured and 
the theoretical profiles of Yuan. It is pointed out that 
the same agreement would have been achieved if the 
comparison had been made with the profiles calculated 


according to Romano.'® 
CONCLUSIONS 


From the experimental investigation reported here 


the following conclusions may be drawn 


(1) Porous plates that satisfactorily approximate 
the ideal or mathematical porous plate—that is, one 
providing continuous normal velocities at the wall 
The microstructure of such plates 
rates of 


can be realized. 


will probably be complex, and, at low 
flow, the pressure drop will be proportional to the veloc 
ity. 


(2) Regions of laminar flow with injection can 
exist despite the inherent roughness and the finite 
number of pores of real porous plates. The laminar 
boundary-layer stability theory predicts qualitatively 
the effect of suction and injection on the transition 


Reynolds Number. 


(3) The isothermal laminar boundary layer on a 
flat plate can be stabilized to an indefinitely high 
suction corresponding to ap- 
This is in 
laminar 


Reynolds Number if 
proximately (—%/U) = 10 
good agreement with the 
boundary-layer stability theory of Ulrich,'® of Libby, 
Lew, and Romano,'* and of Romano." 


' is applied. 
predictions of 


(4) The laminar boundary-layer analysis of Yuan’ 
and Romano" for the flow over partially porous flat 
plates can be agree- 
ment between theory and experiment is found. The 
Yuan therefore, be con 


verified experimentally. Good 


can, 


data of 


heat-transfer 
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Fic. 12. Comparison of theoretical and experimental velocity 
profiles with suction for — = 0.025 
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sidered to have been verified, at least, for low rates of 


heat transfer. 
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Effect of Wind-Tunnel Nozzle on Steady-Flow M2 (*)(* ) - if ; (*)] 5/p ‘- 


7 Ape : ~ 
Nonuniformities Us AIG, 3 A, 1 + (5/p) 
which reduces to Prandtl’s result for 6/p = 0. If the density 
Stanley Corrsin s es t : ae . 
irregularity is due to temperature irregularity in a_ perfect 
Associate Professor, Aeronautics Department, The Johns Hopkins gas—the usual case—then, 


University, Baltimore 


November 12, 1951 wa - “)(*) 4 l = (*)'] : (4) 
U2 A, U; 2 A, 7 


INTRODUCTION where J is undisturbed absolute temperature and # is the 
irregularity. There is no restriction to })< T. 

A numerical example is rather startling: Suppose A2/A; = 
0.05, m/U, = 0.03, and 8/T = 0.005. Then Eq. (4) gives 
uo/U, = 0.000075 + 0.0025, and the 1.5°C. temperature irregu- 
larity causes 33 times as much trouble as the 3 per cent upstream 


pas PRANDTL’S DISCUSSION on the attainment of uniform 
flow in a low-speed, constant density, wind tunnel,! de- 
signers have been acutely conscious of the benefit of high con- 
traction ratio as a means of reducing the dimensionless measure 
of steady-flow irregularity in the working section. 

Two simple generalizations are of practical interest: (a) the 
> ~ ys 2 densi y ~ » er: a i > ‘ ities on ‘ 7 
effect of uy stre am de isity ( temperature ) irregularities o1 Cospansemer Fuow 
working section velocity irregularities in a low-speed tunnel; 
and (b) extension to compressible flow of Prandtl’s rough estimate 
for the effect of upstream irregularities, now in both mean velocity 
and stagnation temperature. CpTo = CpZ1 + (1/2)U i? = cpT2 + (1 

The analyses are restricted to small irregularities in velocity, ad Ba 2 '" 

. f ‘ ‘ . * . ‘ tae Cy To + do) = cf Ti + A) + (1/2)(U1 + u,)* = 
with resulting linearization. Velocity irregularities are con- be : , 
: ‘ ° ° cp T2 +- J) + (1 2)( U2 +- U)* (6) 
sidered only in the mean-flow direction. . 


velocity irregularity! 


The energy equations for unperturbed and perturbed stream- 
lines are, with cp assumed constant, 


2)U2? (5) 


The simplifying assumptions are (1) inviscid (isentropic) flow where 
in the nozzle and (2) constant static pressure across the duct a = stagnation temperature 
both before and after the nozzle. Io = stagnation temperature irregularity 
Low-SPEED FLow (M —~ 0) = ioe peaneandomeomninay ie 
3, 3. = flow temperature irregularities 
We write two Bernoulli equations between Section No. 1, well Cy = specific heat at constant pressure 


ahead of the contraction, in the ‘“‘pressure box,’’ and Section - ; : 
With the assumptions of isentropy and constant pressure 


No. 2 in the working section. Mises ‘ : 
across Section Nos. 1 and 2, 


pi + (1/2)pU2? = p2 + (1/2)p U2? (1) pr a) (y-1) (7 + sy _ 
2)? (2) p2 7 T2 8 T2 + de. 


l : l 
Aat+o(ot+6(Ui 4 m1)? = p2 + 5 (ep + 6)(U2 + ua)? 
i i: where 7 is the ratio of specific heat. Thus, 


where a - - 
de T2 = d; Ti (7) 
p = static pressure : don : , 
' . Linearization of the difference between the right-hand 
p = unperturbed density i a - ’ : : 
. ae , equalities of Eqs. (5) and (6) for (for /Ui, u2/Us << 1) gives, 
6 = density irregularity 7 z oa : : ‘ 
, . with Eq. (7), an expression for the working section velocity 
U = unperturbed velccity : oN 
irregularity 


u = velocity irregularity 


" : . 2 1\? eo T>» 
Eq. (1) applied to an unperturbed streamline, Eq. (2) to a a ({ ) (* ) re —, (: fk (8) 
U2 Us} \ Ui U;? T\ 


perturbed one. 
We are interested in finding u2/U2 as a function of m:/U,, A convenient final form is in terms of Mach Numbers (M,, 


5/p and the contraction area ratio, A;/A2 = U2/U;. With M:) and stagnation temperature. For this we note that 
m/U;, u2/U, << 1, the difference between Eqs. (1) and (2) a ae aes 
' . (U,/Us)? = (Mi/Mz)XT2/T)) (9) 


yields the linearized result 
From Eq. (5) it is observed that 


* Acknowledgement is made to R. Smelt, of Sverdrup and Parcel, for 
his helpful criticism U2? = 2c,(T>. — Te 
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and from the (velocity perturbation) linearized 


between the left-hand equalities of Eqs. (5) and (6), 


Cpt C po ait My ) 
- : — (10) 
U.? U,? U2 U; 


Using Eqs. (9), (10), and (5a) in Eq. (8), we arrive at a simple 


form of the final result 


Gi) +2 Ge) at) ae 


which reduces to Eq. (4) as MW, Mo — 0. 
A nonlinearized result, involving only somewhat more algebra, 
is perhaps not justifiable in view of the isentropy approximation. 
With the additional restriction of 3;/7), d0/7T.< 1 (but not 
negligible), Eq. (11) can be transformed to a simple relation 


mo/ Ms. = (M,/M2)2(m/M, (12) 


where m is Mach Number irregularity. 

It is noteworthy that, for heat removal by rods or vanes in 
Section No. 1 (% <0; d& < 0), the effects in Eqs. (4) and (11) 
are additive, while for heat addition rods (#; < 0; &%& > 0) they 


are opposite. 


REFERENCE 
' Prandtl, I Herstellung Einwandfreier Luftstrome (Windkandle), Hand 


buch der Exp. Phys., Vol. IV, Part 2, Heraus. v. L. Schiller, Akad. Verlag 
M.B.H. Leipzig 1932 (Translation as N.A.C.A. T.M. No. 726.) 


On Nonplanar Surfaces of Very Low Aspect Ratio 


H. M. Voss 
Graduate Student, Aeronautical Engineering Department, 
Massachusetts Institute of Technology 


October 12, 1951 


—— MOMENTUM APPROACH tothe problem of determination of 
the air forces on a wing of very low aspect ratio was first 
suggested by Jones.' Since that time this so-called zero aspect 
ratio theory has been extended in various directions by Spreiter,? 
for the slender wing-body combination, by Lomax and Heaslet,* 
for specific plan forms, and by Miles‘ and others® in unsteady 
flow applications. 

It is the purpose of this note to point out that the cases of a 
nonplanar surface in steady flow or in arbitrary unsteady motion 
may be treated as well, as long as the deflection, and thus the 
downwash, is developable in a finite Fourier series in y/S(x) 
(see Fig. 1). Formulation of this problem by the usual momen- 
tum approach appeared unsatisfactory because of the complica 
tion of the two-dimensional virtual mass of a deformed plate, 
the deformation varying streamwise. However, the reduction 
of the linearized lifting-surface integral, requiring only that the 
span-to-chord ratio be small, 

(y — n) << (x — &) (1) 
over most of the plan form, has been indicated by Lomax and 
Sluder.® 


5, although it contains one unnecessary restriction. 
A solution of the three-dimensional Laplace equation is of the 


Details of such a reduction are considered in reference 


form, 


Z | | ¥y(é, n, t) 
4n J J [ly — 9)? + 2?) 
Rat+ Ru 
j (x — &) 
41 + 2 r Peres 
| (x — &)? + (y — n)? +27] 2f 
where ¢ is the velocity potential, Rg the wing region, Ry» the 
wake region, and y the discontinuity in the streamwise com- 


Px, ¥, 2, t) = 


i) 


>dtdn (2 


AERONAUTICAL 


difference 


SCIENCES—FEBRUARY, 1952 

















; a 
"a 2 be — 
——— 2 S(x) 26, 
x,} | | 
| 
= 


Fic. 1 


ponent of fluid velocity from top to bottom surfaces of the plate 
It should be noted that Eq. (2) holds for both steady (y = () 
in Rw») and unsteady flow. The assumption of Eq. (1) along with 
s negligible compared to (x £) in the bracketed portion of the 


kernel leads to, 


4 . ME y(£, n, t) 
DX, Y, 2, t) = : dé dy 
2m 0 « S(ég \y ae 


Here S(£) is the semispan as a function of the streamwise variable 


rn 


The kernel of the integral of Eq. (3) is seen to be a solution of 
the Laplace equation in two dimensions; the integral of y over 
the chord serves as a distribution function. Thus, this solution 
might have been constructed by assuming the validity of the 
two-dimensional solutions in the present case and superimposing 
them. In this respect, the foregoing formulation serves only as 
a justification of reduction to the quasi-two-dimensional problem, 
as suggested in reference 1 on a physical basis 

A function F(x, y, 4) may be defined, 


x 
/ VE, Y, t) dé (4 
LI 


where LE is the streamwise location of the leading edge. The 
order of integration may be interchanged and Eq. (3) integrated 


by parts to give, 


l "S(x OF y n f 
o(x, y, z,¢) = tan~! dn (5 
2r On Z 


« x) 


Eq. (5) is then differentiated with respect to s, 


Od l te OF 7 n) 
(% 9,23, 0) = dn (6 
ra) 2r © On |(y n)*? + 2? 
and the boundary condition is satisfied in the s = 0 plane in the 


usual manner of thin airfoil theory, 


Od l a3 OF dn - 
(x, y, O, t) = wx, y,t) = (7 
Oz 2r S(x) On (y n 


In this, w(x, y, ¢) is the vertical velocity on the airfoil, positive 
downward 
Eq. (7) is in a form suitable for the familiar Schwarz!! inversion, 


becoming, 


: ~2 pte “S(x)2 2 
inne f£ w(x, 9, t) V Se? — ¥ ag 
 l V S(x)? — y? S(a (y — 9) 


(8) 
The required form of the downwash is, 
\ 
w(x, n, t) = f(x, t) Zz. a, cos né 
@ = cos™! [n/S(x)] yw = cos —[y/S(x)] (9) 


The integrals are of the Glauert type and can be evaluated ex- 


plicitly, 
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OF — f(x, t) 
(x, y,!) = 
dy V S(x)? — y? 
A second integration gives, 
y 


F(x, y, t) = f(x, t) vt 44 Day + a— 
. J V S(x) J ) Sx 


Differentiation with respect to x then yields y and leads directly 
to the pressure distribution over the plan form as a function of 


xand y, 


(x, y, t) (12) 


Ap(x, y, t) = p| Uy(x, y¥,t) +° 
f ot 


The applicability of Eq. (2) to steady or unsteady flow follows 
from the assumption of Eq. (1), leading to the conclusion that 
the influence of the wake in unsteady flow is negligible. Law- 
rence’ has pointed out the fact that no low aspect ratio analog 
exists for the two-dimensional theory for high aspect ratio, and 
that, in fact, the zero aspect ratio theory of Jones corresponds 
to lifting-line theory for high aspect ratio. In that the large 
aspect ratio theory of Reissner reduces to lifting-line theory for 
steady flow, that the assumption of Eq. (1) is the opposite extreme 
of that made by Reissner, and that the result reduces to the 
Jones result for steady flow, the suggested zero aspect ratio theory 
for unsteady flow becomes the low aspect ratio analog to the 
Reissner theory for large aspect ratio. 

In connection with the use of zero aspect ratio theory in un 
steady flow, Miles has made the well-taken point that the Kutta 
condition remains unsatisfied at the trailing edge in steady flow, 
This 


point was discussed in reference 5, and Halfman and Ashley* 


so that extension to unsteady subsonic flow is dubious 


have reported on such an application. 
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Note About Secondary Flow in Cascades 


Hans P. Eichenberger 
Gas Turbine Laboratory, Massachusetts Institute of Technology 
October 99, 1951 


’ I ‘HE NOTE OF STEPHENSON! raises an interesting question about 


the trailing vorticity behind a cascade ot airfoils in shear 
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\ 
) 2aoy + a,S(x) cos 2y — 2S(x) sin y¥ : o a,, sin nv (10) 
n=2 


(11) 





flow. Assuming that the number of blades is infinite, Stephen- 
son adds the secondary vorticity to the trailing vorticity to ob- 
tain the resulting effect. The question arises whether this cor- 
responds to physical reality. 

The change of vorticity for inviscid flow is given by the Helm- 
holz equation (reference 2, page 205) (q’'V)2 = (2V)q, where q 
is the velocity vector and =V-g. If u,v, and w are the com- 
ponents of g and &, , and ¢ are the components of & in direction 
of r, y, z, respectively, this equation can be written: 


v n 
uty + - bo + wk = iu; + Uy + CU, 

r r 

v vé ” nu 
un, + Ne + wn, + = tu; + 9 + $i. + 

r r r r 

ve n 
us, + -S$o + wy, = tw, + - we + fu, 

r r 


They show that we have to know the streamlines in order to 
compute the variation of 2 which does not stay constant, in general 
Assuming, as Squire and Winter® did, a basic free vortex flow 
and u, v, w small as compared to v% = c/r, we can linearize the 


above equations and obtain immediately: 


(1) & = O, — stays constant 
(2) » = —2& as obtained by Squire and Winter* 
(3) f~ = O, ¢ stays constant 


Eq. (2) is in disagreement with Eq. (1) of Stephenson’s article 
because he assumed the vorticity constant. Fig. 1 illustrates 
qualitatively the variation of the vorticity; we assume that up 
stream of the cascade the axial velocity component only has a 
peaked profile. The line AB represents then a vortex tube 
downstream of the cascade this line A’B’ is no longer parallel to 
AB, but A’ lags behind B’ since the flow (and with it the vor 
ticity) goes more slowly on the pressure side than on the suction 
side of the blade. 

Knowing the vorticity 7 in direction of the streamlines in the 
whole r-s plane at the outlet of the cascade, the problem of find 
ing the secondary velocities has a unique solution, and the stream 
lines may look as sketched in Fig. 2. From this result we may, 
if we wish, calculate directly the so-called trailing vortices that 
behave like needle bearings between two rotating masses and 
form a vortex sheet. A vortex sheet is a name for the physical 
fact that we have a discontinuity’in velocity over the plane of the 
This mathematical concept has proved to be useful in the 




















sheet 
— es 
ALLL LLLh$. 
P “ 
: _\g' __ 
AL’ — 
: z 
A\, + 
a W\. = 
VA TITTTTITITITITITITTTIT? 
A Ps 
A 4 Fig. | 











138 JOURNAL OF THE 


AERONAUTICAL 























Fig. 2 





theory of the finite wing in a potential flow, where on the pres- 
sure (suction) side the streamlines bend toward the wing tip 
(midspan) such that we have at the trailing edge a surface of dis- 
continuity. From the magnitude of this discontinuity, we can 
calculate the induced velocities everywhere because the vorticity 
is zero everywhere except on the trailing vortex sheet. But ina 
cascade where the vorticity in the stream is different from zero, 
the knowledge of the trailing vortex sheet is useless for the cal- 
culation of the induced velocities. The knowledge of the dis- 
tributed vorticity is necessary and sufficient to calculate the secon- 
dary velocity. The jump of velocity at the trailing vortex sheet 
then follows immediately. To add this trailing vorticity to 
the distributed vorticity seems not to have any physical mean- 


ing. 
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The Secondary Flow in a Cascade of 


On Chang’s Function for Nonstationary Flow 


John W. Miles* 


Fulbright Lecturer, 
New Zealand 


September 14, 1951 


Auckland University College, Auckland, 


I’ A RECENT REPORT on nonstationary airfoil theory,! Chang 
introduces (and recommends the calculation of) the function 


1 (°* 
C(s, M) = f exp [—78M(M — cos 6)~'] do (1) 
rJ 0 


where M is the Mach Number and £8 is twice the usual reduced 
frequency. It should be pointed out that this function may be 
expressed in terms of Schwarz’s function,? fo(A, 17), according to 


/2f)[8M?2(M?2 — 1)-1, M] 


(2) 


C(B, M) = 1 — i8M(M? — 1) 


Similarly, it is expected that his C(8, 4, 6;) can be related to his 
“incomplete Bessel function’’ Jo(z, 6;). 

To demonstrate Eq. (2), take the Laplace transform of Eq. (1) 
with respect to 8, whence 


* On leave from the University of California, Los Angeles 
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F(s, M) = f e~ S®C(B, M) dB (3a 
J0 
1 {x ; 
a J [s + 1M(M — cos 6)~')~!d@ 3b 
xr J 


= s-! —14s-1M(M? — 1)7~ /*{[s + iM? x 
(M? — 1)-*]? + M%(M2 — 1)-2}-"2 
Taking the inverse transform of Eq. (3c), we find 


1 


C(B, M) = 1 — iM(M? — 1) * a 


op 
f exp [—i2/?(M? — 1)~!7] Jo[M(M? — 1) r]dr (4 


Introducing « = 7/8 as the variable of integration and noting the 
definition 


be 
fo(A, M) -{ e~*" Jo(x\u/M) du (5 
0 
Eq. (4) reduces to Eq. (2), Q.E.D. 
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On the Stability of the Boundary Layer with 
Respect to Disturbances of Large Wave 
Velocity 


C. &. in 

Associate Professor of Mathematics, Massachusetts Institute of 
Technology 

October 26, 1951 


| geetnapee THERE HAS BEEN considerable interest in the sta- 
bility of the boundary layer with a supersonic free stream."* 
In particular, the stabilizing effect of cooling at high Mach Num- 


bers has been the center of discussion. Since ‘subsonic’ dis- 
turbances only are considered proper to the stability pioblem, 


the wave speed c of the disturbance must satisfy the relation 
c/U>1 — (1/M) 


where U is the free-stream velocity and .)/ is the free-stream Mach 
Number. It is clear that for reasonably high Mach Numbers, 
the value of c/U is large. 

It has been argued that the methods of evaluation of the in- 
tegrals associated with the ‘‘inviscid”’ solutions are inadequate 
for large values of the wave speed. This is indeed the case, 
since they were developed analogous to those used for an incom- 
pressible fluid,? where the wave speed is usually small. This 
fact has been used explicitly in the original development of the 
method, and the procedure has been adopted in the compres- 
sible case‘ without essential modification. For wider applica- 
tions, this restriction must be removed as has now been gener- 
ally recognized. 

There is, however, a much more basic restriction inherent in 
the theory developed by Lees and the present writer,> on which 
In Eq. (72) of reference 


5, the velocity and density distributions are expressed as power 


the calculations of Lees* were based. 


series, and only the linear terms are retained by an examination 
of the order of magnitude. This step restricts the ‘‘viscous” 
solutions to the neighborhood of the critical layer and makes it 
unreliable, in the case of moderately large values of c/U, for the 
purpose of satisfying the boundary conditions at the wall. One 
might attempt to improve the accuracy by going to higher terms 
in Eq. (72), but the amount of labor required would be excessive 
if the density distribution across the boundary layer requires 
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ral terms of the series for its accurate representation. Basi- 
cally, the method associated with the use of Eq. (72) is also not 
satisfactory. It leads to the conditions [Eq. (95)] which, though 
adequate for the purpose at hand, have been recently shown by 


seve 


Cathleen S. Morawetz' to be inexact. 

Since all the present calculations for large c/U are based on a 
procedure beyond the validity of the existing theory, attempts 
are now being made to revise the theory considerably to remove 
the restrictions mentioned above. 

It should perhaps be pointed out in passing that all the devel- 
opments in Section II of reference 5 are not influenced by the 


restrictions discussed above. 
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A Method of Solving the Linear Potential 
Equation for Axially Symmetric Flow 


R. Probstein and J. V. Charyk 


Department of Aeronautical Engineering, Princeton University 
October 29, 1951 


— LINEARIZED POTENTIAL EQUATION for supersonic flow 
over a body of revolution whose axis is in the direction of the 


free stream velocity is, 
0’ 0° 1 Op 


) = () (1) 
Ox? or? ror 


(1 — M;? 


where ¢ is the disturbance potential that is a function only of 
The general, solution of Eq. (1) can be shown to be! 


1 fm = f(&) dé 
(2) 
tr JO V(x — §)? — ar? 


xand r. 


oe = 


When the above equations are integrated, we obtain, 
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1 and /(£) is the unknown source distribu- 
Both the position and variation 


where a = V M,;? — 
tion function along the x-axis. 
of the source strength, which is expressed by /(£), are determined 
by the body shape. 

At present, the general method used for solving this integral 
equation for the flow about an arbitrary pointed body of revolu- 
tion is to break up the body into a series of conical flow sections 
and then superpose the readily soluble linearized solution for 
conical flow. For most body shapes, this procedure is extremely 
long, and, in many instances, the only practical method of solu- 
tion would be to resort to computing machines. An alternate 
approach, which does not appear to have been used previously 
and which yields an apparently excellent approximation, is pro 


posed. 
If one recognizes that the source distribution function for linear 
ized conical flow is given by f(t) = Aé, then it would be reason 


able to assume that the source distribution in the general case 
could be represented by a polynomial, such that, 
p 


f(t) = DO Ant (3) 


n= 1 


~ 


Here the coefficients 4, are determined by the boundary condi- 
tion that the flow must everywhere be tangential to the body 
that is, that 
m = (dr/dx) = (0¢/0r)/|U; + (0¢/0x)] (4) 
To illustrate the simplicity of this method, let us calculate the 
pressure distribution over a body of revolution at zero angle of 


attack. Since Eq. (3) is rapidly convergent for reasonable body 
shapes, we shall chose p = 3. If we now make the transforma 
tion 

— = x — arcosh B (5) 


then the disturbance potential is given by, 


1 cosh ~! (x/ar) 
q=- +B A,(x — ar cosh B)" | ds (6) 
tr JO 
_n 1 
and it follows that, 
3 
Op l cosh ~! (x/ar) 
= — oe nA,(x — ar cosh 8)"~! |] dB 
Ox tr JO a 
(7) 
and 
3 
a Pcosh ~! (x/ar) : 
n= = | > ® nA,(x — ar cosh B)""! 
or tr JO 2 1 


cosh B dp (8 


fe.) l 3 x / (n — 1) (n — 2) (ar)? | x | s 
—._ — A, —1 eosh-! - a n—-2 VWV/x2 ee 1 + cosh™! 
oe . 2» n yx"! co h on (n 1) x x (ar) = 2 i s | = one 
(9) 
and 
op ‘ E lf x \2 (n l)x*-*ar | x | s\ x 
=— > nA, \ x" @/ 1 — = | — 1+ cosh" + 
or Wino s \ ar 4 n ar \ ar ar 
(n — 1)(n — 2) J(ar)? I/ x \2 x \? { 
. | — 1 + 2 ( (10) 
2 | 3 ar ar 


The evaluation of the integrals in Eqs. (7) and (8) simply re- 
duces through integration by parts to the evaluation of integrals 


of the hyperbolic cosine to various powers. The following re- 


currence relation permits the result to be readily written down. 


m— | . 
sinh B cosh”! 8 + JS cosh” *3dB 
m m 


S cosh” 6B dB = 
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X= DISTANCE ALONG BODY AXIS 


Fic. 1. Comparison of calculated pressure coefficients for a 
parabolic body of revolution described by r = (x/100) (20 — x) 
at a free-stream Mach Number = 2.6. 


Eqs. (9) and (10), since they are limited to n = 3, involve only 
the results of the integration of terms of the order cosh? 6 and 
cosh® £8, respectively. It will be noted later that the effect of 
angle of attack only demands the evaluation of integrals of cosh* 
8 and cosh* 8 for the case n = 3—that is, the only complication 
is the raising of the degree of the cosh £6 in the integrals by one. 
The complexity of other methods when angle of attack is intro- 
duced is appreciably increased. 

The boundary condition (4) applied at three points on the sur 
face of the body can be expressed in the fori: 

3 


7 Aa, (x, ar) 


n= 1 
My = (11) 


3 
U;, + > A,b,(x, ar) 

n= 1 
where, the a,’s and b,’s represent the functions of x and ar asso- 
ciated with the corresponding A,,’s. The A,,’s are now readily 
calculated from the following set of three equations. 


3 3 3 
U; > m = >» Z A,(ant — mibyr) (12) 
l= 1 l= 1 n= | 


Here / represents the three points along the body at which the 
slope istaken. Since, 


3 

Od 
= p Anb, (x, ar) (13) 

Ox n 1 

we can calculate the pressure coefficient immediately from, 
a 3 
o=> Zz. A ,b,(x, ar) (14) 

l ln 1 


It should be pointed out that the problem is now solved once 
and for all for any body of revolution for which p = 3 yields a 
sufficiently convergent series in Eq. (3). It might be necessary 
to take a larger value of p for bodies in which the describing 
curve is of extremely high order; however, this does not modify 
the method of approach. If, for example, p were chosen to be 
4, this would only necessitate integrating Eqs. (7) and (8) one 
more time, yielding _a fourth-order determinant for Eq. (12) in- 
stead of a third-order one. 

It should also be noted that the same method can be applied 
to bodies of revolution at an angle of yaw as well. One could 
then obtain simple series expressions for the stability derivatives 
At present, this is done by resorting to the von Karman-Moore 
approximation of extremely thin bodies* because of the cumber- 
someness in the handling of the method of conical flows first 
described. 

To show the accuracy of the method, calculations of the pres 
sure coefficient were made for a parabolic body of revolution 
whose shape was given by r = (x/100) (20 — x) at a free-stream 


Mach Number equal to 2.6. The results obtained are shown in 
Fig. 1 and are compared with the method of calculation by Conical 
flows, as well as with the von Karman-Moore approximation. 
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On Slender Wing-Body Theory 


John W. Miles* 


Fullbright Lecturer, Auckland University College, Auckland, 
New Zealand 


November 14, 1951 


TT FREQUENT REFERENCE to, and use of, Spreiter’s analysis 

of the slender wing-body problem! (e. g., numerous N.A.C.A 
reports) implies that Ward’s more complete analysis of the same 
problem* * has not been fully appreciated in this country. In 
particular, Spreiter’s results are incorrect for a winged body of 
varying cross section because of his failure to include properly 
the logarithmic term in the complex potential and the second- 
order terms in Bernoulli's equation A detailed criticism has 
been communicated privately to Spreiter (1950) 

The errors in Spreiter’s analysis explain the disagreement 
between his results and those of Ward for a winged body of 
revolution, the latter obtaining a lift coefficient of similar form 
for all configurations for which the body cross section either 
remains constant or terminates aft of the maximum wing span 
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A Note on the Design of Supersonic Wind- 
Tunnel Nozzies 


R. E. Meyer 
Institute for Fluid Dynamics and Applied Mathematics, University 
of Maryland, College Park, Md 


November 1, 1951 


OR THE DESIGN OF SYMMETRICAL NOZZLES for two-dimensional 

supersonic flow, it is customary to divide the field of flow 
into the subsonic and the supersonic part and to treat the latter 
by the numerical method of characteristics and the former by 
some other approximate method. It is then convenient not to 
extend the characteristics computation to the sonic line but to 
make the transition at some other isobar, where the velocity 1s 
slightly supersonic. One of the main purposes of the calculation 
for the subsonic part is then to find the position of, and the 
distribution of the stream direction on, this isobar 

Since test sections with a uniform distribution of velocity have 
been found desirable for some purposes, interest has been shown 


recently in obtaining a good approximation to the initial condi- 
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tions for the characteristics computation. It may be worth 


noting that the maximum value of the stream direction on such 


an isobar depends only on the velocity magnitude. This pro- 


vides a useful check on the accuracy of the approximation used. 

In a symmetrical nozzle, the two Mach lines through the sonic 
point on the axis are branch lines.! That is, the sense of varia- 
tion on an isobar of the stream direction is reversed where the 
isobar crosses either of these Mach lines. Conversely, it can 
be shown? that such a reversal can occur on an isobar only across 
a branch line, if the wall curvature is continuous, upstream of 
the isobar 

The characteristic equations for the two-dimensional steady 
isentropic irrotational supersonic flow of a perfect gas are® 


dy/dx = tan(@ — yw); @+¢ = const. = a 


and 
dy/dx = tan(6+ 4); @—t = const. = B 


where » is the Mach angle, @ is the stream direction (measured 
counterclockwise from the downstream direction on the axis, 


which is the direction of x increasing), and 


t=pn+Acot '(Atan pw) — (1/2)r 


with 
A= Vi(v + LDA — 1)) 
y being the ratio of the specific heats. Bernoulli’s equation 


relates » to the velocity magnitude, qg, so ¢ is a function of g only. 


At the sonic point on the axis, @ = t = 0, so the branch lines are 
the Mach lines a = 0 and B = 0 and the maximum value of 
@ on the isobar gq = qo is equal to (qo). 
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Parameters Predicting Failure of Partial 
Diagonal Tension Bcams 


Bertram Klein 
Structures Methods Group, Bell Aircraft Corp., Buffalo, N.Y. 
October 4, 1951 


— DIAGONAL TENSION field beams can fail in various ways. 
Two distinct modes of failure of particular importance are 
those due to failure of the web and forced crippling of the vertical 
stiffeners (uprights). For these types of failure, rigorous theoret- 
ical analyses apparently do not exist. Much testing has been 
performed by which semiempirical analyses were evolved and 
criteria established for predicting the allowable loads and stresses 
connected with these failures. ! Much scatter exists, however, 
when the test data are plotted according to those criteria for 
forced crippling. The author presents what he believes to be 
more rational parameters that reduce the scatter in this case. 
Also included are parameters that seem to give best results for 


predicting web failure 


(1) Wes FAILURE 


According to reference 1, the main parameters are: A, /dt, 


é 


Ter, /Tult., ANd Tur, /or,. (For the meaning of these and following 
1 and 2.) 
>» 


pear in the methods of reference 2 


symbols, see Figs The first two parameters also ap- 


They are connected with 
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Fic. 1. Web failure in partial diagonal tension field beams 


the angle and degree of diagonal tension developed in the web. 
To account tor differences in web and upright material, the factor 
E,/Ew is incorporated in the first parameter. Experimental 
data taken trom references 2 to 5 are plotted in Fig. 1 
son the data of reference 5 have been plotted separately is due to 


The rea- 


the fact that it was not known whether single or double uprights 
were used in the tests. 
be kept in mind that other factors that have to be considered be- 


Although there is little scatter, it should 
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EGREE OF DIAGONAL TENSION 


ANGLE OF DIAGONA) TENSION 


he t 








e tw 
| REDUCED) mM FOR PLATE BUCKLING 
FOR ULTIMATE TIEFNER TRE 
PRIGHT LENGTH MEASURED BETWEEN 
iGnT 7 FL ANGE RIVET 
¥ d 
7 
tw” fal 30 
1 x [357 
395 
Qt 
Kot + 
Joy Lok 4 } 
Fic. 2. Forced crippling of single, bent-up uprights in partial 


diagonal tension field beams. 
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fore the picture is complete are: rivet spacing, flange flexibility, 


and type of web to flange attachment. 


(2) Forcep CRIPPLING OF UPRIGHT 


Only single, bent-up uprights are considered here. From the 
work of reference 8, it appears that two factors to be considered 
are the distance between centerline of rivet and heel of upright 
and the stiffener to web rivet spacing. Since these quantities 
were not known for the test data used, the total length of the at- 
tached flange, b, was used. The effect of the rivet spacing was 
neglected. Based on the results of reference 8, one parameter is 
taken as (b/t)rz(t/d)wk tan a, the factor k tan a accounting for 
the wrinkled state of the sheet. A second parameter is h,/d, 
the aspect ratio of the panel, which controls the movability of the 
wrinkled sheet. 

The experimental data of references 2 and 4 are plotted in Fig. 
2. All specimens had bent-up uprights. It should be noted 
that an extruded section is more efficient than a bent-up section 
because the corner radius is eliminated. When the maximum 
stiffener compressive stress at failure, computed according to the 
methods of reference 2, was definitely in the plastic range, a re- 
duced modulus was used. 
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The data appear to follow a reasonable trend. However, More 
test data should be incorporated before it is concluded that the 
chosen parameters are correct for this type of failure, as well as 


for web failure. Such unpublished data, no doubt, exist, and j 


would be of great interest to the author to learn how these data 
plot and compare, using the suggested parameters. 
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